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Quelques remarques 


LIV 


VOLUME XLI 


sur [électromécanique 


quasi-stationnaire des milieux continus 


RESUMO 


O artigo mostra que, em electromecânica quasi- 
-estacionária, é preferivel adoptar, logo à partida, a 
relatividade qgalileana, À razão é que uma teoria mais 
geral a partir da qual se pudesse, de maneira consis- 
tente, deduzir a electromecânica quasi-estacionária, só 
poderá existir a um nivel demasiado complexo para as 
aplicações. 

Estas conclusões são obtidas a partir de uma nova 
formulação da mecânica dos fluidos incompressiveis, 
análoga à teoria do electromagnetismo, que permite um 
tratamento paralelo das duas disciplinas, 


1 — INTRODUCTION 


Les notions de |électromécanique des milieux con- 
tinus utilisées en conversion d'energia (machines élec- 
triques, MHD, pompes électromagnétiques, etc.) se pla- 
cent à un niveau oú les phénoménes de propagation sont 
ignorês, parce-qu'ils jouent en rôle négligeable, 

Par analogie avec la désignation adoptée en élec- 
tromagnétisme, nous parlerons d électromécanique quasi- 
-Stationnaire. 

Cette branche de l'électromécanique, qui doit se 
placer au niveau de la physique dite prê-relativiste ou 
galilêenne pose un certain nombre de problêmes en ce 
qui concerne la cohérence des modéles mathêmatiques 
utilisês. 

En effet on a d'abord mis ensemble les équations 
de la mécanique de Newton (qui sont du niveau galiléen) 
et les équations de I'électromagnétisme de Maxwell (qui 
sont du niveau relativiste), ce qui constitue une incohé- 
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M. 5. GARRIDO 
Centro de Automática da Universidade Técnica de Lisboa 


Faculté des Sciences Appliquées de Louvain-la-Neuve 


SUMMARY 


kt is shown that, in formulating quasistatic electro- 
mechanices, galilean relativity is preferably used from the 
start. This is because a more general theory from which 
electromechanics would have to be consistently derived 
could only exist at a level of excessive complexy in 
comparison with the requirements of applications. 

This conclusion is reached by setting up a new 
formulation, similar to the electromagnetic theory, for 
incompressible fluid mechanics, This allows a parallel 
treatment of both theories. 


rence. Certes, on éliminait normalement le courant de 
déplacement dans les conducteurs, mais les équations 
électromagnétiques approchées ainsi obtenues nétaient 
quand même pas compatibles avec la transformation de 
Galilée ni donc avec la mécanique de Newton. 

C'est pourquoi un certain effort a été fait ensuite 
dans le but de formuler un électromagnétisme quasi- 
-stationnaire vraiment galiléen, 

On peut trouver l'essentiel de cette formulation de 
l'électromagnétisme dans les ouvrages didactiques [1], [2] 
et dans les travaux de recherche [3], [4]. 

Malheureusement, un certain nombre de nouvelles 
difficultés vont apparaitre aprês cette reformulation de 
[électromagnétisme quasi-stationnaire, 

En effet, si on fait une étude critique de la mécanique 
newtonienne à la lumiêre des nouveaux résultats obtenus 
pour l'êlectromagnétisme quasi-stationnaire, on est três 
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vite conduit à la conclusion que cette mécanique ne vas 
pas sans problêmes dans le cadre galiléen, certaines 
difficultés apparaissant au niveau des lois de transfor- 
mation des grandeurs mécaniques, 

Il en résulte qu'on est arrivé à se retrouver aujourd' 
hui avec le même type de problêmes qu'il y a quelques 
années, mais maintenant avec un électromagnétisme qui 
se prête bien à une formulation galiléenne et une méca- 
nique qui, finalement, s'y intêgre assez difficilement, 

Ce três curieux renversement de situation est claire- 
mant montré dans la reférence [3], à propos des formules 
de transformation de l'impulsion et de l'energie. 

Dans ce travail, nous verrons que les difficultés ne 
se limitent pas aux formules de transformation de [impul- 
sion et de |'energie, mais qu'elles existent deja pour les 
formules de transformation des variables mécaniques 
elles-mémes. 

Ainsi, nous commencerons par développer les équa- 
tions quasi-stationnaires de la mécanique du continu, 
d'une façon originale, à notre connaissance, que les 
rends directement comparables aux equations quasi-sta- 
tionnaires de ['électromagnétisme. 

Ensuite, nous étudierons les formules galiléennes de 
transformation des variables électriques et mécaniques et 
nous identifierons les difficultês qui apparaissent au 
niveau de la mécanique, 

Finalement, la nature de ces difficultés sera analysée 
et une solution sera proposée et comparée à d'autres 
solutions envisageables, 


2— LES EQUATIONS QUASI-STATIONNAIRES 
DE L'ELECTROMAGNETISME ET DE LA MECANIQUE 


Les équations fondamentales de l'électromécanique 
des miligux continus sont les équations de |'électromagné- 
tisme et les équations de la mécanique du continu, 
dúment completées avec des termes provenant du pro- 
cessus d'interaction électromécanique., 

Nous allons nous occuper d'abord et essentiellement 
des équations purement électriques et des équations 
purement mécaniques et c'est seulement ensuite que nous 
introduirons le processus d'interaction électromécanique. 

Pour allêger l'exposé, tout en gardant les aspects 
fondamentaux qui nous interessent, nous allons considérer 
des milieux relativement simples, 

áinsi, pour les équations électromagnétiques, nous 
admettrons l'existence de miligux conducteurs continus, 
mais nous n'envisagerons pas de matériaux magnétiques 
ou diélectriques, 

Pour les équations mécaniques, nous allons nous 
intéresser à un milieu fluide incompressible et en régime 
isentropique. Ce milieu sera supposé continu, indéfini et 
au repos à [infini, comme on le fait souvent en méca- 
nique des fluides [5]. 

Nous allons commencer par voir que, dans le cadre 
quasi-stationnaire, les lois de ['électromagnétisme et les 
lois de la mécanique peuvent en fait, être décrites par 
des équations identiques, malgré les três grandes diffé- 
rences qui apparaissent dans les formulations tradition- 
nelles de ces deux disciplines, 
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Nous commencerons par poser une condition fonda- 
mentale, caractérisant tous les modéles galiléens: 


> (1) 


ou c répresente la vitesse de propagation, aussi bien en 
électromagnétisme qu en mécanique., 

Hemarquons que c'est cette condition qui permet de 
passer de la transformation de Lorentz à la transforma- 
tion de Galilée [3] et non pas la condition Vº /cº— 0 que 
on utilise souvent, mais qui est insuffisante, 

La condition (1), caractérisant l'espace-temps de 
Galilée, assure à elle seule le passage au niveau quasi- 
-stationnaire, et est ainsi la premiére équation de ['électro- 
mécanique quasi-stationnaire. 

Dans la formulation que nous allons proposer par 
la suite, les équations électromagnétiques et les équations 
mécaniques auront la même structure, qui est celle de 
la théorie classique de Maxwell, mais simplifite par la 
condition fondamentale (1). 

Pour ne pas devoir écrire deux fois les mêémes 
equations, nous allons adopter des lettres de |'alphabet 
manuscrit pour représenter, indifféremment, des champs 
electromagnétiques ou mécaniques et décrire en une fois, 
le formalisme., 

Dans la théorie de Maxwelle on dispose d'un champ, 
décrit par quatre vecteurs &, D «5, 2, et ayant comme 
sources deux grandeurs ;; Y. que nous supposerons ten- 
dre réguliêrement vers zéro à Vinfini: 


gay =p gy» 
D, cb, 8, j ç O (2) 


SE 
Les équations d evolution du systême sont les équa- 
tons de Maxwell 


= 
Ec o sn 
Vxã =— — (3-a) 
i dt 
—+ 
Ve 8) == (3-b) 
mm 
+ >» dd 
ET: ii 
— + 
V dA = 0 (3-d) 


jp mi 

OD = E 6 (4-a) 
B = Ho cd (4-b) 
p= 0(8+Vxars) + 9 MV (4-c) 


ou E, H, sont deux constantes universelles, electro- 
magnétiques ou mécaniques; O est un paramêtre, élec- 
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trigque ou mécanique, caractéristique du milieu matériel; 


ga est un champ appliqué, d interaction; Vest la vitesse 
du milieu matériel. 

Nous pouvons remarquer que la loi d'Ohm (4-c) est 
plus simple qu'elle ne le serait sans la condition (1) qui 
élimine les corrections relativistes du second ordre 
Vi-vi/e. 

—s> 

On peut associer au champ deux potentiels DP, get 
obtenir une autre formulation qui utilise comme variables 

=as = =) e) =p 
Va, 9, 4 et comme sources &, “2, ces grandeurs étant 


supposées tendre reguliêrement vers zéro à [infini: 


> 3995 
DM, 9, j 88,20 (5) 


Dans cette formulation les équations d'évolution sont 


> o 
V xa (6-a) 
pa. * 
da -—. 
VY + =— é (6-b) 
dt 
—» 
V.g=0 (6-c) 
af 
V.) + —— = 0 
à dt (6-d) 


ou figure la condition de Maxwell (6-c) équivalente à 
celle de Lorentz, en conséquence de la condition (1). 

Les équations (6) doivent être completées par les 
relations suivantes: 


e 4 e 
VºQ=—m) (7-a) 
g 
a o —— Dt (7-b) 
Ú 


o = —>» —» —s .. 

p=0(8+VxX81458)+9V (7-c) 

ou les équations d'Alembert apparaissent degenérées en 
équations de Poisson, grâce à la condition (1). 

Les équations précedentes ne peuvent pas determiner 

complétement l'évolution du champ, électromagnetique ou 

mécanique, parce qu'il manque encore une relation sup- 


plémentaire pour la vitesse Vº du milieu. 

C'est à ce niveau que les formulations de I'électro- 
magnétisme et de la mécanique quasi-stationnaires vont 
différer. 


ES 

Pour |'électromagnétisme la vitesse V est une donnée 
extérieure, qui doit être connue explicitement en foncition 
de I'espace et du temps. Nous avons donc besoin d'une 
équation supplémentaire de la forme 


—+ 
Ve = f(r.t) (8-a) 
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ou nous avons introduit lindice e pour signaler qu'il 
s'agit de la vitesse figurant dans des équations (1) à (7) 
concernant un champ électromagnétique. 

Pour la mécanique, par contre, la vitesse du milieu 
n'est pas une donnée extérieure, mais doit étre fournie 
par les équations du champ mécanique lui-même, C'est 
ainsi que se manifeste le caractêre non-linéaire de la 
mécanique, ou on aura au lieu de (8-a) une nouvelle 
équation mécanique, qui est 


< 
= 

I 

( 
as 


(8-b) 


ou nous avons utilisé indice m pour indiquer qu'il s agit 
d'une relation valable lorsque les équations (1) à (7) 
concernent un champ mécanique, 

Pour terminer cette introduction générale aux équa- 
tions de l['électromécanique quasi-stationnaire signalons 
encore la maniêre d'introduire le processus d'interaction 
dans les équations électriques et mécaniques que I'on 
vient de voir. 


——— 
D'abord, on fait I'hypothêse que la vitesse V. du 
milieu électromagnétique est, simplement, la vitesse du 
Rs 


milieu matériel mécanique V,, soit 


Vo = Va (9) 


ce qui élimine ['équation (8-a). 


—» 

Finalment, on introduit dans les champs &º des ter- 
mes d'interaction, comme par exemple, pour ce qui est 
des équations mécaniques, le champ de la force de 
Lorentz. 

Nous reviendrons sur la question aprês avoir analysé 
plus en détail les équations électriques et mécaniques. 

Compte tenu du formalisme quasi-stationnaire pre- 
senté pour |'électromagnétisme et pour la mécanique, 
nous pouvons préciser maintenant que la vitesse de 
propagation c introduite au début est donnée par 


c=1/Vebt (10) 


ou E, HW sont les constantes universelles, électriques 
ou mécaniques, introduites dans les équtions (4) et qui 
seront identifiées plus loin. 

La condition fondamentale c —» 2 est donc équiva- 
lente à annuler une des constantes universelles £E,, Ly —>0, 
ce qui nous permet de simplifier davantage les équations 
présentées. 

Nous allons voir en quoi consistent ces simplifica- 
tions supplémentaires des équations mécaniques et élec- 
triques et constater, en particulier, que les équations 
mécaniques sont effectivement équivalentes au équations 
classiques de la mécanique des fluides, 


2.1 — Analyse des équations quasi-stationnaires 
de I'électromagnétisme 


Pour analyser les équations (1) à (8) dans le cas 
de I'électromagnétisme nous devons, en principe, com- 
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mencer par faire correspondre les variables utilisées aux 
variables classiques de |'Geletromagnétisme. Il sagit ici 
d'une identification, puisque les équation sont justement 
celles, bien connues de |électromagnétisme. 

C'est pourquoi nous nous limiterons à changer les 


notations pour adopter ci aujourd'hui conventionnelles 
=p =p 3 =p = > —s 
en eéelectromagnétisme, & = =D, A=H g8=B, 


représentant les champs et fis inductions, électriques ou 
magnétiques; 4 Sel q =p representant la densité de 
courant et la densitê de charge eléctrique; O! = &, 
e — 

ct = A représentant les potentiels scalaire et vecteur du 
champ électromagnétique, 

Les constantes universelles E, et |, sont, respecti- 
vement, la constante diélectrique et la permeabilité du 
vide qui, dans le modéle quasi-stationnaire, seront alter- 
nativement annulées. 

Le paramêtre q est la conductibilité électrique, qui 
sera considerée comme une proprieté locale de la ma- 
tigre, transportée par celle-ci Le déplacement mécanique, 


—+ 
obtenu à partir de la vitesse V, nous permettra donc de 
connaitre à tout instant la distribution spatiale de 9 


—-— 
Finalement, le champ appliqué &º = Ea represente 
le champ électrique provennant d'interactions avec 


dautres systêmes de nature physique différente, comme 
les champs d'origine chimique ou comme le champ 
piezoélectrique, dorigine mécanique (pour des milieux 
plus compliqués que celui consideré). 


Pour l'electromagnétisme le champ Es sera donné 
en chaque point fixe de la matiére, ce qui nous permetira, 
la vitesse étant aussi donnée, de le considérer comme 
une fonction connue de l'espace et du temps: 


E = g(r t) (11) 


La condition c -—> * ou E, —> O permet la décom- 
position du systême d'équations en deux systémes plus 
simples, obtenus lorsque I'on fait tendre vers zéro cha- 
cune des constantes E,, |l,. 

Ainsi, si on annule la constante E, on obtient une 
limite magnétique des équations, qui se prête bien à 
la description des systêmes inductifs; quand on annule 
la constante |, on obtient une limite électrique, adaptée 
à la description des systêmes capacitifs. 

Pour voir comment s'écrivent explicitement les équa- 
tions correspondantes aux deux limites, considérons la 
formulation du type champ, donnée par les équations 
(3), (4). 


Compte tenu des équations (4-a, b), les deux limites 
ij 


E + 
signifient l'annulation d'une des induetions D ou B, de 
sorte que nous obtenons les équations, 


— dB => 
VXxE=— — VxE=0 (12-a) 
dt 
VxH=) ou VV. D=0 (12-b) 
— = e DD 
v. Bl VxH= Ji Go) 
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auxquelles s'ajoutent les relations: 


—s —% e mp 
B = Ho H D=8E 


ou 


(13-a) 


j=c(E+VxB+E?) ji=c(E+E9+p.V (13-b) 
La vitesse est supposée connue s'il s'agit de I'étude 
d'un systême purement électrique, ou donnée par les 
éêquations mécaniques s'il s'agit de lun systême électro- 
mécanique. 
Analysons plus en détail les deux limites obtenues. 


a) Limite magnétique 


Nous commencerons par faire remarquer que dans 
la limite magnétique la densité de charge p est nula, E] 


la suite de I'annulation de Iinduction électrique D. 


ll en résulte que le champ électrique É devient un 
simple auxiliaire de calcul, non determinê en dehors des 
conducteurs, mais fixé à lintérigur de ceux-ci à travers 
la loi d'Ohm. 

On pourra done l'éliminer de la formulation en intro- 
duisant la loi d'Ohm dans la premiêre équation de 
Maxwell et en ne gardant dans la formulation que les 
grandeurs vraiment caractéristiques de la limite magné- 

= o. —s 
tique, qui sont H| B et j, comme en magnétostatique 

Le potentiel vecteur A est determinê, dans cette 
formulation, mais le potentiel scalaire & subit le même 
sort que le champ électrique: il n'est determiné qu'á 
Vintérieur des conducteurs par la loi dOhm et cela à 
une constante prês. 

RFemarquons que cette indetermination du potentiel 
scalaire n'est que la traduction mathéêmatique d'un fait 
bien connu dans ['étude des circuits inductifs: l'impossi- 
bilité pratique de connaitre leur potentiel, qui doit être 
fixê en un point de chaque conducteur par des conne- 
xions à la terre. 


bj) Limite electrique 


= 
Pour la limite électrique, le champ magnétique H 


devient un auxiliaire de calcul, qui peut être éliminé 
en prenant la divergence de la troisiême équation de 


Maxwell; de même, la densité de courant j peut être 
eliminée en introduisant la loi d'Ohm dans la même 
equation. 

Les grandeurs caractéristiques de la limite électrique 
sont ainsi E D et p, comme en électrostatique, 

Pour des conditions frontiêres nulles, le potentiel 
vecteur À s'annule avec |'induction B. Le potentiel carac- 
téristique de cette limite est donc le potentiel scalaire &, 
complétement determiné par les êquations et conditions 
frontiéres. 

ll est intéressant de considérer le cas des condue- 
teurs parfaits, cestâ-dire, €-—s O, ce qui annule les 
temps de relaxation, 

Avec la condition 2 —> & la loi d'Ohm (13-b) donne 
immédiatement 


E+E1=0 (14) 
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ce qui montre un équilibre parfait entre le champ élec- 
trique et le champ appliqué, supposé connu, 
-—» 


La densité de courant j est alors indeterminée, ce 
qui n'affecte pas la connaissance de la distribution de 


aii 

charge électrique p (et donc de V.j) complêtement 
= 

determinée par la connaissance de E en fonction du seul 


e 
champ appliqué E*. 


2.2 — Analyse des équations quasi-stationnaires 
de la mécanique 


Considérons maintenant l'application des équations 
(1) à (8) à un systêéme mécanique, constitué par un 
fluide incompréssible dont les caractéristiques ont été 
precisées plus haut, 

Avant de passer à I'dentification des variables, nous 
allons faire sur la mécanique galiléenne, quelques remar- 
ques de caractére général, dont Iintérêt apparaitra 
plus loin. 

Ainsi, nous rappellerons qu'au niveau de la méca- 
nique galiléenne il y a une distinction, à priori, entre 
masse inertielle n et masse gravitationnelle ng, 'identi- 
fication des deux types de masse n'ayant lieu qu'áã un 
niveau plus avancé de la physique. 

Nous aurons besoin de cette distinction dans la 
formulation à developper par la suite. 

La deuxiême remarque concerne la notion classique 
de systêéme mécanique. Nous allons considérer comme 
intrinseque à cette notion la possibilité d'identifier une 
certaine quantité de matiére, ayant une certaine masse 
propre et de la suivre dans son mouvement. 

Il en résulte que le flux de masse est dú uniquement 
à un processus de convection (contrairement à ce qui 
se passe en électromagnétisme, ou I'existance de charges 
positives et negatives permet aussi un flux de charge par 
conduction), ce qui implique, finalement, la loi classique 


v.nV+ dn/dt=0 pour la conservation de la masse 
inertielle ou gravitationnelle. 

Nous verrons plus loin comment cette contrainte, qui 
n'existait pas en électromagnétisme ,va être incluse dans 
le formalisme de Maxwell et comment elle s'articule avec 
l'équation de conservation intrinseque à ce formalisme. 

Nous allons maintenant identifier les variables méca- 
niques intervenant dans les équations (1) à (8), ce qui 
doit être fait avec soin, puisque la formulation proposée 
est assez différente de la formulation classique: I'identi- 
fication n'est donc pas immédiats. 

Considérons d'abord I'équation (8-b), selon laquelle 


pa 
Q=Vv. 
ms 
Dans cette équation V est la vitesse, au sens ciné- 
matique, associée au choix du reférentiel galiléen et qui 


apparait aussi dans la loi d'Ohm électromagnétique. 
L'équation (8-b) transpose ce concept sur le champ 


Em. S 
physique et définit cy comme la vitesse, au sens méca- 


— + 
nique (on verra que € pourrait être consideré comme la 


Sa v2 
limite, pour c—>» %, de IV 1- , ce qui I'associe 
Cc? 


à la quadrivitesse). 
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Á partir d'ici, I'équation (6-a) nous Roma d'identifier 


le vecteur “2 comme étant la vorticité o du fluide, donc 
—», -—» 


RB = 0. 

La variable y peut être identifiée comme la densité 
de masse gravitationnelle, pour laquelle nous avons dejá 
adopté le symbole nx. On a doncg = 


EE 
Finalement, le champ .h a dejá été utilisé sous le 
nom anglais de «protational field» [6]. 


g —) —) mp : 
Quant aux autres variables 5, D, ;, 2 elles nont 


pas été utilisées, sauf pour la limite statique, à voir plus 
loin, ou * devient le potentiel gravifique de Newton, 


que nous representerons par Wy et ou E devient donc 
le champ gravitationnel. 

Avant le passage à la limite qui nous conduira à 
annuler I'une d'entre elles, les constantes universelles 
£o, Hoy sont, respectivement, I'inverse G"! de la constante 
(negative) de gravitation et la constante Gc*. Nous 
continuerons à utiliser, par raison de simplicité, les nota- 
tions £;, Ho. 

Le paramêtre & constitue une mesure de la viscosité 
du fluide et il sera défini par € = (,V)"! ou v est la 
viscosité cinématique. 

Remarquons que la notion de masse inertielle inter- 
vient dans ce paramêtre, puisque I'on a v = Y/m, ou Y est 
le coefficient de viscosité. 

Pour un fluide incompressible, la viscosité cinéma- 
tique reste constante sur chaque point matériel du fluide. 
L'évolution du paramêtre q, dans l'espace et le temps, 
est donc determinée par le déplacement mécanique, 
connu à partir de la vitesse, comme c'est le cas en 
électromagnétisme pour la conductibilité electrique, 

Compte tenu de la nature du paramêtre q, faisant 
intervenir la densité de masse inertielle m,, son évolution 
est, en fait, une consequence de la conservation de la 
masse inertielle, qui intervient ainsi dans la processus. 


Considérons, finalement, le champ appliqué ga, qui 
doit traduire interaction du systéme mécanique avec 
d'autres systêémes de nature physique différente. 

On y trouvera des termes comme, par exemple, le 
champ de la force de Lorentz intervenant dans une 
interaction électromécanique, 

ESA 


Mais le champ &º joue un rôle plus fondamental, en 
rapport avec la contrainte supplementaire dejã signalée 
concernant le flux de masse gravitationnelle. 

En effet, on doit inclure un terme capable d'assurer 
la nature purement convective du flux de masse gravita- 
tionnelle, c'est-â-dire, d'assurer que dans lI'équation de 


— 
conservation de Maxwell (6-d) on ait V.;j= V.(ngV), 
soit, compte tenu de la loi d'Ohm (7-c): 

+ > 3 Ss 


(58 +Vx8B+E)=0 (15) 
Nous postulerons que ce terme existe et est le 
champ de pression—V w =— V p/n,, ou w est I'enthal- 


pie et p la pression, Ce terme n'est pas donné à priori, 
comme les autres termes possibles d'interaction, mais 
doit être obtenu simultannement avec I'évolution méca- 
nique du systéme, Ceci introduit justement un degré de 
liberté supplémentaire dont on a besoin pour obtenir 
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avec les équations de Maxwell un flux de masse gravi- 
tationnelle purement convectif. 

Aprês avoir procédé à l'identification des variables, 
nous allons étudier les deux limites imposés par la 
condition E, Hy > 0, 

Pour écrire explicitement les équations qui corres- 
pondent à ces deux limites, nous allons utiliser la formu- 
lation en fonction des potentiels données par les équa- 
tions (5) à (7), car elle conduit plus directement aux 
résultats souhaités. 

Nous aurons ainsi une limite cinétique, obtenue en 
annulant la constante £,, donc la densité de masse gravi- 
tationnelle ng; et une limite statique, obtenue en annu- 


pn 
lant la constante |, donc la vorticité O. 
Pour abréger, nous utiliserons tout de suite I'équa- 


Rad, 
tion (8-b) en remplaçant le potentiel vecteur c] par la 


——» 
vitesse V Cela nous donne les équations: 


-—+> o —s — 
V XV =n V XV = (16-a) 
—+ 
NV “+ — 
VOO +— =— & ou VDP=— 3 (16-b) 
t 
— —», 
VV =0 VV = (16-c) 
—> —, ó Tg 
auxquelles s'ajoutent les équations 
em er 2 1] 
VV =— o j ViP=— Do (7a) 
, E, 
ou 


j = :(8 + VXO+8%) = c(E+E)+nGV (17-b) 
Analysons un peu plus en détail ces deux systêmes 
d'équations. 


a) Limite cinétique 


Cette limite, obtenue en annulant la constante €,, 
correspond à la limite magnétique de I'électromagnétisme, 
ce qui est en accord avec l'analogie classique de Maxwell 
pour les circuits. 

À cette limite la masse gravitationnelle est annulée, 
mais la masse inertielle subsiste dans q et dans w. 

Le potentiel scalaire 9) ias plus determiné par 


'équation (7-b). En effet, 2? et é sont devenus des auxi- 
liares de calcul, determinés par la loi d'Ohm (17-b), mais 
que I'on peut éliminer de la formulation, si on veut seule- 


ment connaitre les grandeurs caractéristiques du systêéme 
—s —» ra 


qui sont V =4, je o. 
E. 7 
En réalité, classiquement, le vecteur j nest pas 
utilisé en mécanique des fluides. Il peut être mis en 


rapport avec la conservation du moment angulaire, mais 
cette question est sans importance pour les fluides iso- 
tropes que nous considérons ici, 

Nous allons donc éliminer de la formulation non 


—», —», 
seulement 2? et &, mais aussi le vecteur ;. 
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Pour cela, on introduit (16-b) et (17-a) dans la loi 
d'Ohm (18-b), ce qui donne, en prenant le rotationnel: 


— >» 


—»  —s>» —s 
+Vx(Vx0)+Vxsº=0 (18) 


qui est I'équation de diffusion de la vorticité, 
-—» 
Si on réduit & * au terme de pression —V w on a 


V XE =0, ce qui conduit à lI'équation de diffusion 
classique. 
Le comportement du systême est décrit par I'équa- 
tion (18), plus les équations (16-a, c), non encore utilisées. 
Remarquons, finalement qu'il y a dans la formulation 
une équation de trop, comme en électromagnétisme: 
l'équation (16-d) est superflue et peut s'obtenir à partir 
de (16-c) et de (17-a). 


Une fois o et V determinés par les équations prece- 
dentes, on peut essayer d'exprimer O à partir de la 
loi d'Ohm. Pour cela, on élimine de cette derniére é 
et j, par les équations (16-b) et (17-a), et on  SbHont 
l'équation de diffusion du potentiel vecteur = 

MY 


+ 
: D=—>Vw+vV2V (19) 
t 


On reconnait immédiatement Iéquation classique de 
Navier-Stokes. 


La connaissance de v nous fournit donc V (9 + w) 
soit, à une constante prês à determiner par les conditions 
frontiêres, la somme des fonctions 2º + w, 

La séparation entre ces deux fonctions n'est pas 
necéssaire, au niveau des fluides incompréssibles, puis- 
qu'elles ne figurent dans aucune équation autre que 
l'êquation (19). 

Mais si nous voulons garder pour I'enthalpie w la 
signification habituelle, nous devouns poser, en accord 
avec les relations classiques: 


Vº + Quo (20) 


ou Yno sera le potentiel de Newton d'origine extérieur, 
c'est-à-dire, determiné par les conditions frontiêres. Mais 
ce dernier est nul dans notre cas, puisque nous avons 


—s» 
adopté comme conditions frontiêres 9', V-50, 


É À cus 
Remarquons que l'expression 2º = — V? qui provient 


de l'approche classique de l'étude des fluides et qui 
doit être postulée dans notre description des fluides 
incompressibles, possêde, cependant, une signification 
assez générale. En effet, on peut montrer que dans une 
couche de transition entre deux vitesses, par exemple 


—>» —. ——+ - 
Vo = 0 et V, le vecetur j crée un effet de double-couche 
responsable d'une variation de pontetiel scalaire qui est, 


1 
justement, Y— 9) = —v2 
9 2 
En résumé, nous pouvons dire que la limite cinétique 


correspond bien à |'hydrodynamique classique, Les équa- 
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tions sont, simplement, écrites autrement, grâce à l'intro- 
e 
duction de deux grandeurs auxiliaire 5, 2) qui peuvent 


-—» 
étre éliminées et d'un vecteur intermédiaire f tel que 
cc) “ 


; — 
) =—vyV2V represente le champ des forces de vis- 


G 


cosité. 
b) Limite statique 


Cette limite, qui est l'analogue de l'électrostatique, 
sobtient en annulant la constante |, 


La masse gravitationnelle subsiste dans le modéle 
et on n'aura plus besoin de la distinguer de la masse 
inertielle: Ng = Ni. 


—s» 

Compte tenu des conditions frontiêres, annulant V à 

Vinfini, les équations (16-a) et (16-b) nous donnent comme 
% 


solution V = O, 


Nous avons donc un systême au repos, 


Les équations fondamentales sont (16-b) et (17-a), 


o 


qui nous montrent que le potentiel scalaire '“ devient 


a) 
le potentiel de Newton Wy et que le champ & devient 
le champ gravitationnel, 


Pour analyser la loi d'Ohm (17-b), nous devons com- 
mencer par constater que, d'aprês l'expression dO = 
= (11, v)"!, le paramêtre o tend vers linfini, puisque la 
limite statique correspond à I'annulation de la cons- 
tante H,. 

Il en résulte Iéquilibre entre le champ gravitationnel 
et le champ des forces d'interaction: 


—, = 
E+ 8º = (21) 


my 
À cause de la nature du terme figurant en &º, V w, 
qui nest pas connu à priori et que I'on ne peut même 
pas expliciter en fonction des variables mécaniques du 
fluide incompréssible, la répartition de mn; ne s'obtient 
pas à partir des équations precédentes, comme c'était 
le cas en electromagnétisme. Elle est obtenue directe- 
PR. * 
ment à partir de la condition (15) qui donne V. j = 
me) . 
= V.(no V) = 0, donc: 
d 
Ee = 0 (22) 
t 
comme il se doit pour un systême parfaitement statique. 
E 
Ce résultat montre que j est, comme pour la limite 


analogue en électromagnétisme, une variable intermédiaire 
qui peut être eliminée de la formulation, les variables 


e 
fondamentales étant 9? & et ng. 
Signalons, comme remarque finale, que nous n'avons 


ao q 
pas voulu déterminer ;, mais celã serais possible si |'on 
avait posé, au lieu de (15), la condition plus restrictive: 


— -—», er dA — + -—3 
hO(E+VXB+E)=—V?2V (15-a) 


E = eli 
ce qui nous donnerait alors j=V = 0, 
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Mais cette condition (15-a) serait difficile à justifier, 
au niveau quasi-stationnaire, ne devenant claire qu'à un 
niveau beaucoup plus avancé, 


2.3 — L'interaction électromécanique 


Dans les paragraphes precédents nous avons analysé 
les équations électriques et mécaniques, sans nous préo- 
cuper du fait que les systêémes électrique et mécanique 
aient été isolés ou en interaction, 

Precisons maintenant ce qu'il y a lieu de faire si les 
deux systêmes sont en interaction, 

La premiére chose sera l'identification des vitesses 
figurant dans les équations électriques et mécaniques, 
comme dejáã signalé, de sorte que I'on posera 


—), —», 
Vo=VYm=V (23) 


ou les indices e et m signifient qu'il s'agit des vitesses 
intervenant, respectivement, dans les équations électri- 
ques et mécaniques. 

Le deuxiême aspect de |'interaction concerne |'intro- 


duction, dans les champs appliqués sê des termes 
d'interaction. 

Pour les milieux simples considerés, il ny a pas de 
champ appliqué électrique en provenance de |'interaction 


électromécanique. Donc: 


= 0 (24-a) 


Par contre, du coté mécanique il y a un champ de 
force d'origine électromagnétique, que nous postulerons 
comme étant dú à la force de Lorentz, 


e) =) = ep 
La densité de force de Lorentz est f=;jxB+pçE 
ou les deux termes ne coexistent pas lors du passage à 


EPA. * 
la limite, Il en résulte un champ appliqué &2. tel que: 
e e A à 
Mm 8:=)xB+pE-V (24-b) 


ou ny répresente la densité de masse inertielle, 

L'évolution du systêéme électromécanique est complê- 
tement determinée par les équations électriques et méca- 
niques, completées par les relations (23) et (24). 

Nous n'allons pas nous étendre sur |'étude de 
interaction électromécanique, dans la mesure ou les 
conclusions sont classiques et independantes de la formu- 
lation particuliêre que nous avons adoptée. 


3 — FORMULES DE TRANSFORMATION 


Les équations precédentes, comme toutes les équa- 
tions de la physique, doivent être invariantes lors d'un 
changement de référentiel, par translation uniforme. 

Considérons, par exemple, les équations (3-a) et 
(3-b) supposées écrites pour un certain reférentiel R: 


Ra. 
-—) O 38 
Vxg nesponaiA (25-a) 
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V. 8=0 (25-b) 


Dans un autre référentiel R', en translation uniforme 
par rapport à R, on doit avoir: 


— + 
v' e « o sa" 26. 
X 8 =— a 
N ER (26-a) 
— 
VW. OD =0 (26-b) 


ES eis 4 
ou s' s3' sont les champs et V ET les operateurs 
différentiels dans le référentiel R', 


L'invariance des équations impliquera normalement 
= my 


que les grandeurs E, &B se transforment elles-mêmes, 
lors du changement de référentiel, c'est-a-dire, que 
rd =— = 
Ss FSe a EM. 

Pour trouver les lois de transformation galiléennes, 
qui nous interessent, nous partons des lois de transfor- 
mation de l'espace-temps de Galilée, qui sont 


” 


ff = 


os. — 
= 


= 
Vat (27-a) 
f=t (27-b) 


' +» = = P 
ou r, r est le vecteur position; t, t le temps; et V, 


la vitesse de translation du référentiel Rº par rapport à R. 
À partir des relations (27) nous pouvons établir les 
lois de transformation des opérateurs différentiels V et 


———. lors du changement de reférentiel, On eétablit, sans 
t 
difficulté [3]: 
Vt= V (28-a) 


+Vo Y (28-b) 


Si on introduit ces résultats dans les équations (25) 
et que Ion compare à (26), on obtient 


RE FR NE SE = a: 
V'x[5'—s + Vo X88)] =— a [8'—s] (29-a) 
t 
ey —% 
Ve. [8—s]=0 (29-b) 


Les relations (29) étant valables dans tout [espace 
et à tout instant, il est habituel de postuler tout de suite 
que les lois de transformation sont celles qui annulent 
les termes entre crochets, c'est-á-dire, B'- 2, etc, 

Mais nous devons analyser plus soigneusement les 
lois de transformation, à cause des difficultés qui vont 
apparaitre au niveau de la mécanique. 

Nous poursuivrons donc [étude precédente, pour les 
équations (3-b) et (3-c), ce qui donnera: 


E E 
Vi [O -0] = [94] (30-a) 
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VOX [5 —(b—Vo X =; Ls hig Vo) ]+ 


É = =s 
+— [9'— 9] [80-b) 
dt 


La loi d'Ohm nous permet d'obtenir la relation 
-—» —+ —> er a A 
Ly 4 9V)] = (LE —(S+Vo x 38) + 


my E. Es | 


= EA F 5 E 
ou V'=WV-V, et ou on a utilisé I'hypothêse habituelle 


Finalement, les relations (4-a) et (4-b) nous donnent 


[D'—o]=s[8—6] (32-a) 
[5838] = 1, [56] (32-b) 


Compte tenu de l'approximation c — “2, nous savons 


— 


4 —& ih à 
dejá que les champs 9 Det 8 8 ne coexistent pas. 
On peut donc remplacer (32) par: 


['—D] = [8—(6 + Vo, X88)] (33-a) 
[88 —88 | = Bo [25 — [Vo XD)] (33-b) 


Si nous observons les équations (29) à (33), nous 
constatons qu'elles forment un ensemble d'équations de 


Eai 

Maxwell, pour les champs entre crochets LE-(E+ Was: 
= = —% 

xe8)]. [D'—D], ete. 

Si nous faisons maintenant |hypothêse heuristique 
que les conditions initiales ainsi que les conditions fron- 
tiêres pour ce nouveau champ sont nulles, on conclut 
que la solution sera nulle, puisqu'il ny a pas de sources 
dans le systême (supposé stable), 

|| en résulte, donc, les relations 


RR E | 
E=E+V, x 48 (34-a) 
mm -— 

N'= q (34-b) 

= = cd e 

do =—V, XD (34-c) 
e —+ 

SB =48 (34-d) 
P=5+ 9, (34-e) 
gq = q (34-f) 


qui sont les lois de transformation lors du changement 
de reéférentiel. 

Les lois precédentes sont en accord avec |'expé- 
rience, aussi bien en électromagnétisme qu'en mécanique 
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quasi-stationnaires et peuvent donc être acceptées comme 
les formules de transformation valables dans le cadre 
galiléen. 


Les difficultés annoncées au début vont apparaitre 
avec les formules de transformation des potentiels 9º 

Sa. + 
et cl, que nous allons établir maintenant. 

Pour cela, appliquons la méthode d'êtude precé- 
dente aux équations (6-a) et (6-b). Aprês des calculs 
relativement simples on arrive aux équations: 


E PR = po 
V'x[a—a] = [38'—s8] (35-a) 
A DE q É A 
VLW -v,.9] = LE—8] s- 
—+ —p my —— 
— [8'—(& + Vo X68)] (35-b) 


qui sont les êquations des potentiels associês aux champs 
des équations (30) à (33). 

Si on fait Ihypothéêse que les conditions initiales et 
les conditions fontigres sont nulles pour ces potentiels, 
on arrivera aux formules de transformation cherchées en 
annulant les termes entre crochets des equations (35). 


Les résultats obtenus sont corrects pour lI'eléctro- 
magnétisme, ou nous aurons, en effet: 


—> e 

e =Qe (36-a) 
pe ERR 

9 e = We—NVg “Ele (36-b) 


Par contre, en mécanique, les équations (36) sont tout 

à fait fausses. En effet, il suffit de rappeler I'équation 
—> Eee 

mécanique Ym = V et expression générale du potentiel 

scalaire 9º =—V* + qy, pour conclure que l'on doit 


avoir les formules de transformation 


> a 

El rm =E m—Va (37-a) 
em ul 

Ps ja q mn Voa ri: (37-b) 


— —s 
Cela signifie que les potentiels []— a] et 
Do amido 


[9"—(2—V, .G)] ne sont pas nuls, mais constants et 
qu'ils ont les valeurs respectives —V, et = 


L'origine de ces difficultes, analogue à celles rencon- 
trées quand on étudie les équations de conservation de 
lenergie et de ['impulsion [3], est ici immédiatement 
detectée. 

Elles viennent de l'hypothêse que les conditions 
frontiéres et initialles sont nulles pour les potentiels 

——h 
[a] et [ar-(9P —v,. A]. 
> — 

Or, compte tenu des équations Ym=V et Wm= 

add + Px un changement de référentiel doit être 
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acompagné d'un changement de conditions frontiêres sur 


les potentiels 9" et e! 
Cela veut dire que les conditions frontiêres 


e 

Elm > O DO (38) 
ne sont valables que pour le référentiel ou la matiére 
se trouve au repos à |'infini, 


E + 

Pour tout autre référentiel ,animé d'une vitesse V, 

par rapport au premier, les conditions frontiêres doivent 
être 


= 1 
Elm —Va Vos —Vs (39) 
e 
ou —V, représente, finalement, la vitesse de la matiére 
à Vinfini. 
Si nous postulons les conditions frontiéres (39) pour 
les potentiels, au lieu des conditions particuliêres (38) 
initialement adoptées, les formules de transformation (37) 
sont alors justifiées vis à vis des formules (36) obtenues 
pour l'électromagnétisme, Les termes supplémentaires 
constants qu'elles contiennent proviennent du changement 
des conditions frontiêres, que ajoutent justement aux 


o e ” E mp mma 
potentiels euApInicas La'—g] et [9 —Vo-O)] les 
—+ ; 


valeurs —V, et a ve 


Les champs ne sont pas affectés par des constantes 
figurant dans les potentiels, ce qui explique pourquoi les 
formules de transformation correspondantes étaient cor- 
rectes. 

La solution qui consiste à postuler les conditions 
frontiêres (39) nous parait être la seule possible, dans 
un cadre strictement galiléen, mais nous pourrions envi- 
sager d'autres maniéres de procéder, que nous verrons 
par la suite. 


4 — DIFFICULTES D'UNE EXTENSION RELATIVISTE 
RESTREINTE 


L'utilisation des conditions frontiêres (39) n'est pas 
tout à fait satisfaisante, en particulier parce que les 
conditions apparaissent comme arbitraires, à priori, et ne 
trouvent une justification qu'aprês avoir determiné pour 
les potentiels des expressions en fonction de la vitesse. 

| est possible de retrouver les mêémes résultats par 
une autre voile, qui a l'avantage de rendre plus claire la 
question des conditions frontiéres, mais qui a |'inconve- 
nient de sortir du cadre strict de la physique galiléenne 
et devient donc plus compliquée. 

Le modéle que nous avons presenté pour |'électro- 
mécanique, constitue le passage à la limite c-—> YU du 
modéle de Maxwell pour |'êlectromagnétisme, qui est au 
dêpart invariant pour la transformation de Lorentz. 

Ainsi, au ligu détudier les formules de transforma- 
tion aprês le passage à la limite, comme nous |'avons 
fait, il serait possible d'établir les formules de transfor- 
mation d'abord, en travaillant avec la transformation de 
Lorentz au liceu de celle de Galilée et de faire ensuite 
le passage à la limite, 


6/ 


En adoptant cette vole, on retrouverait sans difficulte 
les formules de transformation des champs (34) qui sont 
simplement les limites pour c -—» “O des formules relati- 
vistes correspondantes [7]. 

Voyons en quoi cette méthode pourrait éclaircir la 
question des conditions frontigres pour les pontentiels, 

Pour ce qui est des potentiels de nature électro- 
magnétique, on continuerait à poser les conditions fron- 
tigres 


—+ ” 


et les formules relativistes nous donneraient, aprês le 
passage à la limite 


=p — 
D=D-V. A (41-a) 
—» — 
E = Ele (41-b) 


qui sont les résultats delã connus. 

Pour ce qui est des potentiels de nature mécanique, 
on utiliserait comme conditions frontiêres dans le réfén- 
tiel ou la matiére se trouve au repos à infini, non plus 

Avec ces conditions frontieres, la solution galiléenne 


na . ar 
um 


0 (42 


Avec ces conditions frontiêres, la solution galiléenne 
pour le potentiel scalaire deviendrait On =YPN+ 


: 
+—Vº + e, Lors d'un changement de reéférentiel, les 
Ê 


formules relativistes donneraient en passant à la limite 


/ 
e V E 
et en utilisant lim c [V 1— — =ce+—M: 
c—s E Cc” 2 
=> 1 
E Tm “= m Va > 
Es q E r 4 
Tmn= lim 2 | E 
A 6 =3500 = Va Le E Va 
cê cê 
—, Ea. ” 
= Ca tea 
(43-a) 
Ps RR 
Dm E Em a Wa 
Va — lim j ve vê mg, 
ASA fra —— 
CS cr F, 
, —,  —% É cá 
ate Vs +—VS 
2 
(43-b) 


qui sont bien les résultats corrects. 

Lors d'un changement de référentiel les conditions 
frontiêres (42) se transformeraient selon les mêmes for- 
mules (43) et on obtiendrait: 


' 1 1 z a = 
Va E + Mo En —V, (44) 
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On retrouverait ainsi directement les changements 
de conditions frontigres que nous avons dú postuler au 
niveau galiléen, 

On voit, donc, quil est possible de rendre plus 
claire, la question des conditions frontiéres et des for- 
mules de transformation pour les potentiels, 

Cependant, malgré cet avantage, la méthode nous 
parait avoir des inconvénients majeurs, au niveau galilêen, 

Tout d'abord, à cause de son excessive lourdeur et 
de la necessité de sortir du cadre simple de la transfor- 
mation de Galilée, en passant au cadre plus complexe 
de la transformation de Lorentz, 

En second lieu, parce que on doit admettre un poten- 


| E i md! 
tiel scalaire “2 


1 = = 
= py + Ei + c qui, à la limite, tend 


vers |'infini avec le terme c* (même si une telle proprieté 
disparait dans les équations du systéme, eliminée par 
les opperateurs differentiels qui agissent sur (4). 

Par contre, la méthode serait intéressante si on 
pouvait formuler une électromécanique du niveau de la 
relativité restreinte ou on n'iimposerait plus la condition 
c-— % et, en particulier, ou on aurait plus besoin de 
distinguer la masse inertielle de la masse gravitationnelle, 
dans la mesure oú cette derniêre ne serait jamais annulée. 

Dans ce cadre, l'électromécanique quasi-stationnaire 
apparaitrait comme une simple approximation de cette 
formulation plus générale, point de vue adopté vis-á-vis 
de lélectromagnétisme par la majorité des auteurs, qui 
néglige tout souci de formuler un électromagnétisme 
galiléen. 

Malheureusement, il nest pas possible d'établir une 
formulation satisfaisante de | électromécanique au niveau 
relativiste restreint. 

En effet, le modéle proposé dans le paragraphe 2, 
qui peut être geênéralisé en électromagnétisme par sup- 
pression de la condition c-—» * (ce qui est d'ailleurs la 
forme habituelle de I'électromagnétisme), ne peut pas 
subsister en mécanique sans cette condition. 

La difficulté vient de la seule équation qui distingue 
les formulations galiléennes de ['électromagnétisme et de 
la mécanique, c'est-á-dire de la condition supplementaire 


de la mécanique a. = V laquelle ne peut pas être géné- 
ralisêe au niveau de la relativitê restreint. 

En effet, si on traite les équations de la mécanique 
d'un point de vue relativiste, par analogie avec |'électro- 
magnétisme, on arrive à la conclusion que le potentiel 
vecteur devrait être de la forme 


—+ 
—s —»+ V 
Elm = di na a (45) 
Vo 
Tom ; a) 
Er 


ou le terme 2, est un potentiel retardé analogue au 
potentiel vecteur de ['électromagnétisme et est associé 
à un potentiel gravitique de Newton qu. 
Or le terme “y nest pas lié à la vitesse par une 
= 
relation locale, ce qui ne permet pas de définir V en 


= 
fonction de Ela: 
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La relation (45) ne peut donc pas constituer I'êqua- 
tion supplêmentaire dont.-on-a besoin en mécanique. 

Pour obtenir cette équation à partir de (45) il nous 
faut, en fait, annuler le potentiel an. 

Une façon de le faire est justement introduction de 
la condition c-—> 2, qui ne permet pas la coexistance 
du courant de masse garvitationnelle (qui est ici la masse 


potentiante) et du pontetiel EN, dans la mesure oú cette 
condition annule, soit la masse, soit la constante L.. 

Il en résulte que l'équation (45) peut être réduite 
au second terme lequel donne à la limite 


== 


qui est justement ['équation utilisée. 

Une deuxiême façon de procéder, mais qui ne béné- 
ficie pas de la cohérence de la precédente, est de 
négliger de façon arbitraire la partie des potentiels repre- 

—% 


sentée par Qy et 9y: 


Py = O Uy = O (47) 


ce qui transforme (45) en une relation locale permettant 
Hétude de la mécanique. 

Dans ce cas on peut retrover à partir de la théorie 
de Maxwell les équations de la mécanique relativiste 
habituelle (restreinte) qui compte tenu de l'approximation 
(47) ne nous parait pas satisfaisante pour une généra- 
lisation cohérente de l'électromécanique galiléenne, 

Une fois exclue la possibilité d'une électromécanique 
satisfaisante au niveau relativiste restreint, il nous reste 
une troisiême voie qui serait d'avancer dans le sens de 


= 
la relativitê général, en absorbant les potentiels as, Px 


dans la métrique. 
Pour montrer que cette démarche serait possible, 
fL 
introduisons un quadricourant ] tel que 


[a M 


a mx: 
J=n,u —E,A 8 (48) 


ou mn, est la densité de masse propre; uo la quadrivitesse; 
et le demier terme, une généralisation du «courant de 
conduction», défini par (15-a). 

Le problême est de trouver une métrique qui rende 
compatibles: — |'équation des potentiels, qui s'écrit [8] 


A” sm A ed 49 


[Ed ti 
ou A est le quadripotentiel et R, le tenseur de Ricci; 
— et la condition 


AE it (50) 


qui est la généralisation de l'équation (45), mais avec 
= 


dq = O et px = 0. 
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Si l'on introduit (48) en (49) et I'on tient compte de 
expression du tenseur impulsion-énergie: 


E a 


rd 
Ta=n,€ u ug (51-a) 


(51-b) 


on obtient immédiatement, par utilisation de (50) e de 


Uu ==] 
a 8 E vê E) 
RÁ e dy udft) A 
i -— E, €º p e É 
(52) 
qui est vérifié si 
a É eo o 1.A 
Ro= Pepe 
É — E, C1 ( Er str) 
(53) 


Cette équation est I'équation d'Einstein, déterminant 
la métrique de la relativité générale, qui devrait donc 
permettre d'étendre notre modéla d'une façon cohérente, 
ce qui ne serait pas possible dans le cadre relativiste 
restreint. 

En conclusion, une fois exclue la possibilité d'une 
electromécanique satisfaisante au niveau relativiste res- 
treint, on se trouve devant Ialternative, soit de se con- 
tenter dune formulation située dês le départ au niveau 
galiléen soit d'aller plus loin que la relativité restreinte, 
dans le sens de la relativité générale, 

Cette derniêre voie parait injustifiée, si le but pour- 
suivi est de revenir, de toute maniêre, au niveau quasi- 
-stationnaire, en vue des applications. 

C'est pourquoi nous préferons la formulation complêé- 
tement galiléenne, qui reste cohérente et commode, mal- 
grê les quelques inconvenients qui sont apparus plus haut. 

En tout cas, et cela nous parait une conclusion 
importante, le niveau relativiste restreint, normalment 
adopté en électromagnétisme, est inadequat pour |'étude 
de I'électromécanique et ne devrait pas être adoptê, 
notamment pour |'enseignement de cette branche. 


5 — CONCLUSION 


Le travail a examinê la structure des équations 
quasi-stationnaires de lI'électromécanique, 

Nous avons montré qu'il est préférable de définir 
le modeéle quasi-stationnaire directement au niveau gali- 
léen, plutôt que de le considérer comme une approxima- 
tion d'un modéle plus gênéral. 

Pour arriver à cette conclusion nous avons developpé 
un modéle de la mécanique qui permet une comparaison 
directe avec l[êlectromagnétisme et donc un traitement 
en paralléle des deux branches, 

Certaines difficultés concernant les lois de transfor- 
mation en mécanique galiléenne, et qui avaient deja été 
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mises en evidence pour lenergie et [impulsion [3], sont 
ici abordées et expliquées au niveau plus fondamental 
des champs et des potentiels. Lenergie et [impulsion 
étant, au niveau pré-relativiste, de simples fonctions d état 
associées à ces demiêres grandeurs, leurs formules de 
transformation peuvent être considerées comme justifiées 
à partir de celles que nous avons obtenues, 

Nous nous sommes limitêés à considerer des milieux 
três simples et dans des conditions peu complexes. 

En particulier, nous avons consideré, pour la meca- 
nique, un milieu indefini et au repos à [infini dans un 
referentiel, qui est, en ce sens, un referentiel privilegiê. 
On aurait pu se libérer de cette contingence, en consi- 
derant des milieux bornés, avec des conditions fron- 
tiêres quelconques. 

Mais ces complications supplêmentaires nauraient 
changé ni les résultats ni les conclusions auxquelles nous 
sommes arrivês, 
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Etude de la commutation 
dans une génératrice polyphasée à redresseurs 


Premiére partie: 


Analyse de la génératrice à couplages sinusoidaux 


RESUMO 


Este artigo é a primeira parte de um estudo da 
comutação dum gerador polifásico sem colector À partir 
do estudo de dois casos limite, o «sinusoidal» apresen- 
tado nesta primeira parte e o «triangular» apresentado 
na segunda parte, mostra-se que a forma das indutâncias 
mútuas, em função da posição, pode afectar considera- 
velmente o processo de comutação e, por conseguinte, 
a caracteristica externa do gerador. 

Resultados experimentais confirmam as previsões 
teóricas. 


INTRODUCTION 


Le développement des semiconducteurs de puissance 
a conduit à remplacer, dans certaines applications, le 
collecteur classique par un pont de redresseurs ou collec- 
teur électronique. Les génératrices à diodes tournantes 
des systêmes d excitation brushless constituent un exem- 
ple typique de cette technologie. 

Si le mécanisme de la commutation a été largement 
étudié dans le cas des génératrices à redresseurs tripha- 
sées [1,2], il nen va pas de méme pour les genératrices 
à nombre plus élevé de phases qu'on trouve actuelle- 
ment dans certains systêmes d'excitation [3]. 

C'est à l'analyse d'une telle génératrice qu'est con- 
sacrêe la présente étude. On considére le cas d'une 
machine à pôles lisses, sans amortisseurs, à induit poly- 
gonal et nombre pair de phases (fig. 1). Le nombre de 
phases doit étre supérieur ou égal à 6 pour que la 
presente étude s'applique sans aucune restriction, La 
connexion en polygone des phases de |'induit est imposée 
par des raisons d'utilisation optimale des enroulements. 
Afin de montrer [influence sur le comportement, de la 
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ABSTRACT 


This paper is the first part of a study on the commu- 
tation processes in a polyphase commutatorless D.C. 
generator, From the study of two limit cases, the «sinu- 
soidal» one presented here and the «triangular» one 
presented in the second part, it will be shown that the 
form of the variation of the mutual inductances as a func- 
tion of the relative positions of the different coils can 
affect considerably the commutation processes and hence 
the generator output characteristic. 

Experimental results on laboratory test machines 
confirm the theoritical analysis, 


forme des lois de couplage entre les enroulements, on 
considêre deux cas extrêmes. 


La premiére partie présentêe ici est consacrêe à 
VFetude du premier cas: 


— couplage sinusoidal «parfait» * 


Finduit 


entre les phases de 


— couplage sinusoidal entre les phases de Induit 
et Iinducteur 
Dans une seconde partie, on examinera le second 
cas: 


— couplage triangulaire entre les phases de |'induit 
et Inducteur 


* Par couplage «parfaits, on entend un couplage dont la loit est 
telle que la mutuelle entre deux phases occupant la méme position 
géométrique est égale ,en valeur absolue, à la self propre de 
chacune de ces phases. 
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— couplage triangulaire «parfaits entre les phases 


de l'induit. 


Le premier cas revient à considérer une machine 
avant des enroulements de type alternateur et le second 
une machine possédant un induit de type dynamo et une 
force magnétomotrice d'inducteur quasi rectangulaire, 


L'étude présentée est effectuée dans le cas du 
régime permanent, Une solution analytique du fonction- 
nement depuis la marche à vide jusqu au court-cireuit 
est obtenue en faisant les hypothêses et approximations 
suivantes: 


— |inducteur est alimentê par une source de courant 


— [8 courant de sortie du pont est parfaitement lissé 


— la résistance interne des phases de l'induit est 
négligeable 


— la chute directe et le courant inverse des diodes 
sont négligeables, 


1. EQUATIONS DE DÉPART POUR LA MACHINE 


A COUPLAGES SINUSOIDAUX 


Avec les notations de la figure 1 et un courant 
dinducteur constant, la force électromotrice engendrée 
par l'inducteur dans chacune des n phases de I'induit 


axe direch 
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peut seécrire si [on prend comme origine linstant de 
couplage maximum entre la premiére phase et [inducteur 
(passage de cette phase par l'axe neutre relatif à [in- 
ducteur): 


ex = ÔMçirsin [0-2 (K— 1) m/n) ISKSn 


Dans ces équations, M, représente la valeur maximum 
de la mutuelle entre linducteur et une phase de [induit, 
0 la position angulaire du rotor, O sa dérivée par rapport 
au temps et n le nombre de phases, Les mutuelles entre 
les phases de l'induit sont de la forme: 


M;=M,cos[2(i-)=)n] 1<i,j<n 


(2) 


Dans ces équations M, reprêésente la valeur maximum 
de la mutuelle entre deux phases. Conformément à la 
définition du couplage parfait, cette valeur coincide avec 
celle de la self propre, L, de chaque phase, 


2. ANALYSE DU FONCTIONNEMENT 


En passant de la marche à vide au court-circuit, on 
rencontre plusieurs modes de fonctionnement, chacun 
correspond à un mécanisme de commutation différent. 


axe neufre. 


FIG. 1 
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C'est à l'étude des différents modes qui se succêdent 
et aux conditions de passage dun mode au suivant 
qu'est consacré le présent paragraphe. 


21. Fonctionnement à faible charge (mode 1) 


Pour une telle charge, le comportement du redres- 
seur est le suivant (fig. 2a). A un moment donné, la 
borne de sortie positive du redresseur est relige, par la 
conduction d'une des diodes à cathode commune, au 
sommet du polygone d'induit dont le potentiel est le plus 
élevé. D'une maniére analogue, la borne de sortie néga- 
tive du redresseur est reliée, par la conduction d'une 
diode à anode commune ,au sommet du polygone d'induit 
dont le potentiel est le plus bas. La connexion ainsi 
réalisée partage |'induit en deux voies, mises en parallele, 
une (voie positive) comportant toutes les phases dont 
les forces électromotrices sont positives, lautre (voie 
négative) toutes les phases dont les forces électromo- 
trices sont négatives, Dans le cas du polygone à nombre 
pair de phases, considéré ici, chaque voie comporte 
n/2 phases et le courant se partage de maniére symé- 
trique entre les deux voies, Avec les sens de références 
choisis, on a+l,,/2 dans chaque phase de la voie posi- 


tive et —|.,/2 dans chaque phase de la voie négative. 


du courant qui les traverse et le déplacement, dun 


sommet au suivant, des liaisons du polygone d'induit 
avec le circuit extérieur. Compte tenu des inductances 
propres et mutuelles des phases considérées, l'inversion 
du courant ne peut s'y faire instantanément. Il y a un 
court-circuit transitoire de ces phases à travers les diodes 
qui doivent cesser de conduire et celles qui entrent en 
conduction (fig, 2b). Ce phénoméne constitue |'empiête- 
ment. Sa durée correspond à une rotation du rotor d'un 
angle, appelé angle d'empiétement, 

Comme on se limite à êtudier un fonctionnement en 
régime permanent, il est évident que [evolution qui vient 
détre décrite, nest possible que si langla dempiétement, 
soit Y, est inférieur à I'intervalle angulaire séparant deux 
commutations successives, soit 2 7/n. 

Durant l'empiétement, par suite du couplage existant, 
les tensions aux bornes des phases non commutantes 
sécartent des forces électromotrices correspondantes 
induites par lI'inducteur. Toutefois, comme on le verra, 
cette modification, due à la variation du courant dans 
les phases en commutation, nentraine lentrée en con- 
duction daucune diode supplémentaire, 

Afin de mettre le phénomêéne en équations, consi- 
dérons la commutation qui transfêre la phase 1 de la 
voie négative à la voie positive et la phase (n/2+1) de 


a) «Ee y<9<0 


b) 


0<0<y Cc) p<0<+r 


FIG. 2 


Comme le courant de charge est supposé parfaite- 
ment lissé, et que la résistance des phases de [induit 
est négligée, le passage du courant dans l['induit n'y 
entraine aucune chute de tension et la force électro- 
motrice totale de chaque phase coincide avec la force 
électromotrice d'inducteur qui y est engendrée, 

A la suite de la rotation de I'inducteur, tous les 2 x/n, 
ly a inversion de la polarité de la force électromotrice 
pour deux phases diamétralement opposées, Ces deux 
phases doivent changer de voie, ce qui signifie l'inversion 
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la voie positive à la voie négative. Compte tenu de 
origine qui a été choisie, on peut écrire durant ['empiê- 
tement le systêéme d'équation suivant (fig. 3): 


di di 
E: ) Vig [2+1) 
Meg sin O = EM, + OM, cos (7) ag 
dê dê 
(3) 
di di 
) Ka 
OMrig sin (0— x) = OM, cos (1) —— +OM, ee 
dê do 
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FIG. 3 


Au début de la commutation, soit en 0=0, les 
courants I, et tn/241) valent respectivement = JR et 
+lê. A la fin, soit en 0=%, leurs valeurs sont inver- 
sées. L'intégration du systême (3) conduit aux expres- 
sions suivantes: 


Mir 
u=——(I—cos 8) — 1,,/2 
2M, (4) 


— 


El | 
(n/2+1) 2M (I—cos 0) + 1.4/2 


L'angle d'empiéêtement se deduit immédiatement de 
[une ou l'autre des équations (4) et des valeurs finales 
du courant dans les phases commutantes: 


2M ly 
Y=arcos ( e) (5) 
| Meir 
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Durant l'empiêtement, les forces électromotrices tota- 
les aux bornes des phases non commutantes, soient e'p, 
s'expriment par: 


. di 
eK=ex—M,f cos [2(K—1) x/n] 


— M,É cos [2 [K—(n/2+1)] x/] -— (6) 
d 


Compte tenu de (1) et (4), I'équation (6) peut encore 
se mettre sous la forme: 


ep=— Q Miig sin [2 (K—1) x/n] cos é (7) 


Comme les phases entrent en commutation au pas- 
sage par zéro de leurs forces électromotrices totales, 
cette équation montre que la commutation des phases 1 
et (n/24+1) déforme la tension aux bornes des autres 
phases en déplaçant leur passage par zéro de maniêre 
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telle que leur entrée en commutation est impossible tant 
que l'angle de fin d'empiétement est inférieur à 7/2. 
Pour le mode de fonctionnement et les nombres de 
phases considérés, cette condition est verifiée puisque Y 
reste inférieur à 27/n. 

La tension de sortie du pont redresseur soit U.,, 
présente une ondulation rédiduelle de 27%/n. Sur l'inter- 
valle considéré, elle vaut: 


—de 0 à Yy 
n/2 (8) 
Un=— Lie 
2 
de Y à 27/n 
n/2+1 (9) 


Un= = dus ex 
2 


La valeur moyenne, Kas de la tension de sortie 
sobtient en intégrant les équation (8) et (9) sur les 
intervalles correspondants. Elle est égale à: 


gr n ÔMAi 
Ucp= edaimíiis ( 
T 


i+cosY ) 
2 (10) 


En éliminant l'angle d'empiétement Y entre les équa- 
tions (5) et (10), on obtient I'équation de la caractéris- 
tique externe durant le premier mode de fonctionnement. 
C'est un segment de droite d'équation (fig. 4): 

Ugo PME É dt (11) 
T r 


nôMp 
x 


FIG. 4 
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Pour déterminer le courant dans les diodes raccor- 
dées à un sommet du polygone d'induit ,il suffit de faire 
la différence des courants circulant dans les phases qui 
y sont raccordées, Lorsque cette différence est positive, 
il y a passage d'un courant dans la diode à cathode 
commune. Lorsque cette différence est négative, il y a 
passage d'un courant dans la diode à anode commune. 

Durant le premier mode, il y a successivement deux 
puis une diode en conduction à chaque borne de sortie 
du pont (fig. 2), chaque diode étant passante pendant un 
intervalle 2 %/n+Y (fig 3). Le courant moyen dans chaque 
diode vaut 1.,/n tandis que le courant de crête atteint |,p. 

A la fin du mode 1, lorque Y atteint 2 x/n, le courant 
et la tension de sortie valent respectivement: 


Meir 
lg = ———— sin? x/n (12) 
M, 
n ÔM is 
U= ———— 608º x/jn (13) 
K 


2.2. Fonctionnement jusqu'au court-circuit: 
modes 2 et 3 


Lorsque le courant de charge dépasse la valeur 
donnée par |'équation (12), il apparait un deuxiême mode 
de commutation, à angle d'empiêtement constant et égal 
à 27%=/n, durant lequel l'angle d'entrée en commutation 
de chaque phase se décale progressivement en arriére 
d'un angle f, par rapport à I'angle de passage de la 
phase par l'axe neutre d'inducteur (figs. 5 et 6). 

Ce mode est la conséquence du fait que la commu- 
tation de deux phases diamétralement opposées empéche 
l'entrée en commutation des autres phases pourvu que 
l'angle de fin d'empiêtement ne dépasse pas 7/2, c'est- 
-à-dire, tant que B+Yy < 7/2 [cf équation (7)]. L'angle B 
croit d'une valeur nulle à une valeur maximum, By, 
égale à: 


Bu = /2—-27/n (14) 


Durant le mode 2, il y a à tout instant une phase 
en commutation et deux diodes en conduction à chaque 
borne de sortie. 

Pour écrire les équations relatives à ce fonctionne- 
ment, considérons à nouveau l'intervalle débutant avec 
l'entrée en commutation des phases 1 et (n/2+1). Les 
équations qui régissent I'évolution des courants dans 
ces phases sont les mêmes que les équations (3). Toute- 
fois, elles sont maintenant valables non plus sur I'inter- 
valle [0,7] mais sur l'intervalle [B,8+27%=/n]. On peut 
écrire aprês intégration: 


Meir 
2M, 


j=— = 
1 (n/2+1) 


(cos B— cos 8) — |.,/2 (15) 


De cette équation on peut extraire la liaison existant 
entre le courant |, et I'angle de retard B, soit: 
Mri 
ly= É sin T/n sin (B+7%/n) (16) 


o 
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FIG. 6 
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La tension moyenne de sortie vaut: 


8+2 x/n 

= - n/2 
Un= a e EM ek dê 

2% 2 


ce qui donne, aprês intégration, 


a OM, 
Un= Rin 8 cos t/n cos (B+7-/n) (17) 
Tr 


L'élimination de 8 entre (16) et (17) fournit I'équation 
de la caractéristique externe durant le mode 2 (fig. 4): 


— nôMa M.l , 
Uan= Ea À cos (7/n) 1— a (18) 
T W Meg sin T/n 


c'est un tronçon d'ellipse qui joint le point: 


Mit io - Meir a 
[lh= sin? (x/n); Un= cos” (x/n)] 
My 
pour 8=0, au point: 
Mei — n Mai 
[lh = sin (2 x/n); Un= Rito 8 O (2 =/n)] 
2% 


q 


pour B= /2—27%/n (fig. 4). 


Durant ce mode, chaque diode est passante pen- 
dant 47%/n et conduit un courant moyen |.,/n avec un 
courant de crête 1.p. 
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Au delàã, apparait un troisiéme mode, pour lequel 
il y a successivement deux, puis une phase en commu- 
tation à chaque borne de sortie (figs. 7 et 8). Durant le 
mode 3, l'angle d'entrée en commutation reste fixé 
à 7/22 x/n en retard par rapport à I'axe neutre d'induc- 
teur, tandis que l'angle d'empiétement se remet à croitre 
jusqu'á atteindre 47/n. Comme les forces électromitrices 
totales aux bornes des autres phases sont toutes nulles 
lorsqu'il y a deux phases en commutation à chaque borne 
de sortie, la fin du mode 3 correspond au court-circuit 
permanent de la génératrice. 


Pour mettre en équations le fonctionnement de la 
genératrice durant le mode 3, considérons à nouveau 
instant d'entrée en commutation des phases 1 et 
(n/2+1). Avec l'origine choisie, cette entrée se fait en 
T/2—-2 7/n. A ce moment les phases n et n/2 sont encore 
en empiétement. Il convient de distinguer 3 intervalles 


(fig. 7): 


— de T/2-27%/n à T/2+Y—47/n 


Durant cet intervalle, il y a commutation simultanée 
des phases 1 et n à la borne de sortie positive et des 
phases n/2 et (n/2+1) à la borne de sortie négative. 
On peut écrire, en tenant compte de la symétrie existant 
eentre les courants dans les phases diamétralement 
opposées: 


di di Mai 
Er 
a Ps” rc veebrêsio AD 
do dê 2M, sin(2 x/n) 
di di Mi 
Ff 
E aiii 1a Ts 0040 (20) 
dê de 2M, sin(2 x/n) 


[7 


FIG. B 


— de m/2+y—4 T/n à T/2 


La commutation des phases n et (n/2) est achevée, 
celle des phases 2 et (n/2+2) na pas encore débutée. 
Durant cet intervalle, seules les phases 1 et (n/2+1) 
empiétent. Les équations qui déterminent I'évolution du 
courant dans ces phases sont les équations (3). On a donc 
di, di(n/241) Mae 


do dê 2M, 
— de m/2 à n/2-2 T/n+Y 


sin O (21) 


C'est la fin de I'empiétement des phases 1 e (n/2+1) 
et le début de la commutation des phases 2 et (n/2+42). 
On péut écrire 


di di MA 

n/2+1 Ff 
ga O : ra =" cos(6—2r/n) (22) 
dê dê 2 M,sin(2 x/n) 
di di Mi 

2 (n/2+2) Pt 
PE 7 PP (23) 
dê do 2 M, sin(2 x/n) 
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En intégrant les équations (19), (21) et (22) sur 
leurs intervalles respectifs et en tenant compte de la 
valeur initiale et de la valeur finale du courant dans la 
phase 1, on peut écrire la relation existant entre l'angle 
d'empietement 7 et le courant de charge |,,, soit: 


Meir 
ln= ——— tg (x/n) [1 + cos (4 x/n—*+)] (24) 
2M 


[1] 


La tension moyenne de sortie vaut durant ce mode 


a n Ú Meir 


= -cotg (7/n) [1-—cos (4 1/n—)] (25) 
2% 


En éliminant y entre (24) et (25), on obtient |'équation 
de la caractéristique externe durant le mode 3 (fig. 4): 


O Mir Mola 
Un=-———— cotg (x/n) [|——— cotg (r/n)] (26) 
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Durant ce mode, il y a successivement 3 puis 2 dio- 
des en conduction à chaque borne de sortie. Chaque 
diode est passante pendant un intervalle 27/n+Y. Le 
courant moyen dans chaque diode vaut toujours lp/n 
tandis que sa valeur de crête, soit lym, atteint seulement: 


— Bi Y <37T/n 

sas Meir [sina fi Rs, (27) 
2M, sin (2 x/n) 

—si Y > 3T/n 

EG [sin pa, (28) 
2M, sin (2 x/n) 


Au court-circuit l'angle de conduction des diodes atteint 
6 r/n. 


3. INTERPRÉTATION DU FONCTIONNEMENT AU MOYEN 
D'UNE MACHINE À COURANT CONTINU ÉQUIVA- 
LENTE À LA GENEÉRATRICE 


La caractéristique externe obtenue (fig. 4) se com- 
pose d'un segment d'ellipse auquel se raccordent deux 
segments de droites qui lui sont tangents aux points de 
raccord. Cette caractéristique est, pour n=6, identique 
à celle du pont triphasé alimenté par des sources sans 
couplage mutuel, Ceci résulte du fait que, pour l'applica- 
tion considérée, une source triphasée d'impédance interne 
purement inductive et sans couplages peut se ramener 
à une source hexaphasée à couplage sinusoidal par- 
fait [4]. 

Le comportement de la génératrice n-phasée à 
couplages sinusoidaux qui vient d'être étudiée peut s'inter- 
préter facilement en termes d'une machine à courant 
continu ayant n lames de collecteur, de largeur 27%/n, 
et une paire de balais dont la largeur et la position sont 
fonction du courant de charge débité. La position des 
balais, définie comme celle de leur point milieu et leur 
largeur se déduisent de la position et de la largeur de 
Vintervalle de conduction des diodes. 


Durant le mode 1, on assiste à un élargissement 
progressif des balais depuis une largeur nulle jusqu'à 
une largeur 27/n. L'origine des balais restant calée à la 
ligne neutre d'inducteur, leur déplacement par rapport 
à cette ligne neutre est égal à la moitié de leur largeur 
(fig. 9). 
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FIG. 9 


Durant le mode 2, les balais gardent une largeur 
constante et égale à 27/n, tandis que leur position passe 
de x/n à m/2-—x/n du début à la fin du mode (fig. 10). 
Durant ce mode, pour chaque valeur de courant, les 
balais se placent dans une position telle que I'inducteur 
engendre dans la section en commutation une force 
électromotrice de commutation qui correspond au courant 
à commuter. C'est lã un mécanisme de commutation 
analogue à celui qui a été réalisé dans les machines à 
courant continu sans pôles de commutation en décalant 
les balais par rapport à l'axe neutre, 


FIG. 10 


Durant le mode 3, on assiste à un nouvel élargisse- 
ment des balais dont I'origine reste fixée en 7/2—2'x/n 
(fig. 11), A la fin du mode 3, les balais ont une largeur 
de 47/n et une position égale à 7/2. C'est le court- 
-circuit permanent de la machine, Deux constatations 
peuvent donc être faites lorsque le nombre de phase 
augmente. 
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FIG. 11 


— En premier lieu, cette augmentation nentraine pas 
pour la structure considérée, |'apparition de modes 
supérieurs au mode 3 et cela contrairement au 
cas du redressement sur sources inductives non 
couplées [5]. Pour ce dernier, on distingue n/2 
modes successifs, le nombre de phases en com- 
mutation à chaque borne de sortie croissant d une 
unitê quand on passe d'un mode au suivant et 
atteignant n/2 au court-circuit, 


— En second lieu, laugmentation du nombre n de 
phases va de pair avec un élargissement de la 
zone correspondant au mode 2, La disparition 
progressive des modes 1 et 3 saccompagne dune 
diminution de langle d'empiêtement qui reste 
constant et égal à 47%/n sur pratiquement toute 
la caractéristique. Ceci entraine un écart croissant 
entre le courant de crête et le courant moyen 
vehiculé par les diodes, À la limite, pour n tendant 
vers linfini, le segment d'ellipse constitue la 
caractéristique tout entigre, tandis que l'angle de 
conduction des diodes devient infinitésimal, 


La figure 12 montre clairement que si le nombre de 
phases est suffisant, on peut remplacer avec une bonne 
approximation la caractéristique réelle par la caractéris- 
tique obtenue en prolongeant l'ellipse du mode 2 jusqu'á 
ses intersections avec les axes. Ceci revient à considérer 
que les équations (16) et (17) décrivent le comportement 
de la machine depuis la marche à vide jusqu'au court- 


“circuit, Au niveau de la machine à courant continu 
éêquivalente ceci revient à considérer qu'elle posséde 
depuis la marche à vide jusqu au court-circuit des balais 
de largeur 27/n et que leur position passe progressi- 
vement de l'axe neutre à laxe direct d'inducteur Par 
pasage à la limite pour n tendant vers [infini, on retrouve 
ainsi, en régime, le modéle de la genératrice à couplages 
sinusoidaux établ; à partir des équations des circuits à 
commutation [6]. 
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FIG. 12 


4 NVERIFICATIONS EXPÉERIMENTALES 


Des essais ont éte effectuês sur une machine hexa- 
phasée à póôles lisses et à rotor feuilleté. La f.e.m, 
induite dans chaque phase par l'inducteur est pratique- 
ment sinusoidale (fig. 13a) tandis que les mutuelles 
entre les différentes phases de l'induit, s'êcartent jusqu'á 
20 % des valeurs correspondant à un couplage sinu- 
soidal parfait (fig. 13b). 
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La figure 14 permet de comparer la caractéristique 
externe expérimentale à la caractéristique externe théo- 
rique. L'écart qui apparait entre les deux courbes aux 
environs du court-circuit, provient du défaut de couplage 
entre les phases de l'induit. Cette imperfection de 
couplage entrainant une action plus faible des phases 
en commutation sur les autres phases, il sen suit que 
les zones correspondant aux modes 2 et 3 sont partielle- 
ment réduites. La fin du mode 3 ne coincide plus avec 
le court-circuit de la machine [4]. Un quatriême mode de 
commutation apparait qui correspond au dernier tronçon 
(de pente plus faible) de la caractéristique expérimentale. 
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FIG. 14 


Les figures 15 et 16 permettent de comparer les 
formes d'ondes théoriques et expérimentales pour un 
point de fonctionnement qui correspond théoriquement au 


milieu du mode 2 (8 égal 15º), On remarque une bonne 
concordance entre les diférentes courbes malgré quel- 
ques écarts dus aux défauts de couplage et au lissage 
imparfait du courant. 
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FIG. 15 
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FIG. 16 


5. RÉFÉRENCES 


[1] P. W. Franklin: Theory of the three phase salient 
pole type generator with bridge rectified output. 
IEEE Winter Meeting, New York, Janvier 1972. 


[2] T.W. Stuart & M. W. Tripp: Steady state analysis 
of the rectified superconducting alternator. IEEE, 
AES, Vol. 13 n.º 4. 


[3] J. Debernardi & R. Gillet: Systême d'excitation 
compact pour alternateur nucléaire à 4 póles. AIM 


Liêge Centrales 1974, pp. 75; 1-9. 


[4] FF. Labrique: Analyse de la commutation dans les 
ponts redresseurs alimentés par des sources 


TÉCNICA 454 


polyphasées couplées. Université Catholique de 
Louvain, séminaire d'électronique de puissance, 


mars 1976 (note interne). 


[5] G. Seguier: 
1970. 


Les montages redresseurs. Dunod, 


[6] M. S. Garrido: Etude théorique et expérimentale 
de la génératrice à redresseur. |. Equations de 
la génératrice à redresseur. Etude du régime 


permanent. Revue E, Vol. Vil, n.º 3, 


| MACAE, "ala k 
am PR A 
MAMMA. 


” 
= 
= 


u 
1 
a 


um 


RE 


Ei ê E CENA EN, ENC at 
: aa DR . 
. tt a plata sê A LL E z E 
a , Pa RAN 
aaa nn E x adetatado state lada cara a “ na a ss. 
CEE E E a SE 


PARA: 


ECONOMIZAR ENERGIA 
COMBATER A POLUIÇÃO 
POUPAR TEMPO 
DEVOLVER AO PEÃO OS ESPAÇOS DE QUE CARECE 


METROPOLITANO 


UMA SOLUÇÃO EFICAZ 
| NO DOMINIO DOS TRANSPORTES . 
COLECTIVOS URBANOS 


METROPOLITANO DE LISBOA aaa esa 
TECNICA V 


NUMERO 454 NOVEMBRO 1979 


ANO LIV 


VOLUME XLI 


Etude de la commutation dans une génératrice 
polyphasée à redresseur 


Deuxiême partie 


Etude de la génératrice à couplages triangulaires 


RESUMO 


Este artigo é a segunda parte dum estudo sobre a 
comutação dum gerador polifásico sem colector, A partir 
de dois casos limite, o triangular aqui apresentado e o 
sinusoidal já apresentado na primeira parte, mostra-se que 
a lei de variação das inductâncias mútuas das espiras em 
função da posição pode afectar consideravelmente o pro- 
cesso de comutação e, por conseguinte, a caracteristica 
externa do gerador, 


A análise teórica é confirmada por resultados experi- 
mentais obtidos com máquinas de laboratório. 


INTRODUCTION 


Dans la premiére partie de cet article on a étudié 
le fonctionnement d'une génératrice polyphasée à redres- 
seur présentant un couplage sinusoidal entre I'inducteur 
et les phases de linduit et un couplage sinusoidal par- 
fait entre les phases de I'induit [1]. La seconde partie 
est consacrée à l'analyse du comportement d'une géné- 
ratrice présentant un couplage triangulaire entre |'induc- 
teur et les phases de [induit, et un couplage triangulaire 
parfait entre les phases de [induit. 


Rappelons que, par couplage parfait, on entend un 
couplage dont la loi est telle que la mutuelle entre deux 
phases occupant la même position géométrique est égale, 
en valeur absolue, à la self propre de chacune de ces 
phases. 


On considêre le cas d'une machine à pôles lisses, 
sans amortisseurs, à induit polygonal et nombre pair de 
phases. Le nombre n de phases doit être supérieur ou 
égal à 6 pour que la présente étude s'applique sans 
aucune restriction. 


Comme pour la machine à couplages sinusoidaux, 
I'étude effectuée dans le cas du régime permanent, conduit 
à une solution analytique du comportement depuis la 
marche à vide jusqu'au court-circuit, moyennant les hypo- 
thêses et approximations suivantes: 
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Premier Assistant 
Université Catholique de Louvain 


ABSTRACT 


This paper is the second part of a study on the 
commutation processes in a polyphase commutatorless 
DC generator. From the study of two limit cases, the 
triangular one presented here and the sinusoidal one, 
yet presented in the first part, it is shown that the form 
of the variation of the mutual inductances as a function 
of the relative positions of the different coils can affect 
considerakiy the commutation processes and hence the 
generator output characteristic, 


Experimental results on laboratory test machines 
confirm the theoretical analysis, 


— |'inducteur est alimenté par une source de courant 

— le courant de sortie du pont redresseur est parfai- 
tement «lissé» 

—la résistance interne des phases de l|'induit est 
négligeable 

— les chutes directes et courants inverses des dio- 
des sont négligeables, 


1. EQUATIONS DE DÉPART POUR LA MACHINE 
À COUPLAGES TRIANGULAIRES 


Avec les notations de la figure 1 et une source de 
courant dans I'inducteur, les forces électromotrices engen- 
drées par ['inducteur, dans les phases de lI'induit, s'écri- 
vent, en prenant comme origine le passage de la premiêre 
phase par l'axe neutre relatif à |'inducteur 


ex = ÔMirrect[0-2(K—1)m/)n] 1<K<n (1) 


Dans cette équation Mg représente la valeur maximum 
de la dérivée de la mutuelle entre l'inducteur et une 


phase de I'induit, O la position de I'induit, 0 sa vitesse 


Travail fait en collaboration avec Centro de Automática da Universidade de Lisboa (l. N.1.C.) 


Manuscrito recebido para publicação em 10/7/78. 
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| dxe nevére 


FIG. 1 


et la fonction rect [8-2 (K— 1) x/n] une fonction rectan- 
gulaire antisymétrique, de période 2 7, d'amplitude unitaire 
et décalée de 2 (K—1)7/n par rapport à I'origine (fig. 2). 
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FIG. 2 


Les mutuelles entre les phases de l'induit sont de 
la forme: 
Mij=M, tri [2 (tj) m/n] + Cl ISA (2) 

Dans ces équations M, représente la valeur maximum 
de la mutuelle entre phases et la fonction tri (1) une 
fonction triangulaire, symétrique, d'amplitude unitaire et 
de période 27%, (fig. 3). Les expressions analytiques 
des fonctions rect [0-2 (K-1)m/n] et tri (4) sont don- 
nées à lannexe 1, 
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FIG. 3 


2. ANALYSE DU FONCTIONNEMENT 


2.1. Fonctionnement à faible charge: modes 1 et 2 


Le premier mode de fonctionnement rencontré dans 
la machine à couplages triangulaires correspond à un 
mecanisme de commutation identique à celui du mode 1 
de la machine à couplages sinusoidaux (cfr le para- 
graphe 2.1. de la premiére partie de cet article [1)). 

A tout instant, linduit se partage en deux voies 
mises en parallêle, chacune parcourue par la moitié du 
courant de charge, l'une comportant toutes les phases 
dont les f.e. m. sont positives et l'autre toutes les 
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phases dont les f.e.m. sont négatives (fig. 4.a). La borne 
positive de la charge-est ainsi reliée au sommet du 
polygone d'induit dont le potentiel est le plus élevé et 
la borne négative de sorte au sommet du polygone 
d'induit dont le potentiel est le plus bas. Tous les 2K r/n, 
par suite de la rotation de l'induit, deux phases diamé- 
tralement opposées changent de voie, ce qui entraine 


un court-circuit temporaire de ces phases, pendant que 
le courant s'y inverse (fig. 4.b). 


Afin d'écrire les équations relatives au fonctionne- 
ment de la génératrice à couplages triangulaires durant le 
mode 1, considérons lI'intervalle qui débute par I'entrée en 
commutation des phases 1 et (n/2+1) (fig. 5). Durant 
'empiétement, on a le systême d'équations suivant: 


| o do di 
0Mi=0M, ET E ii 
dO dO (3) 
CC do di 
DMA = aii, 
do dO 


Aprês intégration, compte tenu des valeurs initiales 
et finales du courant dans les phases considérées, on 
obtient I'équation de liaison entre I'angle d'empiêtement Y 
et le courant lp: 


Meir 
2M 


“ 


Durant lempiétement, les forces électromotrices totales 
aux bornes des phases non commutantes valent, en 
tenant compte de (1) et (3) 


ex= U Meir rect [0—2 (K— 1) x/n] — 
— 0 Migtri [2(K— 1) x/n] (5) 


TÉCNICA 454 


Les deuxiêmes termes de ces équations correspondent 
aux tensions induites par la variation de courant dans 
les phases en empiétement Comme ces tensions sont, 
en valeur absolue, strictement inférieures à 0 Mçi,, elles 
ne peuvent pas déplacer les instants de passage par 
zéro des forces électromotrices totales dans ces phases 
par rapport aux instants de passage par zéro des forces 


Lx 
c) Y< B< — 


électromotrices, induites par I'inducteur, correspondantes. 
En conséquence, si, comme pour la machine à couplages 
sinusoidaux, l'empiêtement de deux phases diamétrale- 
ment opposées n'entraine I'entrée en commutation d'aucune 
autre phase, cet empiétement ne peut toutefois pas 
déplacer les positions d'entrée en commutation des autres 
phases. Celles-ci restent fixées au passage de ces phases 
par l'axe neutre relatif à I'inducteur, 


La tension de sortie, qui présente, comme dans le 
cas sinusoidal, une ondulation résiduelle de période 2 7/n, 
a une valeur moyenne égale à: 


CS n 0 Meir 
Uch = e (1—ry/7) (6) 


L'élimination de Y entre (4) et (6) fournit I'équation 
de la caractéritique externe durant le mode 1: 


— nÔMa nOM,. 
FR a a (7) 
2 Tr 


C'est un segment de droite qui joint le point: 


-— n 0 Meir 
(Un= o ln axe 0) 
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Es 


2x 


au point (fig. 6): 


Meir = 


Mo n 


[Ua= E ad Ô Meir b lh= 


nôMy if 


Mp ro 


EM, 


FIG. 6 
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FIG. 5 


Durant ce premier mode, il y a successivement deux 
puis une diode en conduction à chaque borne de sortie. 


2% 
Chaque diode est passante pendant —— -+ Y et conduit un 
n 


courant moyen qui vaut |4/n avec un courant de crête 
égal à lp (fig. 5). 

Lorsque le courant de charge dépasse la valeur 
obtenue en remplaçant Y par 27%/n dans I'équation (4), 
un deuxiême mode de commutation apparait, pendant 
lequel lI'angle d'empiétement croit de 27/n à 47/n. 
Durant ce mode, la fin de l'empiêtement d'une phase 
chevauche le début de l'empiétement de la phase suivante. 
Ceci résulte du fait que l'empiêtement d'une phase à 
chaque borne de sortie (comme signalé plus haut) ou 
de deux phases à chaque borne de sortie (comme on le 
verra plus loin) ne peut déplacer les instants d'entrée en 
commutation des autres phases. Ceux-ci restent fixés au 
passage de ces phases par I'axe neutre relatif à inducteur. 

Pour écrire les équations relatives au mode 2, consi- 
dérons, à nouveau, I'instant d'entrée en commutation des 
phases 1 et (n/2+1). L'empiêtement de ces phases se 
divise en trois intervalles (figs 7 et 8): 


— de 0 à y-27%x/n 


| y a commutation simultanée des phases 1 et n à 
la borne positive de sortie et des phases n/2 et (n/241) 
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; x 
a) ocO<cp.2% b)y-2T <0<rm 
FIG. 7 
eq 
8 
K Rx 
eu 
o) 
x LX 


FIG. 8 
à la borne négative. En tenant compte de la symétrie Pendant cet intervalle, les forces électromotrices dans 
existant entre les courants dans les phases diamétrale- les phases non commutantes valent: 
ment opposées, on obtient: 
: ; di . di E a QUAL Rail a E 
0Mir=20M, —— + 20M, (1—4/n) — (8) e'x= 0 Mir rect [0-2 (K—1) x/n] 
dO do (9) 
. . di d . di . di 
0 Mi=20M, (1-4/n) — + 2M, a —20M, tri [2(K—1) m/n] - —20M,tri [2Kz/n] — 
dê dê de dl 
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En tenant compte de (8), I'équation (9) peut encore 
sécrire: 


ek =0 Mig rect [0—2 (K—1) m/n ] — 
sl Mai tri [2 (K—1) m/n]+tri (2 K x/n) (10) 
1+tri (2 x/n) 


Cette equation appelle les mêmes commentaires que 
léquation (5). | ny a pas de déplacement des instants 
de passage par zéro des forces électromotrices totales 
par rapport aux instants de passage par zéro des forces 
électromotrices correspondantes induites par [inducteur. 


— de y—-27%/n à 27/n 


Seules les phases 1 et (n/2+1) empiétent. Les équa- 
tions (3) et (5) sont applicables durant cet intervalle. 


— de 2%/n à Y 


Il y a commutation simultanée des phases 1 et 2 à 
la borne positive et des phases (n/2+1) et (n/2+2) à 
la borne négative. On peut écrire: 


4 
o 
- 


di 
+ 20M, (1—4/n) — 


0 Mi=28 My 
do do (11) 
di, di 
OMir=20M, (1—4/n)—— + 20M, — 
dê do 


Les tensions aux bornes des phases non commu- 
tantes valent: 


e'«=0 Mir rect [0—2 (K—1) m/n] — 
6 Ma tri [2 (K—2) x/n] +tri [2 (K—1) x/n] (12) 
1 +tri (2 m/n) 


Il convient de remarquer que les équations (11) et 
(12) peuvent se déduire des équations (8) et (10) en y 
incrémentant d'une unité, modulo n, tous les indices des 
phases en commutation. 


Pour déterminer |'évolution du courant dans la phase 1 
durant son empiétement et déterminer l'équation liant Y 
à lp, On résoud les équations (8), (5) et (11) par rapport 

di, 

dê 
dité respectifs. Compte tenu de la valeur de i, en début 
et en fin d'empiétement, on obtient facilement: 


E 


, puis on les intêégre sur leurs intervalles de vali- 


Mei —2%/n 
E E [ 


À + (47/n—9)] (13) 
2 Mo | —2/ n 
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La tension moyenne de sortie vaut, durant le mode 2: 


Y—2 T/n 


Y n Eis furo 
U.p= fe — 1 ) Ô Meig— O Meir ( E nr ) (1 4) 


L'élimination de Y fourmit I'équation du segment de 
droite que constitue la caractéristique externe durant ce 
mode, soit: 


a (n—2) À Molch 
Usn=(n—2)0 Mis (15) 
2T/n 


Durant le mode 2, chaque diode est passante pendant 
2rin+Y. Il y a successivement 2 puis 3 diodes en 
conduction à chaque borne de sortie, Le courant moyen 
dans chaque diode vaut |,/n, tandis que le courant de 
crête, Im, y atteint seulement: 


Meir Y—2 T/n 


ia RR (16) 
2M, 2(1-2/n) 


Im e lg— 


A la fin du mode 2, I'angle de conduction des diodes 
atteint 67%/n et le courant de crête qui les traverse 
vaut 1/2. 


2.2. Fonctionnement jusqu'au court circuit: 
modes 3 à n/2 


Au delã du mode 2, apparaissent une série de modes, 
correspondant chacun à un élargissement de 27%/n de 
'angle d'empiêtement Y, et à l'augmentation d'une unité 
du nombre de phases en commutation à chaque borne 
de sortie. Ainsi durant le mode K(3 <K<n/2), ly a 
succcessivement K puis (K—1) phases en commutation 
à chaque borne de sortie, chaque phase entrant en 
commutation lors de son passage par l'axe neutre 
(cfr. annexe 2). Toutefois, lélargissement de l'angle 
d'empiêtement n'entraine pas, pour la machine considérée, 
un élargissement correspondant de l'angle de conduction 
des diodes, celui-ci ne dépassant en aucun cas 87/n par 
période, Afin de mettre en évidence ce phénomêne, peu 
prévisible à premiére vue, nous allons considérer un 
mode K quelconque, Il convient, pour ce mode, de 
distinguer pour I'empiétement d'une phase (2K-—1) inter- 
valles successifs. lls sont constitués d'une alternance 
dintervalles ou il y a K phases en empiétement et 
d'intervalles ou il y a (K—1) phases en empiêtement 
(fig. 9). 

Si Ion considêre I'empiêtement de la phase 1, on 
peut écrire, pour les deux premiers intervalles (fig. 10). 


— de 0 à Yy-2(K-1)zx/n 


ly a K phases en commutation à chaque borne de 
sortie, les courants dans les phases en commutation à la 
borne positive sont régis par le systêéme suivant d'équa- 
tions: 
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Mtis 2M, 2Mtri(2mn) 2M,tri [2 (K— 1) x/n] E 
Meir 2M, tri (2 T/n) 2M, PRA TEPAR ES SPA PEL AEAeA DORES papos GEES RS d in à 
= e 1 

O AP TS O 
Mais 2M.tri [B(K— mn] 2M tri [2(K=D) mn) 2M, e 

8 

= = j 

1 


| 


al oc0<p2(k-4)7 b) p-tka4)Dçô<cta 


FIG. 10 


Les forces électromotrices totales aux bornes des  e;=0Mrrect [0-2 (j—1) m/n]— 


phases non commutantes valent, durant cet intervalle Dida tri (2 4-1) m/n]+tri [20 +K—2) m/n] (18) 


(cfr. annexe 2): 1 +tri [2 (K— 1) m/n] 
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— de Y—-2(K—1)r/n à 27/n 


ly a (K—1) phases en commutation à chaque borne 
de sortie. Les courants dans les phases en commutation 
à la borne positive sont régis par le systéme suivant: 


Meir 2M, 2M,tri (2 x/n) 
Meir 2M,tri (2 m/n) 2M, 
Meir 2M,tri [2 (K—2) T/n] 2M,tri [2 (K—3) T/n] 


Les forces électromotrices totales aux bornes des 
phases non commutantes valent durant cet intervalle 
(cfr.annexe 2): 


e; =0 Mir rect [8—2 (;—1) x/n]— 
o Mi tri [2 (1-1) m/n] +tri [2 ((+K—3) x/n] (20) 
1 +tri [2(K—2) m/n] 


Pour les (2K—3) intervalles restants, les équations 
obtenues sont analogues et peuvent se déduire des 
équations (17) à (20) moyennant une adaptation des indi- 
ces des phases en commutation, Ces indices doivent 
être augmentés, tous les 27%/n, d'une unité modulo n. 


L'évolution du courant i, dans la phase 1 durant 
son empiétement s'obtient en inversant les équations du 


di 

type (17) et (19) et en intégrant chaque — obtenu sur 
de 

Hintervalle qui lui correspond, L'inversion des équa- 


tions (17) ou (19) conduit à (cfr. annexe 2) 


— pour I'équation (17): 


di, dns Melk 

de de 2M, [1 +tri [2 (K—1)rx/n]] 

di (21) 
di, o dia) = ' (n=K+3) 

do dO TO d8 
— pour I'équation (19): 

di damas o Miele 

do dO 2M,[1+tri[2(K—2)m/n]] 

| (22) 
di, o din) E ditn-K+4) cê 

dO dO TO dO 


Les équations (21) et (22) et les équations analo- 
gues, qu'on obtient pour les intervalles suivants, montrent 
que linversion du courant dans la phase 1 se fait en 
trois étapes, La premiére, de O à 27/n, fait passer le 


Uu= (n/2—K) 0 Meir ; lh= 


courant dans la phase est nul, La troisiême de Yy—2%/n 
à Y est symétrique de la premiére et fait passer le 
courant de O à 1|.w/2 (fig. 10). On voit ainsi clairement 
que les modes 3 à n/2 correspondent tous à un même 
mécanisme de commutation. 


aca 2M,tri [2 (K—2) x/n] ii 
O E d lin = 
q ER =" 
ne 0 |. ( 
E 2M, Ea 


courant dans la phase de —l.,/2 à O. La seconde, de 
27r/n à Y-—27=/n, correspond à un intervalle ou le 

Le courant dans les deux diodes raccordées à un 
même sommet du polygone d'induit se calcule, comme 
pour la machine à couplages sinusoidaux, en faisant la 
différence des courants dans les deux phases qui y 
aboutissent. 

Sauf pour le mode 3, que nous examinerons ensuite, 
le courant dans les diodes a la forme représentée à la 
figure 10. Chaque diode conduit pendat deux intervalles 
dinstincts, de 47%/n chacun, I'un situé de O à 47/n, et 
l'autre de Y—27%/n à Y+27%/n, Le courant moyen d'une 
diode vaut toujours |.,/n tandis que son courant de crête 
atteint seulement |.4/2. 

Durant le mode 3, l'angle de conduction des diodes 
passe progressivement de 67/n à 87%/n, chaque diode 
conduit un courant moyen |.,/n avec un courant de crête 
qui vaut |.,/2 (fig. 11). A la fin du mode 3, le courant 


Ina4 


FIG. 11 


des diodes passe par zéro au milieu de l'intervalle de 
conduction, C'est le début de la scission de lI'intervalle 
de conduction en deux intervalles distincts caractéris- 
tiques des modes supérieurs au mode 3, 


Durant les modes 3 à n/2, la caractéristique externe 
se compose d'une série de segments de droite de pentes 
décroissantes. Le segment de droite qui correspond au 
mode K joint les points suivants: 


U ) Mai, 4 x/n 

Un=(n/2—-K+1) O Meir ; e 

2M, [1 +tri [2 (K-2) x/n]] 
Mai; 4 /n 


2M, [1 +tri [2(K—1) x/n]] 
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3. INTERPRÉTATION DU FONCTIONNEMENT AU MOYEN 
D'UNE MACHINE À COURANT CONTINU ÉQUIVA- 
LENTE À LA GENÉRATRICE 


Le comportement de la génératrice n-phasée à coupla- 
ges triangulaires peut s'interpréter en termes d'une 
machine à courant continu ayant n lames de collecteur, 
de largeur 27/n et deux paires de balais dont la largeur 
et la position sont fonction du courant de charge débite. 
Comme pour la machine à couplages sinusoidaux, la 
position des balais (définie comme celle de leur point 
milieu) et leur largeur se déduisent de la position et de 
la largeur des intervalles de conduction des diodes. 

Durant le mode 1, on assiste au développement 
d'une paire de balais dont l'origine reste calée à la ligne 
neutre d'inducteur et dont la largeur croit avec la charge 


FIG, 12 


jusqu'ã atteindre 27/n à la fin du mode 1 (fig. 12). A ce 
moment la position de cette paire de balais vaut T/n. 
A partir du mode 2, cette paire de balais se dédouble: 
des deux paires de balais ainsi formés, de largeur 2 7/n 
chacune, l'une reste fixée à la position 7/n et I'autre se 
décale progressivement jusqu'ã atteindre une position T— 


FIG 13 
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— T/n au court-circuit (fig. 13). Ce fonctionnement corres- 
pond à une mise «hors service» d'une portion croissante 
de l'induit sans que pour autant il y ait un étalement de 
courant. 

Lorsque le nombre de phases croit, le mode 1 tend 
à disparaitre, A la limite, pour n tendant vers [infini, 
seuls subsistent les modes 2 à n/2, Les balais deviennent 
ponctuels et le mécanisme de commutation se rêéduit à 
la rotation de la deuxiême paire de balais. 

Il est intéressant de remarquer que ce mécanisme 
coincide avec le modéle adopté pour introduire la correc- 
tion spatiale dans l'étude de la génératrice à redres- 
seur [2]. 


4. VERIFICATIONS EXPÉRIMENTALES 


Des essais ont été effectués sur une machine dodé- 
caphasée à pôles lisses et à rotor feuilleté. La fem. 
induite dans chaque phase par |'inducteur est légérement 
trapézoidale (fig. 14.a). On note toutefois la présence 
dharmoniques de denture assez importants, Les mutuel- 
les entre les phases de Iinduit sont três proches des 
valeurs correspondant à un couplage triangulaire par- 
fait (fig. 14.b), (voir page suivante), 

La figure 15 permet de comparer la caractéristique 


LI Fw) 


— Theo rique 
im Th Cort que corria EE 


— — Expúriment ale 


so 


ds 


FIG. 15 


externe expérimentale avec la caractéristique externe 
thêorique. La chute ohmigue dans linduit et la chute 
directe sur les diodes n'étant pas négligeables pour la 
machine utilisée ,on a corrigé la courbe théorique pour 
tenir compte de ces chutes de tension. L'écart qui 
subsiste entre la courbe théorique corrigée et la courbe 
expérimentale provient de la forme légérement trapé- 
zoidale des f.e.m, induites par I'inducteur, On remarque 
toutefois la bonne concordance générale entre les résul- 
tats théoriques et expérimentaux, 
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FIG. 14 
Les figures 16 et 17 permettent de comparer les CONCLUSION 


formes d'ondes théoriques et expérimentales pour un 
point de fonctionnement correspondant thêoriquement à 
la fin du mode 4. 

On remarque le bon accord entre les résultats thêo- 
riques et expérimentaux, notamment en ce qu! concerne 
l'allure des courants de diodes et des courants de phases, 


U:tôv/div; Tifa/div; E: tomsS/div 


FIG. 16 


pas (A) 


Epry 


FIG. 17 
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L'étude effectuée a permis de montrer que le com- 
portement des génératrices polyphasées à redresseurs 
dépend três fortement de la forme des lois de couplages 
existant entre les différents enroulements de la machine. 

Pour chacun des cas extrêmes étudiéês on a pu 
dégager un mécanisme fondamental de commutation. Dans 
les deux cas, ce mécanisme jouit des propriétes suivantes: 


— il s'applique depuis la marche à vide jusqu au 
court-circuit lorsque le nombre de phases tend 
vers Vinfini; 


— il correspond à des modes de commutation ou 
les intervalles de conduction des diodes ont une 
largeur constante, égale à 47/n, et en fonction 
de la charge, une position variable par rapport à 
axe neutre de linducteur, Pour la machine à 
couplages sinusoidaux il y a un intervalle de 
conduction par diode et par période. Pour la 
machine à couplages triangulaires, il y en a deux 
dês la fin du mode 3 et jusquáã la fin du mode 
n/2 


Pour des machines dont les lois de couplages s ecar- 
teraient sensiblement de l'un ou lautre des cas consi- 
dérês il faudrait envisager des mécanismes de commu- 
tation plus complexes ou les intervalles de conduction 
des diodes verraient non seulement leurs positions mais 
aussi leurs largeurs varier de maniére importante avec 
la charge. 
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ANNEXE 1 
Expression analytique des fonctions rect [)—2 (K—1) x/n] et tri (1) 


La fonction rect [0—2(K—1)7=/n] s'exprime de la maniére suivante à partir des fonctions «échelon unité» 
u(O—<) (fig. A.1.1): 


oo 
rect [0-2 (K—1)q/n] = u[8-2(K—-9)zx/n] + 25. (—1)Ku [0-2 (K— 1) m/n—Kz] 
1 


.o 
—u[-0+2(K—-1)x/n)-25.(—1)Ku [—-0+2 (K— 1) x/n— Kx] (A 1.1.) 
1 


2u/-0ralx-1)rn-rc) | rectlo - 2t8-1) 7/0] Lu(e-2(n-1)rmt/n)- 210) 


-V(-0 + 2(M-1)M/n) | -R2U(0-Uk-1)rein)-Te) 


FIG. A 1.1 


Tri(y) 
/ 
(a(yp-n)/nr) u(y-ro) 


= tr 2yru ty) 
Me 


N 
s 
Ns 
E 
n. 


FIG. A 1.2 
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La fonction tri (4) s'exprime de la maniére suivante à partir des mêmes fonctions «échelon unité> (fig, A.1.2.) 


Do 
tri) = (1-2 Y yuh-2 0 (-1k (ES a) ) u (pKa) 
[e] 
+ (14 28 yutoya tk (ERHHO suçoapka) 
| TT 1 T (A. 1.2) 


ANNEXE 2 
Résolution des équations pour les modes 3 à n/2 


Dans les modes 3 à n/2 décrits au paragraphe 2.2, il y a empiêtement simultané de j phases consécutives à 
chaque borne de sortie, chaque phase entrant en commutation à son passage par l'axe neutre. Pour le mode K, le 
nombre de phases en empiétement à chaque borne vaut successivement j=K puis j=K-—1. Les conditions d'existence 
de ces modes sont les suivantes: 


() il y a effectivement j phases consécutives en court-circuit à chaque borne de sortie. Ces court-circuits 
se faisant à travers les diodes correspondantes (ã cathode commune à la borne positive et à anode 
commune à la borne negative) et les diodes ne pouvant laisser passer des courants inverses non nuls, 
ceci revient à admettre que la valeur algébrbique du courant dans chaque phase commutante à la borne 
positive (respectivement négative) est inférieure ou égale (respectivement supérieure ou égale) à la valeur 
algéebrique du courant dans les phases entrées en commutation avant elle. Ainsi, si à la borne positive 
Hordre d'entrée en commutation des phases est "', ”,',... on doit avoir (fig. A ZA: 


E Ss E St SU Siso [12 


FIG. A. 21. 


(ii) les tensions induites dans les phases non commutantes ne déplacent pas les positions auxquelles se pro- 
duisent les changements de polarité des forces électromotrices totales induites dans les différents phases. 
En effet, le non respect de cette condition ferait que les phases n'entreraient pas en commutation à leur 
passage par l'axe neutre, 


Pour vérifier l'existence des modes 3 à n/2, il suffit de résoudre les équations correspondant aux configura- 
tions considérées et de s'assurer que la solution obtenue satisfait aux conditions d'existence citées ci-dessus. 

A un instant O, (0< 0 < 2%), les phases en commutation à la borne positive sont les phases 1, n, n—1, ..., 
n—-j+2 (|=K ou K—1). On peut écrire le systêéme d'équations suivant, compte tenu de la symétrie qui doit exister 
entre les commutations à la borne positive et à la borne négative: 


Meir 1 tri (2 x/n) ORAR)  issssmass trifl20—1) x/n] | p 
Mei tri (2 x/n) 1 til) meme tri [2 (1) m/n] | 
Bê 
Mai =2M. | tri(4x/n) tri (2 T/n 1 A usar X— 
E!F o ( ) do (14) 
Meir | [ti26-Dmn] tmBQO-)r nn] tniBG-3) Tn] a. 1 (n-1+2) 
(A. 21.) 
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Pour résoudre ce systême, en désignant par (1), (2), ..., (1), ..., ()) les lignes de la matrice des coefficients, 
on effectue les substitutions linéaires suivantes 


Os MD-12[(-9) + (1+ 1)] 2<i;€j|-1<n/2 


Compte tenu de la linéarité de la fonction tri (14), le systêéme (A. 2.1.) devient: 


Meir tri2 m/n tr),  cnseamanii ti2(]=—)r/n h 
[O] 1—tri (2 m/n) (à) (à) (0) E 
0 0 tri (2 x/n) 0 0 (n-1) 
Meir tri [2 (j — 2) T/n] tri [2 (|—3) T/n] 1 (n-j+2) 
(A. 2.2.) 
La solution de ce dernier systême est immédiate: 
e E ditn—j+3) = 
do Fr dO 
(A. 2.3.) 
do A Meir 
dê de 2 Ms [1+tri [2 (1-1) 7/n]] 


Comme les dérivées du courant dans toutes les phases supposées en commutation sont 2 0, la condition (i) 
est vérifige. En effet on peut écrire pour le courant dans une phase commutante 


0 
e epi dia 
| (0) = 1 (0,) + — dê 
dê 
Comme - est 20, la fonction i (Bb) est croissante et par conséquent 
d 


(O) = à" (0-2 m/n) <i" (0) 


Les forces électromotrices totales e, aux bornes des phases non commutantes s'expriment par 


o di; 
e1=e—2,9M a ie (tin, a n=j+2) 
d 


(A. 2.4) 
Compte tenu de (A.2.3.), ces forces électromotrices valent 
ti f2(i-yrin]+trif2(i4+=2)rin 
e, = 0Mç [rect [0-2 (1-1) x/n])— Mi si a ii 
1+tri [2 (= 1) m/n] (A. 2.5.) 


Dans ces équations, le deuxiême terme entre crochets correspond aux tensions induites par la variation du 
courant dans les phases commutantes. Comme sa valeur absolue est strictement inférieure à I'unité, la condition (ii) 
sur les instants de passage par zéro des f.é.m, totales aux bornes des phases non commutantes est également 
vérifiée. 


L'existence des modes 3 à n/2 est ainsi justifiêe. 
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Projecto apoiado em computador 
de Sistemas de Controlo 


RESUMO 


O artigo descreve detalhadamente o bloco de mani- 
pulação e tratamento de dados de um conjunto articulado 
de programas para o projecto apoiado em computador de 
sistemas de controlo. Apresentam-se os algoritmos utili- 
zados e os correspondentes fluxogramas para além de 
alguns exemplos que os ilustram. 


O. O COMPUTADOR E O PROJECTO 
DE SISTEMAS DE CONTROLO 


O estudo do controlo de sistemas tem meio século 
de existência, não incluindo trabalhos anteriores sem uma 
sistematização de controlo, muito embora sistemas con- 
trolados existam desde sempre na natureza, O próprio 
homem é um mundo de sistemas controlados; desde a 
regulação do ritmo ou volémia cardiaca até ao segui- 
mento de uma trajectória enquanto caminha, Neste último 
caso, um exemplo muito intuitivo, teremos a visão como 
sensor, o cérebro como controlador/processador de sinais 
e os músculos principalmente das pernas como actuado- 
res. Estes sistemas não têm, no entanto órgãos de coman- 
do construídos pelo próprio homem, Estes são sistemas 
onde pode ter lugar um estudo do sistema controlado e 
da acção de controlo mas onde, apenas em situações 
excepcionais, o homem tem uma intervenção directa. 

Na antiguidade podem-se detectar sistemas, inicial- 
mente em cadeia aberta, em que já se encontra a acção 
de controlo. É o caso do mecanismo de Hero utilizado 
no séc, | para abrir as portas de um templo [1]. O ar 
aquecido pelo fogo sagrado obrigava a água de um 
depósito a passar para um recipiente, cujo peso aumen- 
tava, até compensar a acção de um contrapeso e levar 
à abertura da porta, 

Em 1788 aparece o regulador de James Watt, para 
regulação de velocidade de motores, vulgarmente indicado 
como o primeiro regulador em cadeia fechada, 


ISABEL LOURTIE (1) 
FRANCISCO CARVALHAL (2) 
PEDRO LOURTIE (3) 


ABSTRACT 


The paper describes with detail the data manipulation 
and data preparation block of a control systems computer- 
-aided design package. The algoritms and corresponding 
flow charts are presented together with some application 
examples. 


Em 1922 aparece um estudo por Mimorsky em que 
são introduzidos elementos não-lineares no controlo apli- 
cado ao sistema direccional de navios, 


Só em 1934 aparece um trabalho de Hazen em que 
apresenta alguma da nomenclatura e sistematização de 
sistemas de controlo, em particular de servomecanismos. 
De então para cá há uma aceleração considerável no 
desenvolvimento destas matérias, sobretudo depois da 
|| Guerra Mundial. 

Os métodos clássicos de análise e projecto de siste- 
mas de controlo, em virtude das limitações impostas pelo 
cálculo manual, baseiam-se em aproximações, ábacos e 
outras técnicas simplificativas e sobretudo no domínio da 
frequência. Simultaneamente o desenvolvimento e redes- 
coberta de métodos matemáticos hoje empregues foi 
sendo feita com a motivação do controlo de sistemas, 
embora a sua aplicação tenha sido limitada na prática a 
aplicações simples. 

O advento do computador veio então permitir o 
rápido desenvolvimento de teorias até aí não exploradas; 
começaram então os grandes desenvolvimentos das téc- 
nicas no domínio do tempo. Os grandes investimentos 
realizados na corrida para o espaço e as aplicações 
militares, forneceram a motivação e os meios necessários 
a um rápido desenvolvimento do estudo do controlo 
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automático. No entanto, cada vez mais estes desenvol- 
vimentos vão encontrando aplicação noutros dominios. 

A necessidade de cada vez maior precisão, quali- 
dade, rendimento, leva a considerar no projecto de siste- 
mas de controlo um detalhe cada vez maior tanto na 
descrição do sistema como nas especificações de con- 
trolo. Surgem assim os métodos de tratamento de siste- 
mas multivariáveis, o controlo adaptativo, estocástico, etc. 
O desenvolvimento da modelação e os sucessos obtidos, 
aliados às facilidades de cálculo dos computadores digi- 
tais, levam a estender estes métodos aos campos econó- 
mico, biológico, social, etc. Em todo este desenvolvimento 
a vulgarização do uso do computador é pedra funda- 
mental, na medida em que permite realizar cálculos até 
então impensáveis, 

É assim que os mais variados tipos de projecto se 
apoiam cada vez mais no computador embora nalguns 
casos se limitem a efectuar os cálculos mais complexos. 
A proliferação de terminais e a utilização de sistemas 
operadores cada vez mais sofisticados e protegidos con- 
tra falhas de utilização, permitem a utilização da máquina 
mesmo por operadores menos experientes. No entanto, 
para que o seu trabalho seja proveitoso, é necessário 
que existam programas estruturados de forma a facilitar 
o seu contacto com a máquina, 

Neste sentido, aparecem os sistemas articulados de 
programas para o projecto em computador, usando téc- 
nicas conversacionais, i.e. incluindo no programa as 
perguntas e respostas que possam orientar o trabalho do 
utilizador. É com este espirito que está a ser desenvol- 
vido o trabalho que aqui se descreve e que resulta de 
uma colaboração entre o CAU.T.L e o LN.E.C. 
Este conjunto articulado de programas destina-se ao 
projecto em computador de sistemas de controlo, usando 
técnicas conversacionais e com uma estrutura expandivel 
por forma a integrar todas as técnicas de análise e 
sintese que se considerem relevantes neste dominio. 


1. MANIPULAÇÃO DE DADOS 
DE SISTEMAS UNIVARIÁVEIS 


1.1, INTRODUÇÃO 


O Projecto Apoiado em Computador de Sistemas de 
Controlo Automático — CONPAC — é constituido por dois 
grandes blocos de programação, um relativo a sistemas 
univariáveis e outro a sistemas multivariáveis. 

Cada um destes blocos realiza três funções dife- 
rentes: 


— Tratamento de Dados 
— Projecto de Sistemas de Controlo 
— Simulação no tempo 


Nos sistemas univariáveis, o tratamento de dados é 
realizado em dois blocos distintos: 


— Manipulação de Dados 


— Adaptação da função de transferência a dados 
experimentais («curve fitting») 
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Quanto ao projecto de sistemas univariáveis serão 
implementados os seguintes métodos: 


— Diagrama de Bode 
— Diagrama de Nyquist 
— Diagrama de Evans («Root Locus») 


No projecto de sistemas multivariáveis considera-se: 


— Método baseado nos Teoremas de Gershgorin e 
de Ostrowsky 

— Método de Controlo Modal (deslocamento de 
polos) 

— Método do Espaço de Estados 


1.2. SISTEMAS UNIVARIÁVEIS 


Para os sistemas univariáveis a estrutura utilizada 
tem o seguinte diagrama de blocos: 


| SUPERVISOR + 


(7) 


(1) Manipulação de dados 

(2) Sistema em anel fechado 

(3) Resposta no tempo 

(4) Diagrama de Nyquist 

(5) Diagrama de Bode 

(6) Lugar geométrico das raizes («Root Locus») 

(7) Adaptação da função de transferência a dados experimentais 
(«curve fitting») 


Do ponto de vista da programação, o «supervisor» 
é o programa principal e os blocos numerados de (1) a 
(7) são estruturados como sub-rotinas, pelo que os con- 
tactos entre eles são possiveis através do «supervisor». 
Esta estrutura, usando como linguagem de programação 
o FORTRAN, está a ser desenvolvida no computador 
DEC-10 do L.N.E.C. 

Os sistemas a tratar obedecem à nomenclatura 
apresentada no seguinte diagrama de blocos: 


FIG. 2 


G — sistema controlado 
C — controlador 
H — realimentação 
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Este bloco permite realizar qualquer transformação 


pus O Bloco de.Manipulação de dados aceita a descrição 
dos sistemas C, G e H de três formas diferentes: 
todos os outros, 


— Espaço de Estados 
— Função de transferência na forma racional 


— Função de transferência factorizada 


1.3. FLUXOGRAMA DO BLOCO 
DE MANIPULAÇÃO DE DADOS 


LSLOCO DE MA-| 
NIPULAÇÃO 
DE DADCS) 


ET, RACIONAL 


LEITURA 
DO FICHEIRO 
DOS DADOS 
4 MODIFICAR 


Ê 
= 
ERRADA 


ESPAÇO 
LE EstaDDS 


ENTRADA 


nmEMSaGEM 


EXPLICAÇÃO 
me 
TERMINAL 


DE 


DE ERRD 


dê 
fERMIN AL 


FTRAN 


TRANSP ORMAÇ 


ESPAÇO 
DIScÔ TERMINAL 1? E MNE MÔNICA 
| ERRADA 
LEITURA EXPLICAÇÃO ME NG CE qa ; 
oo NO DE ERRO ESCRITA DE 
FICHEIRO TERMINAL wi TERMINAL 
r “O TEQMINAL 


ENTRADA 
DE DADOS 
PELO 
TERMINAL 


Cons TITUIÇÃO 
DE FIC HEIRO 


COM OS RE- 
SULTADOS 


INÍCIO DE NOVA 
TRANSFORMAÇÃO 


CONSTITUIÇÃO 
Dum FICHEIRO 
NO ISCO | 
COM DADO 


PRETENDE 
SUBADAR OS RE- 


ET [ET [MOD CÓPIA SUPERVISOR NEMO Nica 
RáC. GFACIÁDADOS ERRADA : 
| COSTITUIÇÃO 
, DE FICHEIRO 
(4) (13) () | ae | A RETURN cinta: sic] 5 Md da com ps na- 
TRUT 
bu Teu DE - a, DE ERRO DOS E vbME 
Dado PELO 
TE EHIM AL Nº TERMINAL | DOR 


FICHEIROS 


FIG. 4 
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LEITURA 
DO FICHEIRO 
DOS DADOS 
A MODIFICAR| 


DE DADOS rat 
4) RETURN 
DE 


entre as diferentes representações, servindo de apoio a 


ESPAÇO 


ET. FACTORIZADA li 
ESTADOS 


LEITURA 
DO FICHEIRO 


DOS DADOS 
& MODIFICAR 


ENTRADA 


ENTRADA 


DE 


| TRANSFORMAÇÃO 


ESPAÇO 


= FTFAC 


TRANSFORMAÇÃO 


COPIA 


MODIFICAÇÃO SUPERVISOR. 


“puteyr” 


FIÉHEIROS 


PRE TENDE 
CUARDAR OS RES 
SUuLtaDos sô 
DISCO 


CONSTITUIÇÃO 


DE FICHEIRO 
COM OS RE Suyi- 
TãaD2ã E NOME 
pao PELO 
JL IZADÕER 

FIG. 6 


FERIFERES 
DE ExTaa DA 


TERM NE 


En E Ana 


ur Danos 
PLLO TERA 
Padi 


TE Rm dá, 


cont Pit uaDÃo 
Bum FICHEIRO 
PISCÊ 
COM 05 DADO 


to 


Es+ “DD. COPIA Loira dr 
: GEDOS , 
= HETURN pEPLIA EÇçÃãO 
; EM TES DE E 
] a 
FECHE OA ) “Por IisaL 


ESCMITA DE 
RESULTADOS 


so Cla bna) 


1.4. TRANSFORMAÇÃO DE REPRESENTAÇÕES 
1.4.1. Conceitos Básicos 
1.4.1.1. Espaço de Estados 
[x]=[A].[x]+[B].u 
y =[C].[x]+ D.u 
em que 


[A], [B], [C|—-matrizes de coeficientes reais e 
dimensões nxn, nxl e Ixn, respec- 


tivamente. 
D — coeficiente real 
[x] — vector coluna das variáveis de 
estado 
u — entrada 
— saida 


sendo n o número de variáveis de estado, tem-se, 
como limitação de programação 


1&n<€10 
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a = HOM EA 
Disco HER IN AL » 


Eituaa ENTRA O do DT 
/ E paDos Prá 


 d CHEIRÕ | "CREA 
| 


constituição 


E ua FICHE] 


iBanSEDRAÇãO 


ET 
| - EE 


TE RÉCIO 


ml 


FRÉESFE BisAcÃO | 


E de Fui 
CIO Dáiati| 


Vi 
E E midi 


1.4.1.2. Função de Transferência na forma 


racional 
[E a NR o 
Clero a tactl 
Det 2 sererses + ay 
em que 
= =P ERA 
coeficientes reais 
o PONTO “Di 


oc — atraso no tempo. (só para G(s)). É opcional. 


Tem-se n mio 


1.4.1.3. Função de Transferência factorizada 


B—Z)(9— Za) s—z 
F. T=K ( né) (2 a 1) ( 1) e-ucs 

(P—Dip) « (8 — Papos) 000.0 [> Pi) 

em que 

Zoro e, Ty — zeros do sistema 

Pini sresia |, Pp; — polos do sistema 

K — factor de ganho 

2,0 — atraso no tempo. (Só para G(s)). É 

opcional, 


Tem-se n,, n, € 10 
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1.4.2. Espaços de Estados —* Função 
de Transferência racional 


Esta transformação é realizada pela sub-rotina FTRAN, 
de acordo com o fluxograma seguinte: 


1 S Eiji 
A - — 2 EM) 


=1 


IEA SÃO DE | 


F 
| BI-A ] 


E ti JISELÃS) Eos, 


FIG. 9 


Relativamente à inversão da matriz |sI-A|, tem-se: 


| Adi [SA] | 
AO RETURN 

| sI—A | 
Adj [sI-A]= [Ra] +[Roa)r2+...+[R]s+ [Ro] FIG. 11 
|s-A| =s"+0c,, st+...... +00, s+Oy 


1.4.3. Função de Transferência racional — Fun- 


Os coeficientes oc, e as matrizes [R,| são deter- cão de transferência factorizada 


minados na sub-rotina FADIN, utilizando o algoritmo de 
Faddeev [3]. Esta transformação é realizada pela sub-rotina FTFAC, 


FTFAC 


PRBM 
RAITES DUM POLINOMIO 


[NHUMERADOR = FEROS |) 


[o APROX? | 


TRATAMENTO DE NÚMEROS | 


COMPLEXOS | PEQUENOS 
(ZEROS) 


— PRBM 
RATES DUM POLINDMIO 


| DENDMINADOR = POLOS ] 


RIR o t))+m 
F-E n-E b 


nea+a 
v=3,n | APROX? 


TRATAMENTO DE MUMEROS 


| COMPLEXOS PEQUENOS 
(POLOS ] 


“RETUR 
FIG. 12 


A sub-rotina PRBM, pertencente à biblioteca do 
DEC10 do LNEC, (Scientific Subroutine Package-lBM), 
destina-se a determinar as raizes de qualquer polinômio 
de coeficientes reais. 


TÉCNICA 454 101 


APROX? | | 


biz Mig aja +) b| vin 


RETURN 


FIG. 13 


1.4.4. Função de Transferência factorizada — Fun- 
ção de Transferência racional 


Esta transformação é realizada pela sub-rotina FTRAC. 


FIG. 15 


SEPARE 45 BAÍZES REAIS 
Das COMPLEXAS COmIu 
páDis (POLOS) 


FRODUIO DE POLINÔMIOS 


(DENOMINADOS DE Ef) 


APROXI 
TRATAMENTO DE NÚMEROS 
| REAIS PEQUENOS (DEND- 
HIMSDOR Fr) 


SEPARA &S RÁAIZES REM 
DAS COMPLERAS CONM= | 
FADAS ZEROS] 


PRODUTO DE Plk 
Lima Chá DS 


eifaja pin) (00 4a) 


PRODUTO DE POLinQaiDs 


(MUMPRADOR DE FT) 


e > 5 


TRATAMENTO DE MUMES 
ROS REAIS PEQUENAS 
fE DEFICIENTES mM RADOR EI) 


RETURN. 


FIG. 14 FIG. 16 
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A sub-rotina PMPY, que realiza o produto de polinó- 
mios de coeficientes reais, pertence à biblioteca do 
DEC-10 do LNEC (SSP-IBM). 


IE 


RETURN | 


FIG. 17 


1.4.5. Função de Transferência racional — Espaço 
de Estados 


Esta transformação só é possivel se m É n. 
Sendo dada a função de transferência 


Ea RE a saercasses 
O RR gg 
Stars SI+......+ão 


tem-se, para m<n: 


a 
| = Qua 
[B]J=|b |I[c]-[ o [1] 
b; 
” D E» 0 
Ud) 
O 
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[A] é a matriz companheira que, escrita nesta forma, 
permite obter a matriz [B] por simples substituição dos 
coeficientes do polinómio do mumerador da função de 


transferência. 
Se m=n, tem-se 


4 Di E essas ses PD 
e A ii TO : 
s" ta, sTi+...+a, 


sº ta, S"i+...tas 
gendo 
m=n—1 


by = B, —B, ER 


Neste caso, tem-se: 


— [A], [B] e [C] são as matrizes do Espaço de 
Estados determinados para a função de transferência, 


0 
“em que m<n. 


Esta transformação é realizada pela sub-rotina RMIN. 


| 
| 
mi: | 


FIG. 18 
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8/3] =) a) - ttmia 


RETURN 


FIG. 19 


1.5. ORGANIZAÇÃO DA INFORMAÇÃO 


Toda a informação, quer seja dada pelo utilizador ou 
gerada pelo programa, é organizada em ficheiros no 
disco. Esta forma de tratar os dados equivale a encarar 
o disco como uma memória ilimitada, acessivel a todos 
os blocos. 


1.5.1. Constituição de ficheiros 


Qualquer ficheiro no disco é individualizado por uma 
variável de precisão dupla, com a seguinte forma: 


“NOME. EXT” 


em que 


“NOME” — Variável alfanumérica até 5 caracteres (nome 
do ficheiro. 

— Conjunto de 3 letras impostas pelo programa 
indicando o conteúdo ou a finalidade do 
ficheiro. 


*“EXT" 


Os ficheiros, que permitem transmitir informação 
dentro da estrutura global, são escritos em representação 
binária. 

Nestes ficheiros “EXT” indica a que tipo de sistema 
se referem, tendo-se previsto as seguintes possibilidades: 
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ABG — sistema controlado (G) 


ABC — malha de compensação (C) | "2 forma de espaço 


ABH — malha de retroacção (H) de estados 
GSR— G(s) 

CSR— C(s) + na forma racional 

HSA — H(s) 

GSF — G(s) 

CSF — C(s) | factorizados 

HSF — H(s) 


Quanto à variável “NOME”, esta pode ser imposta 
pelo programa ou escolhida pelo utilizador. 

Quando imposta pelo programa, “NOME” depende 
do conteúdo do ficheiro, tendo-se: 


G — sistema G(s) 

C — malha de compensação C(s) 

H — malha de retroacção H(s) 
ABCD — Espaço de Estados para G, C ou H. 


Estes ficheiros são gerados automaticamente pelo 
programa sempre que surgem novos dados, 

Quando escolhido pelo utilizador, “NOME” representa 
qualquer variável alfanumérica com o máximo de 5 carac- 
teres. Estes ficheiros são gerados, a pedido do utilizador, 
quando se inicia uma nova transformação. 


[im | 


ESPAÇÕ 
DE EXtagos 


“a z | PLA Cipa E, 
ERRADA 


ME ta SA DÊ Ma 
DE ERRO HO 
freBaimaL 


DEFENIÇÃO DO| | DEFENIÇÃO DO 


| DLFERIÇÃãO DE 


NOME E EMT NOME E EM MDME E EmI* 


TIPO 
DE 
FICHEIR 


NT (ND INPOSTO PÊLO ranSnaMa! 


! 


[NOME DADO 
PELO UTILIZADORES 


COMStiTuIção | 
Da DESENaÇÃãO 
| DO FICHEIRO 
| CnDME-ERT 


INDICAÇÃO, 
PELO utiits 
ZADOR, Dt 

“NOM E” 


Comi TITuaçÃão Dus 
| DESIGNAÇÃO DO 
| FICHEIRO 
MOMESE 


LEI URA DO 


FICHEIRO 
Daio 


FIG. 20 
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IMPRESSORA 


TERMINAL 


ESCRITA DOS 
DADOS MO 


TERMINAL 


CONSTITUIÇÃO 
DO FICHEIRO 
PARA ESCREVER 
Ma IMPRESSORA 


RETURN 


1.5.2. «Output= de ficheiros 


FIG. 21 


O «output» de ficheiros, no terminal ou na impres- 
sora, é realizado pela sub-rotina COPIA (ver figs. 20 e 21). 

Quando se pretende o «output» na impressora, é 
constituido um ficheiro de extensão “LPT” e “NOME” 
dado pelo utilizador. Estes ficheiros são automaticamente 
transcritos logo que o utilizador termine a utilização do 
computador, sendo em seguida destruídos. 


1.6. RELAÇÃO COM O UTILIZADOR 


Esta estrutura, programada em modo conversacional, 
permite a interrupção da sequência de execução em 
diversos pontos pré-definidos do programa, A comunica- 


ção utilizador-programa é realizada por um conjunto de 
mnemónicas, reconhecidas por todos os blocos, 


— Descrição dos sistemas: 


ABG — G 
ABC — C na forma de espaço de estados 
ABH—H 
GSR — G(s) 
CSR— C(s) | na forma racional 
HSR — H(s) 
GSF — G(s) 
CSF — C(s) factorizada 
HSF — H(s) 
—. Periféricos 


TTY — Terminal 


DSK — Disco 
LPT — Impressora 
— Ficheiros 


INT — Ficheiros com nome gerado pelo programa 
EXT — Ficheiros com nome dado pelo utilizador 
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—— Ordens de comando: 


COP — «output» de ficheiros 
MOD — Modificação de dados 
INI — Recomeço deste bloco 
SUP — Regresso ao supervisor 


— Outras hipóteses: 


+ — Pedidos de informação 
S — Sim 
N — Não 


Foi prevista a possibilidade de erro por parte do 
utilizador. Nesse sentido, sempre que a mnemónica dada 
não for reconhecida pelo programa, é impressa uma 
mensagem de erro e repetida a pergunta. 

São também testados alguns dos dados numéricos 
indicados pelo utilizador. Está neste caso, por exemplo, 
o teste efectuado aos pólos e zeros de um sistema, com 
vista a garantir que todos os complexos sejam conjugados. 


1.7. EXEMPLOS 


Nota: Nos exemplos que apresentamos, todos os 
caracteres impressos pelo utilizador encontram- 
-se sublinhados. 


1.7.1. Entrada de Dados 


ELOCO DE MANIPULACÃO DE DADOS 


FORMA IHICIAL DO ETSTENA dos 


FERIFERICO DE ENTRADAS IT 


SISTEMA NA FORMA DE FUNCAO DE TRANSTERENCIA RAGTONAL 


bi) Mi-1 
EIN) St BIN-L) 5 tea + * É Atl? D+ TO) 

ET == qe a a pi a a —. EXPi-ALTA E) 
H M=1 
GE + BImM-IDS Eoscas» dk ANIS D+ ALO) 

EM QUE: 

BidIrÃtd) COEFICIENTES REAIS 

ALFA ATRAZO MO TEMPO (50 PARA G(5/).E DPLIONAL. 


GRAU DO NUMERADÓRIO 


IMPRIMA DO COEFICIENTES PO FOLIMONIU EM 5 DO NLUSERADOR» 
FELA SEGUINTE DRUCHMSBUNSeBiNo=Lorec o rRLLocEttos 


a 
=] 


E 


BRAU DO DENOMINADOR SD. 

IMPRIMA OS COEFICIENTES DO FOLINOMIO EM S DO DENOMINADOR » 
PELA SEGUINTE ORDEM: ACN-IDARN-D) reco rACIIFAÇOD 

2d 49 


O SISTEMA APRESENTA ALGUM ATRÁZO NO TEMPOS 


INDIQUE A REFERENCIA DD SISTEMA 
t ATE S CARACTERES JSTESTE 


TERMINOU A ENTRADA DE DADOS 


FIG. 2> 
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Neste exemplo introduziu-se a função de transfe- 
rência. 
G(s) 49 
8) = ————— 
s?+2.85+49 


utilizando como periférico de entrada o terminal. 


1.7.2. Modificação de Dados 


FORMA VINAL DO SISTEMMANOT 


PRETENIH MOUL ToáR od NUMERADOR fãs 
GRau DO Mui RAD: 
IMPRIMA US COEFICIENTES DO POLTIMONIO EM 5 EO NUMERADOR: 


PELA SEGUINTE ORDEM! BON! BAU L dra cor BlLrRtOdo 
PRETENDE MODIFICAR O DENOMIMADORTH 


PRETENDE MODIFICAR O ATRÁZO NO TEMPOTE 


TLERMINOL à TIA Lia DC Dmlrigã 


FORMA FIN PO EGISTEMA uol 


FIG. 23 


Neste exemplo modificou-se a função de transfe- 
rência dada em 1.71. de forma a ter-se 


s—49 


G(s) =——— 
sº+2.8s+49 


1.7.3 F.T. factorizada — Espaço de Estados 
Dada a função de transferência factorizada 


(s+1-=;2) (s+1+]2) 


O] = a 
(s+1—)3) (s+1+]3) (5+5) 


pediu-se, neste exemplo, a determinação do espaço de 
estados 
FORMA FIMAL DO GISTLARSALDO 


DETERMINAÇÃO DA PUNCAO DE TRANSTERÊNCIA Ná FORMA RACIONAL 
DETERMINADAO DE ADC 


FESTI 


MATRIZ A! 


L 13 6,.0000L+00 
lo Ty 1, 0000L 100 
aa Ds DOODEADO 


MATRIZ E: 


|. 14 5,0000101 
L di 2, 0000L +01 


L 3: 1,0000L+401 


MATRIZ Ci: 


L li o Qa00BGE1DO 


106 


0 DOQDE FOú 
O. DOODE | DO 


À 0000L 100 


T. DGODDE GO 


ra DE A] 


= e DLUDE OL 


+ OUCDEA DU 


Io DOOU 


FIG. 24 


1.7.4. «output» de ficheiros 
É AGORATCOE 


COPIA DE FICHEIROS NO TERMINAL DU NÃ IMPRESSORA 


FORMA DO ad 5 TEMA sl 

TIPO DE FICHEIRO: INT 

PERIFERICO DE SAIDASTTY 

FUNCAD DE TRANSFERENCIA G(5) RACIONAL: 


REFERENCIA:! TELETE 


NUMERADOR 


COEF( 0)=-4, POOOCHÓI 


COEF( 1j= 1,0000E+00 
DENOMINADOR : 


COEF( Oj= 4, 7P000E+01 
COEFU Lj= 2. BLO0E+OO 
COEF( 2)= 1,0000E+400 
PRETENDE CRIAR FICHEIRO DE EXTENSAD LPTYN 


FIM DE COFIAs, 


BLOCO DE MANIFULACÃO DE DADOS. 


FIG. 25 


Este exemplo mostra como proceder no caso de se 


pretender escrever no terminal a função de transferência 
dada em 1.7.2, que se encontrava memorizada num 
ficheiro gerado internamente pelo programa. 


1.7.5. Recomeço do Bloco de Manipulação 
de Dados 


E ABORATINI 
PRETENDE GUARDAR A DESCRICAO INICIAL DO SISTEMA NO DISCO?S 
NOME DO FICHEIRO(ATE 5 CARACTERES) TESTE 

E À DESCRICAO FINAL?N 


ELOCO DE MANIPULAÇÃO DE DADOS 


FORMA INICIAL DO SISTEMÁLOSE 
FIG. 26 
Pedido o recomeço do Bloco de Manipulação de 


Dados, constituiu-se um ficheiro de “NOME” = “TESTE” 
com a descrição inicial do sistema, 


1.7.6. Mensagens explicativas 


Nos dois exemplos que a seguir se apresentam, 


mostra-se a resposta do programa ao pedido de expli- 


cação do utilizador (impressão do caractere?),. 
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E AGDRATT 


1.7.6.1. çE qsção. axa ditad 


às SEGUINTES DPCÕES SAO POGSIVETS: 


THI 
SUP 
ROL 
EDEF 


INICIAR à TRANGFORMACÃO DUM NOVO GISTEMA 

REGRESSAR AD GUPERVISOR UND : 

MODIFICAR UM CONJUNTO DE DADOStE RECOMECAR ESTE DLUCO) 
REALIZAR O OUTPUT DE FICHEIROS = 


JT 


E AGORATIS 


FIG. 27 


1.7.6.2. 


FORMACFENAL DO SISTEMAS? 


FODE DFTAR POR UMA DAS SEGUINTES FORMAS 


LSF 
AEE 


FUNCAO DE TRANSFERENCIA D(5) FPACTORIZADA 
SISTEMA N& FORMA DE ESPACO DE ESTADOS 


DU» AINDAs INTERROMPER COM UM DOS GEGUINTES COMANDOS: 


Eur 
MOL 
COF 
INI 


REGRESSAR AD SUPERVIGOR UNIV. 

MODIFICAR UM CONJUNTO DE DADOSt(E RECONECAR ESTE PLOCC) 
REALIZAR OD DUTPUT DE FICHEIROS 

INICIAR À TRANGFORHACÃO DUM NOVO CISTEHA 


FORMA FINAL DO SISTEMA:INT 


FIG. 28 


1.7.7. Detecção de erros 


Ltda, 


FERIFERICO DE ENTRADAS DAS 


à MNEMONICA ESTA CRRADA. 
IMPRIMA T PARA MAIS INFORMAÇÃO. 
TENTE NOVAMENTE 


PERIFERICO DE ENTRADA: DSI 


FIG. 29 | 


O erro cometido pelo utilizador ao imprimir a mnemó- 


nica do 


periférico (disco) foi detectado pelo sistema. 


TÉCNICA 454 


1,7,7.2. 
FORMA INICIAL DO SISTEMA ZOSE 
| 


PERIFERICO DE ENTRADASITY 


—— SESTEMA-NA FORMA DE-FUNCAD DE TRANSFERENCIA FACTORIZADA 


| 
FT = FG === cnc cocummenes EXPI-ALFA 5) 
(CS-P1) CS-PD) a o q + €5-PM) 


EM QUE: 

Za ZEROS DE FT 

PJ POLOS DE FT 

FG FACTOR DE GâáNHO 

ALFA ATRAZO NO TEMPOÇSO PARA G(5)), E DFCIONAL 


NUMERO DE ZEROS!1S 
O SISTEMA NAD ACEITA POLINOMIOS DE GRAU TÃO ELEVADO 
NUMERO DE ZEROS!1 


FIG. 30 


A indicação de um número superior a 10 para o 
n.º de zeros de um sistema, conduz à não aceitação 
da entrada. | 


ITS. 


Este exemplo refere-se à definição do numérador da 
função de transferência factorizada. À indicação de zeros 


“que não sejam complexos conjugados, conduz à não 


aceitação da entrada, 


NUMERO DE ZEROS!Z 
IMPRIMA OS ZEROS DA FUNCAO DE TRANSFERENCIA. 


ZERO 1ii 2 


Dad 


LERU 241 5 


— 


OS ZEROS TEM DE SER COMPLEXOS CONJUGADOS. 
TENTÊ NOVAMENTES O = 


IMPRIMA OS ZEROS DA FUNCAO DE TRANSFERENCIA, 


FIG. 31 
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ANEXO: Demonstração do algoritmo de Faddeev 


Designando por [sl-A]* a matriz obtida por aplicação do algoritmo de Faddeev, 


[Ra] s"1+ [Rg] s"2+...+[R;]s+ [Ro] 


[sI-A]|*= 
| sI—A | 


em que 

(1) |st-Aj=s"+ cy si +... + oc, s + cc, (equação caracteristica da matriz [A|) 
(2.4) [Rgal=[1] 

(28) [Rix]=[A] [Ro-pa] + Cota [1] 


= [AJ + oq [AJ + ops [AJeS +... + Ca [A] + Coy [1], Vk=2n 


Tem-se: 
[Ras] st + [Roo] 62 +... + [R] + [Ro] 
|  |si-A] | 


ei tforA)e (Roo) 624 [Rg] stê 4,4 [Rijo [RD = 
|sI—A | 


[sl-A]- [sl-A]*= [sI-A] - 


E 


Ga (Uaa] st + [Bal 6º 4 + [RA] 0º 4 [A] 6 —[AT Ra] te 


— [A] [Ros] s"*—..— [A] [R]s—[A][R]) 


Atendendo às relações (2.4) e (2.B), tem-se: 


[sI—A] - [si—A |* = RR (1) sr + oca [st +... + os [1]s?+ oc, [1] s—[A] [R5]) 


Particularizando a relação (2.B) para k = n, vem: 


(1) sr + cs [si +... + ce [iJst+ o [1] s— 


(3) [sI-A] s [sl-A]*= | 


— ([A]"+ oc [Aju +... + oco [A]2+ oc, [AT)) 


Atendendo ao teorema de Hamilton-Cayley — qualquer matriz quadrada satisfaz à sua equação caracteristica 
— podemos escrever; 


[AJ + ca [AJ +... + ce, [A] + o [1] = [0] 


(4) [LAJE + Os [AJ +... + 06, [A] =— 00, [1] 
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Substituindo (4) em (3), tem-se: 


[sl-A] - [sl-A]|*= cla Ms" + [st + oco [Ii]? +... + oc [1] s—(— oc, [1])) 


|sI—A | 


” rr (LI)? + oera [1] st + 06gog [1] ont +... + ces [Jo + oco [1]) 


meo 


(84 + Cos Lp 6, d+) De ee e) 
1 | l S 


| si—A | 


a a 


O cc. O (Mt. si+.. +, s+ cc) 
Atendendo à relação (1), vem: 
ta Pta to BHO Qercccs 0 E 


[si—A] - [s-AJ*= 
sh + Co gr RR o | 5 + e 


4 ; 0 
Dea dreitia O M+r-X sit... +0C, s+ Oo, 


sS1+tC o sPl+..+O, 8+ € 
Atendendo a que, por definição 
[sI-A] - [sl-=A]-* = [1] 
então 


[si—A [*= [sIi—A]— 
c.d.d, 


[1] C. Limpen, E. Holmyard, A. R. Hall, T. S. Wil- [2] D'Azzo, C. H. Houpis — «Feedback Control System 
liam — «A History of Technology», Oxford Univer- Analysis & Synthesis», Int. Student ed., 1966. 
sity Press, 1956. [3] H. H. Rosenbrock — «Computer-Aided Control Sys- 
tem Design», Academic Press, 1974. 
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ANO LIV: 


VOLUME XLI 


Modelisation et stabilité de I'association onduleur 
de courant-machine asynchrone * 


RESUMO 


A associação ondulador de corrente-máquina assin- 
crona constitui um sistema para o qual as equações dos 
circuitos de contactos deslizantes são de aplicação 
imediata, Era, 

A complexidade do estudo da estabilidade do sistema 
leva-nos a introduzir hipóteses de simplificação que o 
tornam viável analiticamente. Desta maneira é possivel 
considerar dois tipos de instabilidades: eléctricas e me- 
cânicas. 

A modificação dos parâmetros para assegurar a 
estabilidade de associação ondulador de corrente-máquina 
assincrona não pode ser vista como uma solução, porque 
assim se é conduzido a sistemas de baixo rendimento. 

Com a introdução de cadeias de retroacção obtém-se 
sistemas estáveis sem modificar as boas performances 
da máquina assincrona. Os controlos da corrente con- 
tinua de alimentação e da frequência de escorregamento, 
com efeitos complementares, permitem eliminar os tipos 
de instabilidade considerados. 


1 — INTRODUCTION 


L'introduction des élêments semiconducteurs de puis- 
sance dans le domaine de la conversion électromécanique 
denergie a rendu possible régulation de vitesse de la 
machine asynchrone par variation statique de la fréquence 
d'alimentation. Les avantages de ce systéme face aux 
machines à courant continu résident dans la fiabilité et 
prix des machines d'induction à cage. Avec I'évolution 
du prix des dispositifs semiconducteurs de puissance, 
on prévoit que les machines à courant altematif rempla- 
ceront les ensembles classiques dans les entrainement à 
vitesse variable, 

Un des systêmes le plus souvent étudié actuelle- 
ment est celui de |'association de I'onduleur de courant 
à la machine asynchrone, Le systême résultant, se 
présente comme une machine à commutation électronique 
à laquelle il est possible d'appliquer la théorie des 


J. SANTANA 
Centro de Automática da Universidade Técnica de Lisboa 
Departamento de Comportamento de Sistemas (1.5.T.) 


SUMMARY 


This paper presents a dynamic model for the induc- 
tion motor driven by current source inverter, which is 
developped as an application of the commutation circuits 
theory. 

An analytical study of the system stability is possible 
with the introduction of simplifying hypothesis in the 
mathematical model. We will consider two types of insta- 
bilities: electrical and mechanical. 

k is shown that the actions on the parameters to 
obtain the stability induce to systems with poor efficiency. 
The solution is the introduction of the feed-back loops. 

The controls of the direct current and of the slip 
frequency with complementary effects leave out the two 
types of instabilities: the first one is equivalent to the 
increase of the stator resistance and so it is able to 
eliminate the electrical instability, the second one elimi- 
nates the steady instability bringing an electromechanical 
decrescent characteristics in the system. 


circuits à contacts glissants [1] [2] Le modéle mathéma- 
tique obtenu représente assez bien le comportement du 
systême, qui s'écarte beaucoup du fonctionnement de la 
machine asynchrone classique., 

Comme les équations decrivant le systême sont 
non-linéaires, létude de la stabilité electromécanique 
se fait par linêarisation autour du point de fonctionnement, 
pour des petites perturbations [3], [4], [5]. Vu sa com- 
plexitê I'étude analytique de la stabilitê n'est pas possible, 
que si on introduit des hipothêses simplificatrices, clas- 
siques dans ['étude des machines électriques. De cette 
façon on vérifie qu'il est possible de considerer deux 
types d'instabilitês: électrique et mécanique, et quon 
peut analyser les mécanismes qui les determinent. Des 
expressions simples, donnent des indications sur quels 
paramétres on doit agir pour obtenir des systêmes 


* Travail fait en collaboration avec |Unité FORT de |'Université de Louvaln. 
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stables, Les conclusions ainsi obtênues sont en accord 
avec celles presentées pour I'étude de la stabilité globale, 
c'est-â-dire, sans hypothêses simplificatrices, faite par 
voie numérique dans la reférence [5]. 


On vérifie que l'adaptation des parâmetres du sys- 
tême pour assurer qu'il soit stable dans toute la zone 
de fonctionnement moteur, conduit à de mauvais ensem- 
bles du point de vue du redement. C'est donc avec 
introduction de boucles de rétroaction qu'on obtient la 
stabilisation dans de bonnes conditions. À cause de leur 
importance, les effets des controles du courant continu 
d'alimentation et de la fréquence de glissement sont 
analysés. On observe qu'ils sont complémentaires du 
point de vue de la stabilité: tandis que le premier est 
équivalent à une augmentation de la resistance stato- 
rique et permet d'éliminer les instabilités électriques, 
limposition d'une fréquence de glissement constante 
donne à la caractéristique électromécanique un caractêre 
de fonction monotone decroissante et elle lêve les insta- 
bilités mécaniques. 


2 — MODELISATION 


La description du fonctionnement de l'association de 
l'onduleur de courant à la machine asynchrone peut se 
résumer comme suite: des signaux decalés entre eux 
de 60º et élaborés par le circuit de commande sont 
envoyés au gachettes des thyristors, de façon que chacun 
conduise pendant 120º dans chaque période, Le passage 
d'un thyristor aux thyristor adjacent se fait par le déclen- 
chement de ce dernier puisque la configuration qu'on 
obtient avec cet amorçage fait appliquer au thyristor qui 
était en conduction une tension inverse fournie par les 
condensateurs placés dans le montage. Dans la figure 1, 
les thyristors A et F sont en conduction et les conden- 
sateurs sont chargés comme indiqué. Ainsi, quand on 
amorce B, la tension des condensateurs s'applique en 
inverse au thyristor À et provoque son blocage, 


FIG. 1 


La self extérieure qui lisse le courant continu et la 
présence des condensateurs qui assurent la commutation 
dune phase à lI'autre de façon presque instantanée, 
donnent au courants de phase de la machine la forme 
Suivante: 
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FIG. 2 


Un équivalent mécanique du systêéme est donné par 
la représentation de la FIG. 3, ou nous avons les trois 
phases du stator lies à un collecteur mécanique sur 
lequel tournent des balais qui font un angle 8, et I'enrou- 
lement (x; 8) qui tourne en faisant un angle O, par 
rapport à une position de réference. 


A Ja, 


FIG. 3 


Vu des terminaux A, B lI'association onduleur de 
courant-machine asynchrone se présente comme une 
machine à commutation à laquelle il est possible d'appli- 
quer la théorie des circuits à contacts glissants, || est 
interessant de noter qu'en 1968 dans la réference [6] on 
fait déjã I'étude théorique du regime permanent pour un 
systême semblable à celui de la FIG. 3. 

La condition rigoureuse d'applicabilité de la théorie 
des circuit du type à contacts glissants est d'avoir un 
nombre infini de branchements à commuter de façon 
discontinue. En vérifiant que le transfert d'une phase à 
l'autre se fait de façon presque instantanée et que la 
correction non-lineique [7] due au nombre fini de bran- 
chements (six) n'apporte guêre des modifications, on peut 
dire que se systême est une illustration directe de 
application des équations des circuits à contacts glis- 
sants. 

Les équations de l|'association onduleur de courant- 
-machine asynchrone sont: 


dy, 0-0, q 
U, nc a la + + Via (1 ) 
de à 90, 
0. ()= d 0, (2) 
dt 
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dy, 


“oi PU (3) 
; db 
A cj4— 4 
O a gt = (4) 
dO 
T—T,=] =D; (5) 
dt 


(1) est l'éêquation du circuit du stator qui est du 
type de circuit à contacts glissants, L'équation (2) donne 
la position des balais, definie par la logique de commande 
de Ionduleur. (3) et (4) sont les équations du rotor. 
(5) traduit la dynamique de la machine avec la charge. 

Dans ce modéle mathématique les variables i, et ig 
(courants du rotor) ne sont pas constantes en régime 
permanent, Aussi on change la representation des gran- 
deurs du rotor ou l|'enroulement biphasé. (a, B) sera 
remplacé par un autre placé dans un referentiel qui 
tourne avec les balais du stator. 


De cette façon on obtient des systémes d'équations 
dans lesquels les variables sont constantes en regime 
permanent, ceci étant important, en particulier dans I'étude 
de la stabilitéê quand on fait la linéarisation autour d'un 
point de fonctionnement, 


Les équations du systêmes, aprés la transformation 
(x ,B)— (d,q) des grandeurs du rotor, sont: 


U, Y=/R+Lp Mp Mm 4/1, (O 
O M p A+L p d. L, lar 
O | — br M — O L, A +L p lar 
(U) = (2) (1) 
q d 
o.()=— 0 2 
s (b) o (2) 


"TR 
Te-—Tp=] + D 0. avec Te=M R lar (7) 
dt 


oú 0. =0.—0 


Un autre modéle ou les grandeurs du rotor sont 
constantes est obtenu en assimilant la force magnéto- 
-motrice du rotor, exprimée dans les axes dq à une 
f.m.m, ni, avec une position Y par rapport à la ligne 
des balais du stator. 


Les équations obtenues avec cette transformation 
sont: 


di di El 

U,=Rilatli-—+McosY—+(0,-Y)Msin Yi, (8) 
dt dt 

di, di: à 

O=Mcosy-— +R.L+L,—— —0 Msinri, (9) 

dt dt 


dl * a 
O=Msiny — +80, Mcos Yi, +(0,—Y)Lçi (10) 
dt 
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L'équation de I'équilibre mécanique est bien sur la 
même dans tous les modéles. Dans ce cas, le couple 
électromagnétique est donné par |'expression: 


Te=Mii siny (11) 


|| est interessant de Voir la similitude de ces équa- 
tions et de celles qui décrivent la génératrice à redres- 
seur [7]. fi | 


(8) est |'équation de Vinducteur, que peut être 
presentêe sous la forme suivante: 


la 


dt 


U,—8.MsinYyi=A,i+L, + 


dic | 
+Mcosy-—— — YMsinY% ii, (12) 
dt 


C'est-â-dire, que la tension appliquée à l'inducteur de 
la génératrice équivalente à notre systême, premier 
membre de lI'équation, est le résultat d'une boucle de 
retroaction du à la force contre-électromotrice 


— B. Mi sinYy sur la tension appliquée 
(9) est I'équation de Iinduit de la génératrice équi- 
valente, celle-ci tournant avec une vitesse 0, et étant 


mise en court-circuit, puisque le premier membre de 
l'êquations est égal à zéro. 


(10) est I'équation qui nous donne la position des 
balais de la génératrice (appelée équation électromé- 
canique). 


L'equivalence qui vient d'être décrite est presentée 
sous une forme graphique à la FIG. 4, 


FIG. 4 


A de petites différences prês le deuxiême modéle 
mathématique developpé est identique à celui présenté 
aux réferences [3] et [4]. Cependant comme il est déduit 
d'une struture mathématique d'application géênérale aux 
machines à commutation, son établissement ne depend 
pas de la forme de la loi de variation des inductançes 
mutuelles stator-rotor, contrairement aux modéles deve- 
loppés à partir de la transformation de Park ou il faut 
admettre que cette loi est du type cosinusoidale, Aussi 
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ce modéle pourra être genéralisé à d'autres machines, 
par exemple de grand nombre de phases oú |'hypothêse 
de variation cosinusoidale de la mutuelle n'est pas 
aoceptablea. 


3 — STABILITÊ 


Les problêmes de stabilité de |'association onduleur 
de courant-machine asynchrone en boucle ouverte cons- 
tituent une des principales limitations à son utilisation. 

Dans la mesure ou les équations du systême sont 
non-linéaires, l'étude de la stabilité electromécanique se 
fait par linéarisation, autour du point de fonctionnement, 
pour de petites perturbations. Une étude globale, qui 
est faite par une voie ordinateur dans la reférence [5], 
se revéle impossible de façon analytique, Cependant, si 
on fait des hypothêses simplificatrices, on arrive à des 
structures plus faciles à traiter. 


3.1 — Stabiblité electrique 


Si on suppose constante la vitesse du rotor, I'in- 
fluence de la partie mécanique sur la stabilité du systême 
disparait, ce qui nous reduit l'analyse au «Calcul de la 
Stabilité Electrique», De cette façon Il'étude se réduit à 


l!'analyse de I'éguation (6) avec 8. constant, Elle peut se 
présenter sous la forme matricislle suivante: 


UN—/QMiu= 
O | O) 
O o 
(13) 
= / RtLp Mp (6) la 
M p R,+ Ls Pp Fte Or Ei ldr 
— 0, M =, L, + bs 8) lgr 
Ou 
(U', O 0)! =(Z;) (is lar lar)! (14) 
ou 


Us, = U,—, M gr 


Le schéma-bloc correspondant est: 


Une analogie entre ce systême avec un alternateur 
en court-circuit, à tension d'excitation controlée, est éta- 
blie à partir de I'équation (13). En effet, cette équation 
peut être vue comme I'équation d'un alternateur, qui 


tourne avec la vitesse 0... et ou les enroulements (d q) 
d'induit sont en court-circuit. La tension d'excitation U', 
est le résultat de la tension continue U, et de la chaine 


de retroaction, qui donne une tension — OM ig C'est-ã- 
-dire, on remplace le circuit statorique, par un circuit 
du type filiforme ou le terme de glissement apparait 
comme un effet de retroaction, 

Avec |l'hypothêse d'une vitesse du rotor constante, 
[alternateur en court-circuit est toujours stable, des lors 
la fonction de transfert igr/U', présente ses pôles dans 
le demi-plan de gauche, Aussi on peut dire que les 
instabilitê electriques ont pour origine l'action de la 
boucle de retroaction. Les FIGS. 6a,b répresentent la 
localisation des pôles sans et avec rétroaction pour 
différentes valeurs du glissement. 

Comme on le voit les pôles de lalternateur en c.c. 
se trouvent dans le demi-plan de gauche FIG. 6a), mais 
quand on «ferme» la chaine de retroaction, on a un 
déplacement des pôles 8 jy vers la droite, FIG 6b), 
ce qui peut conduire le systême à linstabilité, pour un 
certain domaine de valeurs du glissements, 

La stabilitê du systême est dautant mieu garantie 
que les pôles en boucle ouverte solent le plus possible 
à gauche et que l'action du gain en boucle fermée soit 
la plus faible possible. On montre que la somme des 
parties réelles des pôles en boucle ouverte et que le 
gain sont donnês par les expressions suivantes: 


FIG. 6 
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a+2Bb=—-. — —2. * 4k2. (15) 
SA L. Es 
E R. 
Kp=s (0,K). (16) 
ou | 
L',=0L L=0L, 
AZ 
g=1-K=1— M 
La Lo 


Ces expressions constituent des «mesures» de la 
stabilité puisqu'elle donnent des indications sur les para- 
métres sur lesqueles on doit agir pour obtenir des sys- 
têmes stables. 

La FIG. 7 donne la localisation des pôles en boucle 
ferméa pour différentes valeurs de la résistance statorique 
(par simplification sont presentés uniquement ceux qui 
peuvent se déplacer vers le demi-plan de droite), 


FIG, 7 


Une augmentation de la résistance statorique est 
favorable du point de vue de la stabilité, elle fait augmen- 
ter [a +28]. 

La FIG. 8 donne la localisation des póles en boucle 
fermée pour de différentes valeurs de la fréquence stato- 
rique (sont presentés uniquement les pôles B = jYy par 
simplification). 


FIG. 8 


La diminution de la fréquence statorique a une action 
favorable sur la stabilitê du systême. Puisque dans la 
définition du gain Kp, la fréquence apparait associé au 
coefficient de couplage K, la conclusion tirée de la figure 
precédente est aussi valable pour la variation de celui-ci. 
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3.2 — Stabilité mécanique 


L'autre hypothêse qu'on fait, néglige les variations 
de flux dans I'équation (6) ce qui conduit au calcul dit 
«Calcul de la Stabilité Mécanique-, Dans ce cas, I'étude 
de la stabilité du systême se réduit à lI'analyse du 
schême-bloc de la figure suivante. 


FIG. 9 


ou G= 


est la pente de la caractéristique électro- 
o, 
mécanique. 
La condition de stabilité mécanique est donnée par: 


G-D< O (17) 


qui est la condition connue de stabilité quand on analyse 
les regimes transitoires avec les équations du regime 
permanent electrique. 

L éêxpression du couple électromagnetique est donnee 


par: 
0. MR U; 
EA (18) 
cu er E E 
A+ UM L R.+ O. 0. M R, 
R, +04 ll 


Aussi la forme de la caractéristique électroméca- 


nique dépend du fait que la condition (Mº I/L, Rj) < 2 
soit verifiée (FIG. 10b) ou non (FIG. 10 a). 


FIG. 10 


Quand le systême présente la caractéristique électro- 
mécanique de la forme b) et quand on a L,/R,= L/R, 
on peut montrer que le systême est stable electriquement 
dans tout le domaine du fonctionnement en moteur. On 
remarque que les paramétres qui determinant la stabilité 
électrique sont les mêémes que ceux qui déterminent la 
forme de la carateristique électromécanique. 


3.3 — Stabilitê electromécanique 


Comme indiqué la complexité du calcul de la stabi- 
litê électromécanique ne permet pas qu'il soit fait analyti- 
quement. Cependant, sous certaines conditions, ce com- 
portement peut étre prévu à partir des études partielles 
des paragraphes (3.1) et (3.2). 
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La méthode de linéarisation autour du point de 
fonctionnement pour de petites perturbations permet 
d'établir toutes les fonctions de transfert du systéme. 
En particulier, la fonction AT./AT, est trés utilisée pour 
lêtude de la stabilitê du systême. Si on néglige le 
coefficient de frottement (D=0O) elle est representé par 
le schême-bolc suivant: 


ou 
£ .2 RE Ui 
Kem= M (RL) n (A; + Or Er) É (1 9) 
] 
p'-+a./a, p'+asla 
GR e 
p (p'+do/d, pº+da/d, p+dild,) 
avec 
d =0LL, 
d=R,L+RLLA+R LL 
dy =2R,R,L+RL,+0 LoL, o 
d, Es E A; + R, UM à, ME + OS RR ia 
et 


a = [LL MRI RE+OLY] (UA Z) 
as = [0,0,L,M A(1-3S)+0 OL LÍLMR,— 


—RM (R,R;L,+0, L/R, +RÇL)] (U,/3 2, 


a =G.á2Z, 


G est la pente de la caractéristique électromécanique 
du schéma-bloc de la FIG. 9 et 4Z, est donné par 
expression: 


A És HR, pó +8.. 0. Mº R, +04 LE BR. 


Pour toute valeur du glissement, on connait la locali- 
sation des pôles de la fonction dT /AT, en boucle 
ouverte, c'est-á-dire, de la fonction G (p), dont le domi- 
nateur est le produit de p par le polynôme caractéristique 
de la matrice (Z) (équation 6). 

Le pôle à ['origine est du à la prise en considération 
de la charge mécanique du systême, tandis que les 
autres sont ceux étudiés dans le calcul de la stabilité 
electrique au paragraphe 3.1. 

Pour Kem 50, ce qu'on vérifie pour ]J=% ou 
U,=0, I'étude de la stabilité electromécanique coincide 
avec celle réalisée en (3.1,). 
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Pour Kem + 0, en plus de la connaissance des pôles 
en boucle ouverte, on doit connaitre la localisation des 
zéros de la fonction 4T,./àT,. Ceux-ci sont les racines 
de |'êquation 


ajp'+a,p+a,=0 (21) 


Comme le signe de a, est négatif, quand a, > O quel 
que soit le signe de a., nous aurons deux zéros réels, 
dont un est positif, FIG, 12, 


FIG. 12 


De cette façon, pour une valeur donnée du glisse- 
ment, quand U,*/]— = le pôle p, se déplace vers le 
zéro Z., et le systême devient instable. Mais a,>0 
signifie que G > 0 et alors on peut conclure que [insta- 
bilitê mécanique determine linstabilitê électromecanique. 

On vérifie que le signe de a.,, sauf quand on est 
três proche de la vitesse de synchronisme, est positif. 
Dans cette condition et avec a, <D eta, < 0, les racines 
de l'êquation (21) sont conjuguées avec partie réelle 
positive En supposant que le systême est stable électri- 
quement pour une valeur donnee du glissement, quand 
U,*/]—> =, on aura la localisation des pôles presentée 
à la FIG. 13. 


FIG. 13 


On voit que si le systême est stable electriquement, 
alors il y a un domaine de U,*/] dans lequel il I'est 
aussi electromécaniquement. On montre que l'importance 
de ce domaine est tel qu'on peut dire que la stabilité elec- 
trique determine la stabilité électromécanique, quand 
a, <o0. 


41 — Stabilisation du systême 


Pour la stabilisation de l'association onduleur de 
courant-machine asynchrone deux solutions se presentent: 


a) Modification des paramétres 
b) Introduction des boucles de rétroaction 
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L'alternative a) conduit à des systêémes assez mau- 
vais du point de vue de rendement, En effet comme 
indiqué en (3.1), les paramêétres dont les variations con- 
duisent à la stabilité electrique, soient la résistance A, 
et le coefficient de couplage K, doivent prendre des 
valeurs qui, comme on peut le voir dans les FIG 7,e 8, 
font augmenter beaucoup les pertes du systêéme. De cette 
façon, l'introduction des boucles de retroaction est la 
solution qu'on recommande et cest dans cette direction 
que les recherches trouvées dans la littérature se sont 
orientées. 

On va présenter deux contrôles, qui ayant des effets 
complémentaires, lêvent les deux types d'instabilitês trou- 
vêes dans le systême., 


4.2 — Effet du contrôle du courant continu 
sur la stabilité 


On agit sur la commande du redresseur, FIG. 14, de 
façon que le courant continu, |, SOit constant, 


FIG. 14 


En donnant à lalimentation un caractére de source 
de courant cette regulation a une action três favorable 
pour la stabilisation du systême. Aussi on trouve dans 
certains articles [inclusion de ce controle dans la struc- 
ture de [association onduleur de courant-machine asyn- 
chrone [9], [10], [11]. 

Pour montrer influence du controle du courant i, sur 
la stabilité, on considére la vitesse du rotor constante. 
Dans ce cas l'analyse du systême se réduit à I'équa- 
tion (6) plus la relation suivant:; 


U,=K. (1º) (22) 
qui traduit action d'une boucle de régulation du type 


proportionnel. Ainsi, on remplace la variable d'entrês U, 
par le courant de reférence |,*. 


La stabilité électrique de ce systême est connue 
par l'analyse de I'équation suivante: 


KAEN = 
d (23) 
0 
is (R4+KJ+Ls Pp Mp 6.M ha 
M p R-+Lrp Be L lar 
— 8,.M — a L. R4Lp e 
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Comme on le voit, l'effet de ce contrôle se traduit 
par une augmentation apparente de la résistance du 
stator. Ainsi, il se présente favorable du point de vue 
de la stabilité. En augmentant la valeur de K. on peut 
lever les instabilitês électriques, 


4.3 — Effet du contrôle de la fréquence 
de glissement, 8,., sur la stabilité 


On aditionne ou signal de la vitesse instantanée du 


rotor, 0, la pvaleur de la fréquence rotorique ou de 


E 
glissement, 8... qu'on veut et on obtient ainsi la valeur 
de la fréquence de commande de I'onduleur, c'est-á-dire, 
la fréquence statorique. 


FIG. 15 


L'action de ce contrôle se revéle importante du point 
de la stabilitê mécanique, En effet avec ce contrôle, la 
stabilité mécanique du systême présenté à la FIG. 15 
est toujours garantie puisque quelle conduit a un gain 
G < O dans le schéme-bloc de la FIG. 9. 

On remarque que G répresente la pente de la carac- 
téristique électromécanique. Le couple électromagnétique 


est donné par: 


& A: as 
Tesfilge — (24) 
(R,+A 8) 
0, M R, 
avec Ass. 
AZ, 


Sa representation graphique est: 


Te 


FIG. 16 
Le systême présente une caractéristique électromé- 


canique assez semblable à celle de la machine c.c. type 
série et n'a pas de problêmes d'instabiblité mécanique. 
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5 — CONCLUSIONS 


L'association onduleur de courant-machine asynchrone 
constitue un systéme ou les équations des circuits à con- 
tacts glissants sont d'application directe. 

Les hypothéses faites pour rendre possible [étude 
analytique de la stabilitê présentent des avantages. 
D'abord on voit qu'il y a deux origines pour les instabi- 
lités trouvés dans le systême, et on obtient des indica- 
tions utiles pour Iobtention des systêmes stables. 

La modification des paramétres du systême pour 
assurer la stabilitê electrique conduit à de mauvaises 
solutions du point de vue du rendement, donc la voie 
qui reste est [introduction de boucle de retroaction. Les 
contrôles du courant continu et de la fréquence de 
glissement, avec des actions complementaires, conduisent 
à la stabilité du systême, puisque si le premier élimine 
les instabilitês electriques, le second fait disparaitre les 
instabilités mécaniques. 


[1] Manuel Garrido — Les liaisons internes des ma- 
chines electriques à commutation. Técnica 
n.º 368, Maio 1967. 


[2] M. Garrido — Les équations générales des ma- 
chines electriques déduites de Ielectromagné- 
tisme. Revue E, Vol. 8, 1971. 


[3] N. Sawayi and MN. Sato — Steady-state and 
stability analysis of induction motor driven by 
corrent source inverter. I.E. E. E. —- Transactions 
on Industry Applications, Vol. IA— 13 nº 3. 
May/June 1977. 


[4] E. Cornell and T. Lipo —- Modeling and design 
of controlled current induction motor driven 
systems. |. E. E. E. — Transactions on Industry 
Applications. Vol. IA-— 13, ne 4, July/August 
1977. 


118 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


N. Sato, M. Ishida and N. Sawaki — Some 


characteristics of induction motor driven by 
current source inverter. |. E. E. E./l. A. S. Semi- 


conductor Power Converter Conference, 1977, 


E. Gudefin— Les machines à commutation 
electronique. Revue Générale de L'Electricité 


n.º 1, Janvier 1968. 


M. Garrido — Equations de la génératrice à 
redresseur. Etude du régime permanent. Revue 
E, Vol. Vil n.º 3, 1972. 


E. Cornell and T. Lipo — Design of controlled 
transfert function 
1. F. A E. 
Symposium on Control in Power Electronics and 


current a.c. drive using 


techniques. Conference Record of 


Electrical Drives, Dusseldorf, October 1974. 


B. Trannoy, M. Boidin and B. de Fornel — 
Investigation of the transient behaviour of a 
current-fed autopiloted asynchronous induction 
machine. PROC. |. E. E. Vol. 122 n.º 10, Octo- 
ber 1975, 


B. de Fornel and D. Noyes — Transient charac- 
terístic of an asynchronous machine fed from 
a static converter. PROC. |. E. E. Vol, 124 n.º 10, 
October 1977. 


T. Lipo and E. Comell — State-variable steady- 
state analysis in a controlled current induction 
motor drive. |. E. E. E. — Transactions on Indus- 
try Applications. Vol. IA— 11 nº 8, November/ 
[December 1975. 


TÉCNICA 454 


NUMERO 454 NOVEMBRO 1979 


ANO LIV 


VOLUME XLI 


Étude analytique et application des équatlons 
variationnelles de la génératrice à redresseur 


RESUMO 


São introduzidas as principais constantes de tempo 
que intervêm no funcionamento variacional dos geradores 
de rectificação e justificados vários níveis possíveis de 
simplificação. Os resultados obtidos são aplicados ao 
estudo da função de transferência das excitatrizes sem 
escovas e ao estudo da estabilidade do motor de indução 
alimentado por ondulador de corrente. 


1 — INTRODUCTION 


Les éqguations de la génératrice à courant alternatif 
suivie d'un redresseur à diodes, que nous appellerons 
génératrice à redresseur, sont aujourd'hui bien connues. 
Elles ont fait Vobjet d'un certain nombre d'études, soit 
dans le cas général des machines polyphasees, à póles 
saillants, en polygone ou en étoile, avec redressement 
simple ou en pont, à coupplage sinusoidal ou non [1], 
[2], soit dans le cas particuliegr d'une machine triphasée, 
à poles lisses, avec couplage sinsoidal et redressement 
en pont [3]. 

Ce cas particulier présente l'intéret théorique de 


mettre en évidence, d'une façon simple, les principales 


proprietés de ce type de systêmes et il est important en 
pratique parce-que beaucoup de réalisations utilisent, 
effectivement, des machines triphasées, à poles lisses, 
avec redressement en pont, 


Ce sera, par exemple, le cas de beaucoup d'excita- 
trices à diodes tournantes et encore, bien qu'indirectement 
(nous le verrons plus loin), le cas des moteurs à induction 
alimentés par un onduleur de courant, 

Pour létude de problêmes de stabilitê ou de régu- 
lation Ffaisant intervenir ces systêmes il est donc utile 
d'établir les équations de petits signaux ou équations 
variationnelles de la génératrice triphasée, à poles lisses, 
suivie d'un pont redresseur. 


Ces équations ont déjã été écrites et confirmées 
experimentalment en vue de son application aux excita- 
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SUMMARY 


This paper introduces the main time constants that 
influence the small signal behaviour of the rectified out- 
put AC generators and presents some reasonable simpli- 
fications. The results so obtained are used to study the 
transfert function of the brushless exciter and to study 
the stability of the induction motor driven by current 
source inverter, 


trices, On a ainsi constaté que, pour les transitoires 
purement electriques, on obtient un systêéme du trosiême 
ordre [4], [5], mais qu'il est raisonnable d'utiliser, pour 
la plupart des applications comme excitatrice, un modéle 
simplifié du second ordre ou même du primier ordre [6], 
[7]. 

Dans les travaux précedents, les modéles simplifiés 
ont été obtenus soit par des considérations à priori 
(négliger des termes diy/dt) soit par des considérations 
à posteriori (critêre des pôles dominants). 

Dans le présent travail, nous écrivons les équations 
de façon à justifier les simplifications possibles, dont 
certaines n'avaient pas été considerées auparavant. Pour 
cela, nous mettrons en évidence un certain nombre de 
constantes de temps caractéristiques du systême et leur 
rôle aux différents niveaux d'approximation, 

En vue d'applications plus générales, nous ajouterons 
aussi I'équation mécanique du systême. 

Les résultats obtenus seront ensuite appliqués à deux 
cas importantes en pratique. 

Le premier sera le cas classique de [excitatrice à 
diodes tournantes, ou on établira un schéma-bloc simple 
du systême. 

Le deuxiême cas, sera le moteur à Induction alimenté 
par onduleur de courant, ou on étudiera les problêmes de 
stabilité, qui difficultent I'utilisation de ce systême d'entrai- 
nement à vitesse variable, 
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2 — EQUATIONS VARIATIONNELLES DE LA GENERA- 
TRICE A REDRESSEUR SUR CHARGE INDUTIVE 


Considérons une gênératrice triphasée, à pôles lisses, 
avec couplage sinusoidal, suivie d'un pont redresseur et 
alimentant un circuit inductif de paramétres RL... 

Un tel systême peut être assimilé, 
po avec une bonne approximation [3], à 


TM : 


| 


une machine à courant continu, dont 
| les balais sont mobiles en fonction 
a de la charge, comme representé a 
la figure 1. 

Compte tenu de I'analogie du sys- 
tême avec une machine à courant 
continu, nous utiliserons, par la suite, 
les notations propres à ces machines, 
plutôt que celles habituelles pour les 
machines synchrones triphasées. 

FIG. 1 L'équation du circuit inducteur est 
léquation classique d'un circuit fili- 


Schéma de la géné- 
ratrice à redresseur forme 


dy 


U; = He Ip + Ei (1) 


ou W, est le flux couplé avec linducteur. 


Pour le circuit induit nous avons besoin de deux 
eêquations: l'une qui est I'équation éléctrique du circuit à 
commutation a e l'autre qui donne la position p des balais. 

L'êquation du circuit à commutation est, en négligeant 
la resistance intemne du circuit ou en la supposant inclue 
dans la resistance de charge R..: 


d Va w—p I Wom 
dt dp 


(2) 


al 


ou Woem est la coénergie magnétique du systême; , est 
le flux couplé avec le circuit induit; et u,=—R.ij— 


di 
Lee est la tension aux terminaux du circuit à com- 
t 


mutation. 


Dans le cas du redresseur à diodes, I'équation de 
position des balais est donnée par la condition: la ten- 
sion u, entre les balais est plus élevées que la tension u', 
entre une autre paire de balais imaginée à n'importe 
quelle position p'. Il en résulte, par application de I'équa- 
ton (2) à la paire de balais imaginée à la position p': 


[ao/ ay o OW, 
o =| Pl dt dai 9p' De, (3) 


ou () concerne les variables du circuit imaginé à la 
position p' et ou on a tenu compte du fait que 


Wo OW 
Lim = | - 
Mo dp' a dp 


Aux équations d'induit (2) et (3) on devrait ajouter les 
conditions d'irreversibilité du pont de diodes u,>0, 
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|, & 0, mais ces conditions ne vont pas intervenir dans 
l'analyse de petits signaux à faire par la suite, 
Finalement, on écrira encore l'êquation mécanique 
du rotor à commutation, qui est, en négligeant le frot- 
tement: 
d ty dW 


É me pommuceççms 4 
dt dp (8) 


ou ] represente le moment d'inertie et C le couple 
appliquê. 

La validitê du systêéme d'equations (1) a (4) a eté 
confirmée experimentalement dans d'autres travaux, y 
compris en régime variationnel. 

C'est pourquoi dans le présent travail nous les 
prendrons comme équations de départ, sans plus de 
soucis pour montrer leur validitê et nous bornerons à les 
transformer ou simplifier en vue de leur utilisation dans 
Vétude de certains problêémes pratiques qui nous inte- 
ressent. 

Introduisons les selfs L,, Lr et la mutuelle inducteur/ 
finduit M=—M, cosp. 

Le flux couplé avec le circuit induit est donné par 
U=Li—M,cospir et le flux couplê avec un circuit 


identique en quadrature (p, =p + a) est donné par 


Va = Mo sinp TR 

Nous allons utiliser Vag comme variable auxiliaire et 
profiter du fait que, pour les couplages sinusoidaux et 
aux fuites homopolaires prés, on a: 


dy dy dW 
— a ut — = ia ua (5) 
dp dp dp 


ou le symbol () concerne toujours la paire de balais à 
la position p”, 


L'ensemble d'équations (1) à (4) devient alors: 


dy 
re Elybo fi 6-a) 
í EE FF (6-a) 
E tu= md (6-b 
=—— + (v-p)y, E 
dt ' 
dy . 
SR sa (6-c) 
Cu ELO. (6-d) 
dt 


Nous allons maintenant écrire les équations varia- 
tionnelles, en linêarisant les équations precédentes pour 
des petites variations autour d'un état d'équilibre quel- 
conque, défini par: 


Ufo = Rrlto (7-a) 
Uso = Wo Vago (7-b) 
O = Vo (7-c) 
C, =— lo Yago (7-d) 
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On obtient, sans difficulté, les équations 


A ur = REA ip+p Ay (8-a) 


á Us, = Pp Á UR + Us à Dag + Vago à W— Wags páp (8-b) 


O = pá Vag, à b, (8-c) 
AC=]paát A Wig—Vago À la (8-d) 
d 


ou on a fait 


dese 


dt 


Nous allons réduire les variables par rapport à l'état 
d'équilibre du systême, en utilisant comme grandeurs de 
base les suivantes: 


Ufo» lo Leito, lao» Vão» Wago: Wo, Mot, C, (9) 
Si, en plus, nous utilisons l'équation (8-c) por élimi- 


ner Ay, en fonction de AW de systême d'équa- 
tions (B) se réduit à: 


ao 


Au=AS+T,p' ay (10-a) 
Auj=(1+pDAabiç+ AwW'—p'Ap (10-b) 
AC =Fp'aw + ap +ai, (10-c) 
u | ! 
ou on a u= É | = - E 
f R 
Ufo Io f 
“ Pp “ Wi * La * Vag 
pes ps o Mme; Nf im 
Wy Ly fo Uão =a Vago 
2 
to CG ] | 
(ts ma Gts, Peso (fo — 
Ud, C, la Vaio has 


Pour décrire complétement le systême, nous devons 
ajouter à ces équations les équations qui relient les flux 
aux courants et à l'angle p ainsi que équation de char- 
ge, soit (toujours en éliminant 41, par I'équation (8-c), 


AV =p Ad): 


Aba = 4 + cotgp,AP (11-a) 


* LO Gai co UR : x É O] 
pP'Aba=—cotgp, A + Ap + cotgp, AL, (11-b) 


Ay; = Aij-—Kcos2p, AÍ, + Ksinp,cosp Ap (11-c) 


Au=(1+Tep)Ai, (11-d) 


ou K est le coeficiente de couplage entre I'inducteur et 
Hinduit et Tº = ww, L./R. est la constante de temps de 
la charge. 


Nous disponsons donc, au total, d'un systême de 


7 équations, avec 7 sorties Lip Vi Vito) Pruç et 
2 entrées u,, Cº. 
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Pour les applications envisagées nous allons retenir 
comme sorties |; et ww" et éliminer les autres variables. 
Néanmoins, par commodité d'écriture, nous garderons 
dans les équations, comme variable intermédiaire, le flux 
en quadrature W'.. 


Dorénavant, nous allons abandoner le sybol (*) pour 
les variables réduites, étan entendu que, sauf indication 
en contraire, nous travaillerons toujours avec des gran- 
deurs réduites, 


21 — Modéle général, du quatriême ordre 


Pour éliminer les variables ir, P nous partons des 
équations (11-a, b) qui donnent 
Ap =sinpocos po [lt + tgro Pp) Sb di] (128) 


dir = (sin?p,—sin p, cos P, p) AV + cos?pçãi, (12-b) 


En introduisant les équations (12) ainsi que les équa- 
tions (11-c,;d) en (10), on obtient le systême: 


Au,=[(1+Tp)-sinZo,(1+T; pJA i,+ (13-a) 
+[sinZe, (1 +Tf p)+TEp O +TE PIA Via 
O = Tp+sino,(1+T. plAi,— (13-b) 
— [sino ,(1+T, pD)+T2pJ] Ab sin? p, Am 
AC=AL+AWItIPAW (13-c) 


Dans ces équations variationnelles figurent les cons- 
tantes de temps éléctriques suivantes, caractéristiques du 
systême et de la charge: 


Us 
T = (1—K2 cos? p,) (14-a) 
E 
vw. L; 
"o a (1 +K2 cos? Pu) (14-b) 
| A, 
ty 
T'=— É (-Ksinte) a 
É 
Ud, E 
T' = - sin? Po=—sin Po COS P, (14-d) 
q] 
to. L. 
T o o sin? p, (14-e) 


E 


ou la deuxiême égalité (14-d) provient des équations de 
régime permanent (7). 


La comparaison des constantes de temps précedentes 
peut fournir des indications utiles pour la simplification 
des équations variationnelles du systême, comme nous 
verrons. 


2.2 — Modéle du troisiême ordre 


Nous pouvons abaisser l'ordre du modele en négli- 


geant la derivée 4, (ainsi que W,) dans les équations 
d'induit, par l'analogie avec ce que est fait dans I'étude 
classique de l'alternateur ou de la machine à courant 
continu, 

Ainsi, si on élimine pAW,, [par I'équation d'induit 
(8-c)] on obtient àw, = 0. L'expression (11-b) devient 
nulle et sa cobinaison avec (11-a) fournit les nouvelles 
relations: 


Ap =sinpocos Po (A Vg— A ij) (15-a) 
A ip = sin?p, Ab + cos? p, Ai, (15-b) 


L'introduction de ces résultats fournit le nouveau 
systême d'équations simplifiées: 


Aup=[0+T p)-sin?e,(+Tr p)JAi+ (16-3) 
+sin po (1+TE PJA Va 
O =[T.p+sino,(1+T.p)JAi, 
— sine, (1+ Ti PDAb—sinp,A (16-b) 
AC=AL+HAba+IpAM (16-c) 


Ce systême d'équations est du troisiême ordre, 
comme annoncé au départ, 

Malgré le fait d'avoir négligé 1), la constante de 
tempes T. figure dans les équations concernant I'induit 
et cela parce que les balais du systême restent mobiles 
dans le modeéle, 

On peut faire disparaitre cette constante de temps 
si, en plus de négliger b,, on néglige » dans les êqua- 
tions d'induit, c'estâ-dire, le terme p dp dans I'équa- 
tion (B-b). 

Cette simplification supplémentaire est équivalente à 
supposer by = O et 4, = 0, dans un modéle de Park de 
la génératrice à redresseur [1]. Cette hypothése est juste- 
ment celle qui a été adoptée dans la référence [6] pour 
simplifier les équations variationnelles de la génératrice 
à redresseur., 

Avec notre écriture des équations, on obtient dans 
ce cas le systême suivant: 


Aup= [(1+Trp)-sinZp, (1+Tf p)] Ai,+ 
+sinZo, (1+TF pJA Va (17-a) 
O =(T.p+sinZo)A i—sin?p, A Wg—sinZpçã to (17-b) 
AC=AL+AVatIpAs (17-c) 


ou intervient une seule constante de temps à I'induit 
quis est T. 
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2.3 — Modeéle du second ordre 


Lorsque les effets de la charge eélectrique sont 
dominantes, on peut simplifier davantage le modéle en 
négligeant, en plus des derivées W,, Vig la derivée du 


flux à I'inducteur Jr. 


Les équations du systême deviennent alors 
Aur=(I-sinZp,) A is+sinZp, A Va (18-a) 
O = [T, p+sin?p,(1+T; DJ Ai,— 
—sin?p,(1+T, PD) A Va—sin?p, A to (18-b) 
A=AL+AbgtipAt (18-c) 


La constante de temps T, pourra encore être eliminé 
en imposant à [induit p= 0 c'está-dire, en faisant 


da =W = 0. 
On obtient alors 
duç=(1—sin?p) Ai, + sin?p, A La (19-a) 
O = (Tep+sin?p,) A i—sin?p,A Lag—sin?pção  (19-b) 
AC=d,+áv,+IpAm (19-c) 


Dans ce modéle la machine elle-même est réduite à 
à un simple gain, les élêments ênergétiques étant, uni- 
quement, la charge électrique et la partie mécanique. 


2.4 — Remarques sur les differents modéles 


La comparaison entre les différents modéles, établis 


par des simples considérations sur les termes W ou p, 
peut-étre faite à partir des valeurs des constantes de 
temps du systême. 


On voit que le passage au modéle du troisiême ordre 
est fait en négligeant des termes ou figure la constante 
de temps T. Dans un premier temps on néglige les 
termes du second ordre des équations (13), ainsi que 
la constante T. comparé à T; T, ; dans un second 
temps, on élimine totalement Es ce qui implique, em 
particulier, qu'elle soit négligeable comparé à T.. Le 
passage au modéle du second ordre implique, en plus 
simplifications précedentes, la possibilité de négliger les 
constantes de temps T, T; de l'inducteur, face à la 
constante de temps T. de la charge électrique et à 


linertie mécanique du systême, 


Cette derniére simplification sera normalment peu 
satisfaisante pour les génératrices à redresseurs, mais 
elle sera acceptable pour les systêmes d'excitation stati- 
que (qui sont assimilables à des génératrices à redres- 
seurs [8]) à condition de pouvoir négliger les fuites du 
transformateur. 
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Par contre, l'approximation du troisiême ordre parait 
três raisonnable pour la plupart des applications, de 
sorte que nous allons la traiter plus en détail et I'utiliser 
ensuite pour deux exemples d'application importants en 
pratique, qui sont |excitation d'alternateurs par diodes 
tournantes et l'alimentation de moteurs à induction par 
onduleur de courant, 


3 — ETUDE ET APPLICATIONS DU MODÉLE 
SIMPLIFIÉ DU TROISIÊME ORDRE 


Dans lI'utilisation de ce modéle nous avons interêt 

à introduire la variable (à b,,;— di.) qui posséde un 

rapport étroit avec la variable 4p par I'équation (15-a), 
que nous réecrivons: 

dp =sinp,cosp, (Abq— Ai) (15-a) 


Si on introduit la variable (AW, — di) les équa- 
tions (16) deviennent 


Su=(1+Tp) Si +sine (1+T/ pn) (AV —Ai) (20-a) 
O =Tpái—( +Tp)(Ab—Si)—Amw (20-b) 
AC=24,+(Abg— di)+HIpão  (20-c) 


La deuxiême équation a été divisée par sin?p,, ce qui 
fait apparaitre explicitement la constanet de temps de la 
WD Lo 


Li 


charge Ta= 


Lorsque T', €T'. on pourra pousser plus loin la 
simplification en négligeant tout à fait la constante de 
temps T',, comme déjá vu dans le paragraphe precédent. 


FIG. 2 


Schéma-bloc de la génératrice à redresseur (modêle d'ordre 3) 
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Le systême d'équations (20) peut-être representé par 
un schéma-bloc relativement simple 

Ce schéma-bloc est susceptible de transformation 
qui le rend plus commode, mais nous preférons le garder 
tel quel, puisque les exemples à considérer par la suite 
permettent certaines simplifications assez immédiates sur 
le schéma de la Fig. 2. 

Ainsi, par la suite, nous poursuivons lanalyse du 
modéle à propos de deux exemples concrets, 


3.1 — Excitatrices à diodes tournantes 


Dans cette application, une génératrice à redresseur, 
à induit tournant monté sur ['arbre d'un alternateur, fonc- 
tionne comme excitatrice, le courant i, étant alors le 
courant d'excitation de |'alternateur, 


cad À 

q 3i(Hmps (CE 
1 ! 
1) 


a O OO 


excitatrice alternateur 


FIG. 3 


Excitation par diodes tournantes 


L'entrée principale de l'excitatrice est la tension ur 
et la sortie principale le courant redressé i,. 

Pour les excitatrices utilisées la constante de temps 
d'inducteur Tr=tw, L/R, est, en général, beaucoup plus 
grande que la constante de temps du circuit induit 
Ti=w,L./R. (ou figure toute la resistence de charge R,), 
de sorte qu'une approximation du troisiême ordre est 
raisonnable. 

Quelques simplifications supplêmentaires sont encore 
possibles. 

Ainsi, dans les excitatrices à diodes tournantes la 
vitesse est imposée, pratiquement, par le seul alternateur 
principal, De la sorte, àw devient une entrée, ce qui 
permet d'ignorer I'équation mécanique. 

D'autre part, T. peut-être assimilé à la constante de 
temps transitoire en charge T',, de l'alternateur principal, 
en général beaucoup plus grande que T',. On se trouve 
donc dans un cas ou on pourra, tout à fait, négliger 
cette derniêre. 

Avec ces deux approximations supplêmentaires, les 
équations (20) se réduisent à: 
du =(1+Tp)di+sine(1+T; (Av, —Ai) (21-8) 
Sto =TapSt(Ab—A io) (21-b) 


ou ne figure pas la constante de temps T', et dou est 
disparue [équation meécanique. 
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On obtient donc un systême d'équations purement 
électriques, qui est alors du second ordre, Son schéma- 
-bloc est particuligrement simple: 


FIG. 4 


Schéma-bloc de |'excitatrice à diodes tournantes 


Dans la réference [6] on étudie un modéle du secon? 
ordre de la génératrice à diodes tournantes, qui est 
équivalent à celui que nous venons de décrire. On y fait 
sa comparaison au modéle complet, dans plusieurs cas 
concrets, et on constate qu'il donne satisfaction dans les 
situations réelles envisageables. 

En poursuivant l'analyse du modéle, nous pouvons 
écrire les fonctions de transfert entre les deux entrées 
et la sortie, On a 


Ai, 1 
= 11 STE (22a) 
du 1+(Ti+sinp,Tu)p+sin?p, Tal p? 
et 
ai, sin? p, (1+Tp p) 
(22-b) 


Ao  1+(T,+sin2 Poa) pAsin?P, Ta TE p? 


La fonction de transfert (22-b) correspond à un terme 
d'amortissement supplémentaire dans I'alternateur princi- 
pal, normalement négligeable, 

En effet, si on utilise des valeurs réduites, de base 
P, = Coto, pour la puissance d'amortissement Pp de 
Halternateur principal, on a: 


A Pp AP, Ai, 


(23) 
At) Ai, AW 


Pour un angle de puissance ô constant, on obtient, 
toujours en valeurs résuites: 


AP 
ça À (24) 
A is 
On peut donc écrire 
A Pp sin Po (1+Tr p) cê 
dw 1+Ti+sinp Ta) p+sinp, Ta Tr p? 
Le gain statique, en valeurs non-réduites, est 
Po 
G =— sin? p, = C sin? p, (26) 
W,, 
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Ce gain est nul pour p, =T (géneratrice à redres- 

TE 

seur proche du court-circuit) et maximum pour a 
(gênératrice à redresseur proche du vide). 

On voit que le terme d'amortissement, de l'ordre 


de C.. est d'un ordre inférieur en aux autres termes 


RR 
E) 


de I'équation mécanique approchée: 


ps e” Ô 21 
dE Te E id (27) 


On pourra donc l'ignorer, normalement. 


Considérons maintenant les transitoires purement 
électriques (Ato = 0) et la fonction de transfert (22-a), 
que nous allons étudier avec un peu plus de détail, en 
vue de la modélisation des boucles de de regulation de 
tension, 

Ainsi, voyons comment la nature de la réponse 
dépend des constantes de |'alternateur principal, en parti- 
culier de T',, et des paramétres du systêmes d'excita- 
tion, en particulier de p, = arctg (—R./w, Lj). 

Commeçons par considérer le fonctionnement pour 
un angle de régime fixé, proche de la puissance maxi- 
male de la génératrice: fp, = 135º, 

Voyons l'influence du paramêtre T'. = sin?p, T'a, sur 
la nature de la réponse électrique., 

L'amortissement critique du systême sera obtenu 
pour les valeurs 


T=2T7[ 1= Fioação —T (28) 


qui sont toujours réelles et positives compte tenu des 
expréssions de T"p et T'. 

La réponse sera oscillante pour des valeurs de T'. 
comprises entre les deux valeurs critiques précédentes 


(To < Tia < (T'o)z (29) 


Pour un angle p, = 135º et en admettant un bon 
couplage de la génératrice KZ1, on a T,=3T;= 
= 1.5T,. Ceci donne 


(Th =0 TJ =5T; (30) 


1 
Donc, compte tenu de T', = sin Po a ET la 


rêponse sera oscillante pour 


Tu <1OTE (Po = 1359) (31) 

Celui-ci sera le cas normal pour les systêmes d'exci- 
tation par diodes tournantes, 

En dehors de cette zone normale d'utilization, ou la 
réponse est oscilante, le systême présentera des répon- 
ses amorties qui pourront, à la limite, être représentées 
par des systêmes du premier ordre. 
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Ainsi, pour T'. <(T'.),) on peut supposer T',.— 0 et 
on obtient alors un systême du premier ordre, avec une 
constante de temps T'. 

Pour T'. > (T'.)) on peut admettre que la charge est 
dominante et poser T;— O. On obtiendra alors un sys- 
têéme du premier ordre, qui correspond simplement au 
circuit inducteur de [alternateur. 

Considérons maintenant une valeur fixé pour T'4, et 
voyons l'influence de l'angle de régime p, sur le compor- 
tement du systême. 

Si le systême est conçu pour travailler avec des 
angles p, proches de T/2 (cest-â-dire, génératrice à 
redresseur proche de la marche à vide), on obtient 
T=STET. 

Dans ces circonstances les deux valeurs critiques 
sont égales à T, et le systême est toujours amorti, pré- 
sentant les constantes de temps Tr et Ta. 

On obtient ainsi, compte tenu de sine, = 1: 


dia 1 
du +HTEPD( + Tap) 


(32) 


Ce modéle est identique à celui d'une excitatrice 
classique à collecteur, ce qui résulte justement du fait que 
les balais sont fixés dans une position proche de T/2. 

Si le systême est conçu pour travailler proche du 
court-circuit (p, = 7), on obtient, pour T',, fini, une cons- 
tante T.=sinfp,T'a =0. Dans ces circonstances le 
svstême devient un systêéme du premier ordre, caracterise 
par la constante de temps Tp=(1—-K9Tp=cT, qui est 
la constante de temps transitoire longitudinale en court- 
«circuit de l'excitatrice 


O ai (33) 
du 1+0Trp 


Dans ce cas il y a une predominance de |'excitatrice 
sur Hinducteur de [alternateur principal, ce qui nest pas 
une situation normale, 

Remarquons, neéanmoins, que la situation p,=T 
pourra se présenter avec les machines criogêniques, mais 
dans ce cas, déjã traité dans la reférence [5], on ne 
pourra pas négliger T', à cause de I'augmentation de T'. 
On a méme constaté |'existance d'une possible instabilité 
électrique du systême, qui ne pourra pas être étudiée 
avec le modéle simplifié que nous sommes en train de 
considérer. 

En conclusion de l'étude, nous dirons que, dans les 
situations normales, le transfert entre la tension d'entrée 
de [excitatrice u, et le courant d'excitation de I'alterna- 
teur i, est donné par une fonction du second ordre, avec 
réponse normalement oscillante. 


L approximation du premier ordre, parfois adoptée 
pour ce type “excitatrice, ne sera correcte que pour 
certains cas extrêmes, moins fréquents en pratique. 
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3.2 — Stabilité du moteur à induction alimenté 
par un onduleur de courant 


Nous allons considérer un moteur à induction alimenté 
par un onduleur de courant, systême représenté dans la 
Fig. 5. 


Moteur à induction alimenté par onduleur de courant 


Le systêéme présente un certain nombre de problêmes 
de stabilité, étudiés avec détail, pour des cas concrets, 
par exemples dans les reférences [9], [10], [11]. 

Les études precédentes, effectuées par voie numéri- 
que sur les équations complétes du systême, ne per- 
mettent pas de dégager des conditions analytiques simples 
de stabilité, ni de voir três clairement l'influence des 
différents paramétres, 

Une premiére tentative d'étude analytique du pro- 
blême est faite dans la reférence [12], moyennant un 
certain nombre d'hypothéses simplificatrices permettant de 
dégager les principaux phénomêénes en jeu, 

Nous allons reprendre le problême de la stabilitê du 
systême à partir des équations variationnelles de la 
génératrice à redresseur, 

Pour cela, nous utiliserons un résultat établi dans la 
reférence [12], qui montre que le moteur à induction 
alimenté par onduleur de courant est équivalent à une 
génératrice à redresseur, mais dont [inducteur est aussi 
un circuit à commutation, 


Pour le voir, nous allons consi- 
dérer le stator du moteur comme 
linducteur et assimiler l'entrée (coté 
continu) de I'onduleur à une paire de 
balais, comme pour le redresseur en 
pont. 

Ces balais tournent avec une fré- 
quence angulaire t),, imposée par la 
fréquence de commande de [onduleur. 

Maintenant nous prendrons comme 
reférentiel celui des balais du stator, 
de sorte que ces balais deviennent 
fixes et c'est le stator qui tourne 
avec une frequence tu. 

Ce referentiel sera celui de l'inducteur de la génê- 
ratrice à redresseur, qui est donc un circuit à commu- 
tation. 

L'induit de la génératrice sera formé par le rotor 
du moteur à induction, tournant par rapport à linducteur 
avec une vitesse ty = gw, ou g est le glissement du 
moteur à induction, 

En supposant le circuit rotorique du moteur comme 
polyphasê et à couplage sinusoidal, nous pouvons le 
réduire à un circuit à commutation, par la transformatión 
«8-6 [13]. 


FIGI 6 


Modéle du moteur à 

induction alimente 

par un onduleur de 
courant 
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Nous obtenons donc un induit de génératrice à 
redresseur, avec la seule particularité qu'il se trouve en 
court-circuit, Ce fait signifie que T', = 0, ce qui permettra 
de simplifier les équations. 

Les équations de la génératrice à redresseur seront 
donc applicables au systême, à condition d'ajouter à 
l'équation d'inducteur un terme u de tension induite par 
commutation, comme le dernier terme de |'équation (2). 
Compte tenu de p. = 0, on a, valeurs non réduites: 


ww, AWom 
f dp, 


(34) 


Avec des couplages sinusoidaux, la mutuelle induc- 
teur/induit est de la forme M=—M,cos (p—p.), de 
sorte que, pour p, = O on obtient 


ur=—tw. M,sinpi, (35) 
La variation du”, est donc, en valeurs non réduites: 
> u'f EE mi), Mo SinP, + la— U), My CosP, lag 4 P (36) 


Ce qui donne, en valeurs réduites: 


duf=— Es Ti(sinp,cosp, di + cos?p do) (37) 


Es 


Cette tension supplêmentaire à ['inducteur peut donc 
se traduire par deux simples réalimentations de la géné- 
ratrice à redresseur à partir des sorties i, et Pp. 

Pour introduir ces réalimentations, nous partirons 
directement du modéle du troisiême ordre, qui est valable 
si on fait |'hypothése que la constante de temps du circuit 
inducteur est plus importante que cella du circuit induit, 
ce qui sera possible grace à la presence de la self 
extérieure de lissage. 

Compte tenu du fait que T, est nulle, nous utiliserons 
le premier modéle du troisiêéme ordre, ou figure la 
contante de temps T'.. 

Le schéma-bloc de la Fig. 2 devient ainsi, en utilisant 
toujours des grandeurs réduites: 


FIG. 7 


Schêma-bloc du moteur à induction alimenté par un onduleur 
de courant 
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Les constantes % et q sont données par 


1 
q = — K2 T sinp, cosp, (38-a) 
Jo 
1 cos'p, : 
& =— KT; — = & cotg?p, (38-b) 
Jo sinp, 


la 
Nous allons nous intéresser au transfert é Une 


transformation simple nous donne alors le schéma-bloc 
suivant, plus simples: 


| sinêP, 
[+ Tapldp-1 


[=> 
- 


FIG. 8 
Schéma-bloc du moteur à induction alimenté par onduleur de courant 


Le polynome caractéristique du systême sobtient 
facilement: 


AI Te LITE = 0) TElpés (39) 
+[2 sino, T;—T; +] (1—2)] p+I2 sino, (1-04)—(1—0)] 
soit, en éliminant & en fonction de dy: 

A =]TTEpê+ITE+(I—a)T]p?+ (40) 
+ [2 sino, T;—T; + J(1—a,))]p—cos20,(1 + «,) 
Il s'agit d'un polynome du troisiême ordre à = a, p'+ 


+ajp?+asp+as que nous allons traiter par le critêre de 
Routh-Hurwitz. 


À fin de faciliter l'étude, nous transformerons 
dabords É, pour avoir une expression plus commode, 
ou figurent les paramétres de |'inducteur et de I|'induit. 


to To 19, 
Si on tient compte de T;= — E = So tia cê on obtient 
Re af 
ty, Mo ã p) p (41) 
& = sin 5 
2L, R 


L'angle p, variant entre 90º et 180º, la fonction q, 
est toujours negative, 


Nous avons intêret à représenter les fonctions 1—ua, 
et 1 +, qui interviennent dans le polynome caractéris- 
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tique, parce que cela nous facilitera lanalyse qui doit 
suivre: 


FIG. 9 


Passons maintenant à l'étude de la stabilité du sys- 
têéme, en analysant le polynome caracteristique, 


3.2.1 — Stabilité apériodique 


Nous remarquons d'abords l'existance de limites de 
stabilité apériodique, imposées par la condition a; > 0, soit 


— cos2p,(1+8)2>0 (42) 
Une limite de stabilité est atteinte pour 
cos2p, =0 (43) 
soit p, = 135º. 


Compte tenu du fait que tgP,=—— =— 
Us 


la limite précedente est obtenue pour 


(44) 


Étudions le comportement du systême à droite et à 
gauche de ce point. 

Nous devons distinguer deux cas, selon le signe 
de (1 + &ias” 


two, Mo? 


————. < 1 
2L, A 


a) 


Dens ce cas (1 + &)ja5º > 0, de sorte que le terme 
ay=—cos2p, (1 + m) est positif si p, <. 135º et néga- 
tif si Pp, > 135º, 

Pour que le systême soit stable il faut donc que 
Hon ait p, < 135º, soit g,<R./wL.. 

Cette situation correspond à avoir une caractéristique 
electromécanique du type suivant 
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Fil 
dl Resulla O 9% 
180" re 135º — 90ºp, 
FIG. 10 
s Ms 


Limite de stabilitê apériodique pour <1 
2 


Ls H 

Avec un couple extérieur C, constant, le fonction- 
nement ne pourra donc être stable qu'au delã du maximum 
de couple. 


vw Mp. 
b) mm > 
2L. Rs 

Dans ce cas, on a (1 + 0); <0, comme on voit 
dans la figure 9, 

Pour que le systême soit stable autour du glissement 
critique g, = R./w L, il faut donc que Pp, > 135º; il sera 
instable pour Pp, <. 135º. 

Cette situation correspond à une caractéristique elec- 
tromécanique du type suivant 


FIG. 11 


w Mt 


Stabilité apériodique pour 


La condition (1 + €) < O reste valable jusqu'ã deux 
limites p;, po ou (1 + €) = O. Pour déterminer ces limites 
nous introduisons |'expression (41) de %, et nous utilisons 
le fait que 


2 cotg P, 2 tg P, 
BiNnZ=—"——— = ——————— 
1 +cotgip, 1 +tg?p, 
Nous obtenons, soit Iéquation: 
Ud, “1 
cotg? p, + “cotgp,+1=0 (45) 
F = 
soit Véquation: 
ty, Mo? 
tg” po + tgp, +1=0 (46) 
f a 


Le systême possêde donc deux solutions &;, py telles 
que tg 2) = cotg Pp; soit, compte tenu de 90º < p, < 1809, 
telles que Pp, + py = 270º, 
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Les glissements correspondants sont donnés par 


= ——————s 
MR. | 2. Re Nº 
2L*A UM, Ma? 


qui est le résultat proposé dans la reférence [14] et 
retrouvé dans les reférences [9], [10], [111]. 

Au delã des limites g; et g; foumies par (47), les 
conditions de stabilitê apériodique changent, comme repre- 
senté sur la figure 11. 


WU, M,? 
Si ——— > 1, comme il arrive souvent en prati- 
2L.R 
que, on obtient 
Mo? Fla T, 
= ss = K— 48 
9 Ja ZA, T, (48) 


Remarquons que dans ce cas g, peut-être supérieur 
à 1 et que les angles p, correspondants à g, et 9, sont 
alors py = 90º et py = 180º, 

Le couple est le même pour les deux limites g, et 92 
et il peut être determiné par l'étude de la caractéristique 
electromécanique [14] On obtient alors 


up 


-— 49 
E (49) 


Nous utiliserons ce résultat plus loin, 


3.2.2 — Stabilité périodique 


Nous allons maintenant analyser les autres conditions 
de stabilitê du systême, qui donnent origine aux instabi- 
litês oscillatoires observées en pratique. 

Compte tenu du fait que x, est négatif, le terme a; 
est toujours positif. 

Il nous restera donc à analyser les conditions a; > O 
et aa—a,a; >0. 

En réalité, il nous suffirait d'analyser la deuxiême 
condition, puisqu'elle implique la condition a, > 0. 

Néanmoins, nous analyserons aussi, directement, la 
condition a; > 0, puisque cela nous donne les résultats 
qui seraient obtenus avec |'approximation supplêmentaire 
du modéle qui consiste à négliger tout à fait la constante 
de temps T,. Parce-que T.= 0, cela nous conduirait à 
un modéle du second ordre. 

Nous allons séparer les deux cas typiques identifiés 


& 


2 
O : ] , 
selon la valeur de ——— et étudier ce qui se passe 


a CM 


dans les zones stables, apériodiquement, de ces deux cas, 


tos Mo? | 
3.2.214—- Cas -————— < 1 
2L,R, 


Ce cas correspond à la caractéristique electromé- 
canique de la Fig. 10 et nous savons déjà que le systême 
ne sera stable aperiodiquement que dans la zone ol 
Do S 135º, 
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Dans cette zone, la condition a; > 0, soit 
2 sin? Ps T'—T, +J(I-2)>0 (50) 


est toujours verifiée, 
En effet, pour p, < 135º on a 2sin?p, > 1. Puisque 
T'p>Tret q < 0 la somme sera done positive. 
Considérons la condition a, ay— a, as > 0, soit, aprês 
division par J: 


T+) T [sine TT + —)] + 
+TTicos2(1+]4)>0 (51) 


Le dernier terme étant négatif dans la zone considerée 
de stabilité aperiodique, la condition pourrait ne pas être 
satisfaite. En réalité nous allons voir qu elle [est toujours 
compte tenu de |'hypothêése T, > T, deja faite. 

Puisque (1 — 64) T', est positif, la condition (51) sera 
verifige si on verifie la condition plus sevêre que Ion 
obtient en négligeant (1 — ey) T'.. 

Si on divise encore par T;, on a: 


2sin?p,Tº--THI(i—a)+T;cos28,(1+%)>0 (53 


Si on developpe T” et T; et on utilise I'égalité 
i—2sin?p, = cos2p,, la relation (52) devient 


cos2 po [sin e, Tilt + )—T| + 
+ Kcos?p, (1 +2sine)Te+rI](i—)>0 (53) 


Compte tenu du fait que 1+a < 1 etdecos2p,<O, 
dans la zone considerée, le premier terme sera toujours 
positif si 

ST (54) 


hypothêse que nous avons dejá admis, 
Les autres termes de (53) étant positifs, nous con- 
cluons que la condition sera donc satisfaite, 


En conclusion, le systême sera stable periodique et 
apériodiquement, pour p, < 135º, si la condition T, < Te 
est satisfaite, ce qui est une hypothése acceptable pour 
le moteur à induction alimenté par onduleur de courant. 


o, M,2 
St 


3252 Os: sarna 
2L,R, 


Dans ce cas nous avons identifié deux zones de 
stabilité apériodique, que nous devons maintenant étudier. 

Dans la zone des petits glissements, ou Pp, <. 135º, 
nous pouvons utiliser les arguments du cas precédent, 
et conclure que le systême sera periodiquement stable, 
si on vérifie la condition 


Tu Te (55) 


Les problémes peuvent apparaitre dans l'autre zone 
de stabilité apériodique, obtenue pour p, > 135º, 
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En ce qui conceme la condition a, > 0, soit 
2sinp, TETE +I(i— a) >0 (56) 


nous voyons qu'elle pourra ne pas étre verifige, parce 
que, maintenant, on a 2sin?p, < 1. 

Cette condition n'est pas três commode, de sorte que 
nous allons chercher une condition plus simple, bien 
qu'approximée. 

Compte tenu de la self extérieur, le couplage K? est 
três petit. Nous ne commettrons donc pas une grande 
erreur en remplaçant T"; par T' dans expression (56). 
ce qui donne une condition plus forte: 


J(1— a) > cos26,T'; (57) 


Cette condition est verifiée pour p, = 135º, 


Si & est important et compte tenu des expressions 
des valeurs réduites, nous pouvons conclure que le pre- 
mier membre croit, au départ, moins vite que le second, 
avec l'angle p, (dans le rapport cos” p,/cos 26). 

Ceci nous ammenera à chercher une condition suffi- 
sante de stabilitê pour toute la zone considerée, que 
sera obtenue en considérant le fonctionnement proche du 
couple maximum C. 

On peut alors faire cospH1, I-,=Z2 et 
TE=0T, ou c=1—K? est le coefficient de fuites. 
La condition (57) devient donc 


o = S Tm (58) 


que nous prendrons comme une condition suffisante pour 
toute la zone de fonctionnement considerêe, 


Celle-ci serait la conclusion d'une étude faite avec 
le modéle plus simplifié, ou on aurait posé T', = 0. 


Si on utilise le modéle plus complet, ou cette simpli- 
fication n'est pas faite, on doit encore considérer la 
condition plus sevêre, a, az—a, a; > 0, dejá écrite en (51): 


Me+(1—-a)T sine TETE (1—)] + 


+ TiTicos?p,(1+%)>0 (59) 


Cette condition est difficile d'analyser et nous cher- 
cherons aussi une autre condition, plus forte, qui la 
remplace. 

Ainsi, compte tenu du fait que Ion a TE <T", est 


que la condition (55) implique T', <T',, on peut écrire: 
(Q—m) 2sin?o, TT + II) + 
+T;cos2p,(1+mÁ)>0 (60) 


Soit, aprés quelques transformations simples: 


1i+a 


J(1—aq) > ( rem ) cos 2p,T' (61) 


— 1 
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Cette condition est plus sevêre que la condition (57), 
puisque le coefficient du second membre est compris 
entre 1 et 2, dans la zone considerée. 

Néanmois, si nous acceptons les arguments utilisés 
pour la condition précedente, selon lequels il fallait consi- 
dérer ce qui se passe du coté du couple maximum, le 
coefficient tends alors vers 1, se qui nous conduit de 
nouveau à la même condition (58), soit 


Im > Tim (62) 


Cette condition parait donc assez forte pour assurer 
la stabilité dans la zone des grands glissements, 


3.2,3 — Verification des résultats 


On peut résumer maintenant les résultats obtenus, 


tu, Mo 


Dans le cas — <. 1, correspondant à la Fig. 11, 


a" 
on arrive à la conclusion que le systême est periodique- 


ment stable, compte tenu du fait que I'on peut admettre 
comme normalement satisfaite la condition suffisante: 


Tas Te (63) 


s Mo 


> 1, correspondant à la Fig. 12, 


Dans le cas 


ar M 
en plus de la condition precédente, assurant la stabilité 
periodique dans la zone des petits glissements, on est 
arrivê à une autre condition dans la zone des grands 
glissements, à laquelle on a donné [expression appro- 
chée (62). 

Compte tenu du fait que cette expression est écrite 
en valeurs réduites, on préferera utilizer une expression 
en valeurs non-réduites, qui est: 


J E” q Tç (64) 
vu”  8RW2gm 


ou 9, est le glissement pour le couple maximum. 


i ) Mo l Ra 
Pour l'approximation g, = —— on pourra 
| E Re 


écrire 


o LL? 
= — ———— (65) 
up 8 t),? Mo? af Ra 


Dans la Fig. 13, nous comparons ce résultat à celul 
obtenu, nmumériquement, dans la reférence [11], pour un 
systême qui satisfait aux conditions de notre étude, en 


2 
M, Mo 


2L,R 


=63 gm=05 et 


particulier, Te/Ta = 49, 


Po (Gm) = 1710 
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FIG. 12 
Limite de stabilité pour un sistéme avec o = 0,9; 
T, = 0,68,; Re =921 wu =27X 30 s-! 
a) Etude numérique compléte du systême; 
b) Condition (64); c) Condition (56); 
dj) Condition (59). 


La condition de stabilité (65) montre que l'augmenta- 
tion de ), AR; R, et to, ont un effet favorable sur la stabi- 
lité periodique, tandis que l'augmentation de up et L, ont 
un effet défavorable. Ces résultats sont, justement, con- 
firmés par I'étude numérique faite dans la reférence [11]. 

La condition montre encore qu'une diminution de 

x 


= 1— 


et de L, sera favorable à la stabilité perio- 
al 
dique. 

Les résultats obtenus montrent que dans le cas con- 
sideré on pourrait même adopter le modéle plus simplifié, 
ou on néglige T',, de sorte que le systême s'étudie, prati- 
quement, comme un systême du second ordre, 

Mais l'etude de la stabilité pourra se compliquer dans 
le cas ou la condition (63) n'est pas satisfaite, I'utilisation 
du modéle du troisiême ordre lui même eétant alors peu 
acceptable, 

Dans la reférence [10] on considêre des cas oú la 
condition (63) nest pas verifige et les résultats obtenus, 
numériquement, laissent prévoir que les conditions de 
stabilité seront plus complexes que celles établies dans 
notre étude, le systême devenant pratiquement inutilisable 
sans boucles de stabilisation, 

Certaines conclusions pour le cas ou TE; =T, sont 
indiquées dans la reférence [12] et elles s'accordent 
encore, d'une façon générale, avec celles de [10]. 


4 — CONCLUSION 


Nous avons effectué une étude analytique des équa- 
tions variationnelles de la géênératrice à redresseur et 
établit plusieurs modéles simplifiês, justifiés soit par des 
approximations à priori (en négligeant I'effet de termes 


comme 


) soit par des approximations numériques à 


posteriori (en négligeant des constantes de temps), 
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D'une façon générale, il est apparu que ['essentiel 
du comportement du systéme est décrit par un modêle 
d'ordre réduit, comme il avait déjá été constaté dans les 
reférences [4], [6]. 

Nous avons ensuite utilisé un modéle simplifié pour 
étudier la fonction de transfert des excitatrices à diodes 
tournantes et la stabilité du moteur à induction alimenté 
par onduleur de courant. 

Pour les excitatrices à diodes tournantes, nous avons 
montré que leur réponse est normalement oscillante et 
que, donc, le modéle classique d'excitatrice inspiré des 
machines à collecteur, n'est pas recommendé. 

Pour le moteur à induction alimenté par onduleur de 
courant, nous avons obtenu certaines conditions suffi- 
santes de stabilitê, analytiquement simples, qui recouvrent 
l'essentiel des phénomêénes et qui donnent une idée de 
influence des différents paramétres du systême. 


REFERENCES 


[1] —M. S. Garrido: Equations de la génératrice à 
redresseur. Etude du régime permanent, Re- 
vue E, vol. Vil, n.º 3, 1972. 

Special effects in rectified output A. EC. gene- 
rators, Int. Conference on Electrical Machines, 
Bruxelles, Sept. 1978. 

[24 — F. Labrique: Commutation processes in poly- 
phase commutatorless DC generators, Int. Conf. 
on Electrical Machines, Bruxelles, Sept. 1978. 

[3] — H. Buyse & M. S. Garrido: 4 new dynamical 
model of rectified output AC generators, Elec- 
tric Machines and Electromechanics, vol, 3, 
nº 1, 1979, 

[4] — L. Baes, H. Buyse & L. Pierrat: Etablissement 
de la fonction de transfert de la génératrice à 
redresseur, Bulletin de D.E.R. de I'EdF, sé- 
rie B, n.º 2, 1976. 

[5] —H. Buyse & L. Pierrat: Limit of intrinsic stabi- 
lity of a synchronous machine supplying an 
inductive load through a rectified, Int. Conf. 
Electrical Machines, Wien, Sep. 1976. 

[6] — Hanson & Thiry: Modélisation mathématique 
de Iexcitatrice sans balais des machines syn- 
chrones, TFE, U. C. Louvain, 1978, 

[7] — H. Buyse: Contribution à "étude de la machine 
à commutation electronique, Thêse de Doctorat, 
U, C. Louvain, 1976. 

[B] —L. Baes & M. S. Garrido: Compléments à 
Pátude du redressements triphasé classique, 
Revue E, vol. VII, n.º 1, 1975. 

[9] — E. Comell & T. Lipo: Modeling & design of 
controlled current induction motor drive systê- 
me, lJEEE Trans. Ind, Ap. vol. 14-13, n.º 4, 
July /August, 1977. 


TECNICA 454 


[10] — N. Sato, M. Ishida & N. Sawaki: Some charac- 
teristics of induction motor driven by current- 
-source inverter, IEEE/IAS Semiconductor 
Power Converter Conference, 1977. 

[11] — N. Sawaki & N. Sato: Steady-state and stability 
analysis of induction motor driven by current 
— source inverter, IEEE Trans. Ind. Ap, 
vol, 1A-13, n.º 3, May/June, 1977. 


REVISTA DE ENGENHARIA 


associação dos estudantes do 


[12] — J. Santana: Etude de la stabilité du moteur à 
induction alimenté par onduleur de courant, 
mémoire de Maitrise, U. C. Louvain, 1979, 

[13] — D. C. White & H. H. Woodson: Electro- 
mechanical energy conversion, Ed. John Willey, 
1959, 

[14] —- E. Gudefin: Les machines à commutation elec- 
tronique R. G. E., nº 1 Janv. 1968. 


instituto superior técnico 


Av. Rovisco Pais — 1000 Lisboa — Portugal — Telefone 88 93 23 — Telegramas AEIST 


TABELA DE PREÇOS DE ANÚNCIOS 
ADVERTISEMENT RATES 


1/4 página 


1 página 1/2 página 

page page page 
2.º ou 3.º capas 
ând & 3rd cover 8000800 6400800 4000800 
4.º capa 
th cover 9600$00 — — 
Última página antes do texto 
Last page before text 6400800 —. — 
Páginas intercaladas no texto 
Pages in text 5600$00 4000$00 3200800 
Qualquer outra página 
Any other page 4800800 3200800 2400800 


Publicidade redigida — Agravamento de 100 % 
Editorial advertising — 100 % more 


Encartes — Por cada folha tamanho A4 será cobrado 
o preço de 1 página. 

Inserts — For each sheet A4 (297 x 210 mm) the price 
of one page. 


Descontos: 5 números 20 % 
Dicounts: 10 insertions 40 % 


Ênúncios a cores: mais 


Colour advertisements: plus 1500800 


por cada cor 
per colour 


Os anúncios são pagos depois da publicação. 
Advertisement are paid after publication. 


1. Tiragem: 
Circulation: 


3000 


2. Linguas utilizadas na revista: Português, Inglês e 
Francês 


Languages used in our journal: Portuguese, English 
and French 


3. Tipo de leitores: Engenheiros e estudantes de enge- 
nharia 
Type of readershipe: Engineers and engineering 
students 


TÉCNICA 454 


4. Periodicidade: Mensal (excepto Agosto e Setembro) 
Frequency: Monthly (except August and September) 


5. Tipo de impressão: Tipografia. Anúncios em offset 
Screen: 120. Printing process: Typography. Adver- 
tisements in offset 


6. Mancha: 
Type area: 16,5xX25 em 
7. MN. de colunas: 
No. of columns: 2 
Altura da coluna: 


Column depth: 25 em 
Largura da coluna: 
Column width: 8 em 
B. Material de impressão usado: Qualquer 
Press maetrial required: Any 
9. Número corrente de páginas: 
Usual number of pages: 60 
10. Cores possiveis: Quaisquer 
Colours available: Any 
11. Comissão de agência: 
Agent comission: 20 % 


131 


HIDRAULICA 
G E ij Ã E Pelo Professor 


Alberto Abecassis Manzanares 


1.º VOLUME 


— FUNDAMENTOS TEÓRICOS Preço 675$00 


2.º VOLUME 


— ESCOAMENTOS LÍQUIDOS 


(A SAIR BREVEMENTE) PEDIDOS A TECNICA 


TUDOR 


Alimentadores de C.C. para Serviço Estacionário 


A TUDOR projecta e fabrica alimentadores e carrega 
dores de bateriais para todos os fins obedecendo as 
mais exigentes especificações. 


Características Eléctricas 


Carregador Principal 


“Tensão de alimentação ..3x380V T 10% 50Hz +5j | 


Potência Nominal. . .. .&3 KW 


* Funcionamento . .... .Regime flutuante estabili- 


Carregador Adicional 
* Potência Nominal. .... 


* Funcionamento . ..... Regime Tampão 


Bateria 


o mm 


Ha * Alcalina de Cádmio Níquel 


BLOCO para alimentação de permanência e Socorro dos equipamentos de Protecção e 
Comando dos Serviços Auxiliares das subestações. 
Fabricado para a EDP e E.P. 


| SOCIEDADE PORTUGUESA DO ACUMULADOR TUDOR, SARL. R. Actor Tasso 1 — Tel. 56 21 


TECNICA Mill 


—— — — — -— — a —— — -—. — 0 + DD —- 0 


Resumo dos Artigos publicados na «Técnica» n.º 454 


GARRIDO, M. S. — Algumas considerações sobre 
a electromecânica quasi-estacionária dos meios 
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O artigo mostra que, em electromecânica quasi- 
-estacionária, é preferível adoptar, logo à partida, a 
realtividade galileana, A razão é que uma teoria mais 
geral a partir da qual se pudesse, de maneira consis- 
tente, deduzir a electromecânica quasi-estacionária, só 
poderá existir a um nível demasiado complexo para as 
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formulação da mecânica dos fluidos incompressíveis, 
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polifásico com rectificação. Parte | — Análise do 
gerador com estrutura sinusoidal, «Técnica», Lis- 
boa, LIV (454) Nov. 1979, p. 71-81. 


Este artigo é a primeira parte de um estudo de 
comutação dum gerador polifásico sem colector, A par- 
tir do estudo de dois casos limite, o «sinusoidal» apre- 
sentado nesta primeira parte e o «triangular» apresen- 
tado na segunda parte, mostra-se que a forma das 
indutâncias mútuas, em função da posição, pode afectar 
consideravelmente o processo de comutação e, por con- 
seguinte, a característica externa do gerador. 

á ri experimentais confirmam as previsões 
teóricas, 


LABRIQUE, F.— Estudo da comutação num gerador 
polifásico com rectificação. Parte Il — Análise do 
gerador com estrutura triangular, «Técnica», Lisboa, 
LIV (454) Nov. 1979, p. 83-95. 


Este artigo é a segunda parte dum estudo sobre a 
comutação dum gerador polifásico sem colector. A par- 
tir de dois casos limite, o «triangular» aqui apresentado 
e o «sinusoidal» já apresentado na primeira parte, 
mostra-se que a lei de variação das inductâncias mú- 
tuas das espiras em função da posição pode afectar 
consideravelmente o processo de comutação e, por con- 
seguinte, a característica externa do gerador. 

A análise teórica é confirmada por resultados expe- 
rimentais obtidos com máquinas de laboratório. 
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LOURTIE, Isabel; CARVALHAL, Francisco e LOUR- 
TIE, Pedro — Projecto apoiado em computador de 
sistemas de controlo, «Técnica», Lisboa, LIV (454) 
Nov. 1979, p. 97-109. 


O artigo descreve detalhadamente o bloco de mani- 
pulação e tratamento de dados de um conjunto arti- 
culado de programas para o projecto apoiado em com- 
putador de sistemas de controlo. Apresentam-se os 
algoritmos utilizados e os correspondentes fluxogramas 
para além de alguns exemplos que os ilustram. 


SANTANA, J. — Modelação e estabilidade da asso- 
ciação ondulador de corrente-máquina assincrona, 
«Técnica», Lisboa, LIV (454) Nov. 1979, p. 111-118 


A associação ondulador de corrente-máquina assin- 
crona constitui um sistema para o qual as equações dos 
circuitos de contactos deslizantes são de aplicação 
imediata. 

A complexidade do estudo da estabilidade do sis- 
tema leva-nos a introduzir hipóteses de simplificação 
que o tornam viável analiticamente. Desta maneira é 
possível considerar dois tipos de instabilidades: eléc- 
tricas e mecânicas. 

A modificação dos parâmetros para assegurar a 
estabilidade de associação ondulador de corrente-má- 
quina assíncrona não pode ser vista como uma solução, 
porque assim se é conduzido a sistemas de baixo ren- 
dimento. 

Com a introdução de cadeias de retroacção obtém-se 
sistemas estáveis sem modificar as boas performances 
da máquina assíncrona, Os controlos da corrente con- 
tinua de alimentação e da frequência de escorrega- 
mento, com efeitos complementares, permitem eliminar 
os tipos de instabilidade considerados. 


GARRIDO, M. S. — Estudo analítico e aplicação 
das equações variacionais de geradores AC com 
saída rectificada, «Técnica», Lisboa, LIV (454) 
Nov. 1979, Pp. 119-131. 
EsEHa) 
São introduzidas as principais constantes de tempo 
que intervêm no funcionamento variacional dos gera- 
dores de rectificação e justificados vários níveis possi- 
veis de simplificação. Os resultados obtidos são apli- 
cados ao estudo da Tunção de transferência das 
excitatrizes sem escovas e ao estudo da estabilidade 


do motor de indução alimentado por ondulador de 
corrente, 


TÉCNICA 454 


Synopsis of articles published in «Técnica» n.º 454 
LIV — November 1979 


LOURTIE, Isabel; CARVALHAL, Francisco e LOUR- 
TIE Pedro — Control systems design computer 
aided, «Técnica», Lisboa, LIV (454) Nov. 1979, 


p. 97-109. 


The paper describes with detail the data manipula- 
tion and data preparation block of a control systems 
computer-alded design package, The algoriitms and 
corresponding flow charts are presented together with 
some application examples. 


SANTANA, J. — Modeling and stability of induc- 
tion motor driven by current source inverter, «Téc- 
nica», Lisboa, LIV (454) Nov. 1979, p. 111-118. 


This paper presents a dynamic model for the induc- 
tton motor driven by current source inverter, which is 
developped as an application of the commutation 
circuits theory. 

An analytical study of the system stability is possi- 
ble with the introduction of simplifying hypothesis in 
the mathematical model. We will consider two types 
of instabilities: electrical and mechanical. 

IH is shown that the actions on the parameters to 
obtain the stability induce to systems with poor effi- 
ciency. The solution is the introduction of the feed- 
“back loops. 

The controls of the direct current and of the slip 
frequency with complementary effects leave out the 
two types of instabilities: the first one is equivalent 
to the increase of the stator resistance and so It is 
able to eliminate the electrical instability, the second 
one eliminates the steady instability bringing an elec- 
tromechanical decrescent characteristics in the system, 


GARRIDO, M. S. — Analytical study and applica- 
tion of the small signal equations of the rectitied 
output AC generators, «Técnica», Lisboa, LIV (454) 
Nov. 1979, p. 119-131. 


This paper introduces the main time constants that 
influence the small signal behaviour of the rectified 
output AU generators and presents some reasonable 
simplifications. The results so obtained are used to 
study the transfert function of the brushless exciter 
and to study the stability of the induction motor 
driven by current source inverter. 
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GARRIDO, M. S. — Some remarks on the quasi- 
-static electromechanics of continuous media, «Téc- 
nica», Lisboa, LIV (454) Nov. 1979, p. 59-70. 


It iz show that, in formulating quasistatic eleciro- 
mechanics, galilean relativity is preferably used from 
the start. This is because a more general theory from 
which electromechanics would have to be consistently 
derived could only exist at a level of excessive comple- 
xity in comparison with the requirements of appli- 
cations. 

This conclusion is reached by setting up a new 
formulation, similar to the electromagnetic theory, for 
incompressible fluid mechanies. This allows a parallel 
treatment of both theories. 


LABRIQUE, F.— Study of the commutation pro- 
cesses in a polyphase rectified output AC gene- 
rator. Part on — Analysis of the generator with 
sinusoidal coupling, «Técnica», Lisboa, LIV (454) 
Nov. 1979, p. 71-81. 


This paper is the first part of a study on the 
commutation processes in a polyphase commutatorless 
D. C. generator. From lhe study of two limit cases, lhe 
«sinusoidal»s one presented here and the «triangular» 
one presented in the second part, it will be shown 
that the form of the variation of the mutual induc- 
tances as a function of the relative positions of the 
different coils can affect considerably the commutation 
processes and hence the generator output characteristic. 

Experimental results on laboratory test machines 
confirm the theoretical analysis. 


LABRIQUE, F.— Study of the commutation pro- 
cesses in a polyphase rectified output AC generator. 
Part two — Analysis of the generator with triangular 
coupling, «Técnica», Lisboa, LIV (454) Nov. 1979, 
p. 83-95. 


This paper is the second part of a study on the 
commutation processes in a polyphase commutatorless 
DC generator. From the study of two limit cases, the 
elriangulars one presented here and the «sinusoidal» 
ome, vet presented in the first part, it is shown that 
the form of the variation of the mutual inductances as 
a function of the relative positions of the different 
coils can affect considerably the commutation processes 
and hence the generator output characteristic. 

Experimental results on laboratory test machines 
confirm the theoretical analysis. 
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da Aguieira e sob licença da 
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o primeiro Alternador Motor tendo como 
principais características: 


Potencia 100 MVA . 

Velocidade de rotação 125 r.p.m. 

Tensão nominal 12 000 volts 
— | Diâmetro exterior 11,5 m 


Diâmetro da parte girante 95 m 


Peso do rotor -— 225 tonelados |) 
Peso total 420 toneladas já 
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